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This  report.  Volume  IV  of  the  series  on  "Modeling  Air  Traffic 
Performance  Measures,"  represents  the  culmination  of  3 years  invest- 
ment on  the  part  of  the  National  Aviation  Facilities  Experimental  Center, 
Federal  Aviation  Administration  and  the  Transportation  Program,  Depart- 
ment of  Civil  Engineering  at  Princeton  University.  The  initial  challenge 
was  to  take  air  traffic  communications  as  recorded  on  voice  tape  over  the 
New  York  Center,  and  using  this  raw  material  to  synthesize  a statistical- 
mathematical  simulation  model  that  could  be  put  to  practical  use. 

Volume  I of  this  series  contains  extensive  dictionaries  and 
analyses  of  message  elements;  Volume  II,  the  initial  statistical  analyses 
and  early  simulation  models;  Volume  III,  the  construction,  validation  and 
initial  exercises  of  the  appropriate  simulation  model. 

This  volume  reviews  the  simulation  model  for  the  New  York  Air 
Traffic  Control  Communications.  The  validation  of  the  model  and  the 
sensitivity  of  the  model  to  various  changes  in  the  input  variables  is 
described.  The  successful  use  of  the  model  for  the  Houston  Air  Traffic 
Control  Communications  system  is  reviewed.  Discussions  of  applications 
of  the  model  are  provided: 

(i)  for  estimating  communication  capacities  For  sector  types  as 

a function  of  the  aircraft  arrival  Intensity,  with  particular 
emphasis  on  determining  when  an  unstable  performance  appears, 

(ii)  for  determining  the  effects  of  reducing  the  number  of  trans- 
missions per  transaction, 

(ill)  for  measuring  consequences  resulting  from  tone  bursts  of 
different  lengths. 

The  model  has  been  used  to  characterize  "general"  sector  functions  (as 
contrasted  with  individual  sectors).  Tt  is  anticipated  the  model  may  be 
appropriate  for  simulating  oceanic  sector  functions.  Communication  queuing, 
a response  not  readily  measured  in  the  real  world,  has  been  Investigated. 

The  many-response  character  of  the  communication  system  has  been  simplified 
by  the  construction  of  indices.  Measures  of  non-stationary  performance  arc 
provided  along  with  estimates  of  the  associated  probabilities  of  occurrence. 
Initial  efforts  to  combine  the  sector  functions  into  networks  are  described 
Those  working  on  the  project  believe  the  present  air  traffic  control  simu- 
lation model  can  prove  of  immediate  use,  and  look  forward  to  its  employment. 

In  the  struggle  to  construct  a model  descriptive  of  real  world  data, 
many  novel  theoretical  questions  arose.  The  consequence,  a happy  synergism 
from  the  point  of  view  of  the  university  researchers,  lias  been  the  introduc- 
tion of  several  novel  statistical  tools.  Of  special  mention  is  the 


ill 


1 


application  of  the  statistic  G,  developed  by  Dr.  Der-Ann  Hsu,  used  to  test 
whether  the  structure  of  two  time  series  is  identical.  The  use  of 
confidence  regions  to  provide  probability  statements  for  the  onset  of 
nonstatlonarlty  in  system  performance  is,  we  believe,  new.  Both  these 
novel  applications  of  statistics  should  prove  widely  useful.  The 
study  of  sector  functions  has  led  naturally  to  the  consideration  of  the  ATC 
Communications  as  a network,  and  Markov  chain  models  have  been  constructed 

by  Mr.  Neil  W.  Polhemus  to  characterize  the  aircraft  flows  in  the 
New  York  ATC  system.  This  work  has  stimulated  additional  theoretical 
investigations  of  considerable  promise  associated  with  non-homogeneous 
point  process.  The  subject  is  of  great  concern  to  many  transportation 
scientists. 

This  project  has  been  a happy  coupling  of  practical  problems, 
real  world  data,  and  both  the  adaptation  and  development  of  theory. 

We  give  special  thanks  to  Mr.  Allen  Busch  of  the  FAA  for  providing 
cordial,  and  frank,  appraisal  and  direction  to  these  efforts. 


J.  Stuart  Hunter 
Principal  Investigator 
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CHAPTER  1 


SIMULATION  MODEL  FOR  AIR  TRAFFIC  CONTROL  COMMUNICATIONS 


1.1  Introduction 


In  the  continuing  effort  to  construct  a useful  model  for  Air  Traffic 
Control  (ATC)  Communications,  the  results  of  data  analyses  and  the  con- 
struction, validation  and  applications  of  a simulation  model  have  been  des- 
cribed in  detail  in  the  reports: 


Modeling  Air  Traffic  Performance  Measures 

Volume  I:  Message  Element  Analyses  and  Dictionaries 
(Report  No.  FAA-RD-73-147 , I) 

Volume  II:  Initial  Data  Analyses  and  Simulations 
(Report  No.  FAA-RD-73-147, II) 

Simulation  Model  for  New  York  Air  Traffic  Control  Communlca t ion s 
(Report  No.  FAA-RD-74-203) 

Data  used  in  these  reports  for  modeling  ATC  communications  were 
collected  from  the  New  York  area,  during  a 2 --hour  peak  traffic  period, 
on  April  30,  1969,  Message  contents,  origins,  frequences  and  time  durations 
of  voice  communications  transmissions  were  analyzed  and  reported  in  full 
detail  in  volume  I.  Extensive  analyses  of  communication  data,  on  each  of 
the  elements  comprising  the  ATC  performance,  were  described  in  volume  II. 

In  the  third  volume  the  structure  of  a simulation  model,  based  on 
the  New  York  data  file,  was  described  and  its  generalization  presented 
through  the  use  of  various  master  equations  and  tables.  Validation  of  the 
simulation  model  was  performed  by  comparing  generated  responses  with  their 
historical  counterparts  through  the  use,  in  part, of  newly  developed  sta- 
tistical techniques.  The  results  of  comparisons  indicated  that  the  simu- 
lation model  did  operate  satisfactorily  not  only  for  individual  sectors 
but  for  sector  types  as  well.  Applications  of  the  simulation  model  were 
also  reported  in  the  third  volume.  One  special  example  concerns  the  deter- 
mination of  communications  capacity  for  HI  sectors,  a study  of  interest  to 
those  Investigating  oceanic  air  traffic.  The  value  of  the  simulation 
model  was  especially  demonstrated  in  a study  of  communications  queuing, 
that  is,  the  frequency  and  length  of  time  a pilot  must  wait  before  he  can 
use  the  communication  channel.  There  is  no  direct  measure  of  communications 
queuing  available  in  the  historical  data.  Other  related  studies  were  also 
presented  on  the  subjects:  (1)  response  coupling;  (il)  air  traffic  control 
network  capacities;  (ill)  sequential  structure  of  the  intercommunication 
gap  lengths;  and  (Iv)  tlie  applications  of  Sichel  distribution  for  number  of 
transmissions  (TR)  per  communications  transaction  (CT). 


1-1 


Although  the  studies  outlined  above  are  available,  to  prepare 
readers  for  a further  discussion  in  this  volume,  we  now  briefly  review 
some  terminology  and  concepts. 


1.2  Elements  in  the  Simulation  of  Air  Traffic  Control  Communications 

It  is  appropriate  to  describe  briefly  the  characteristics  of  the 
ATC  communications  system.  The  air  space  within  an  airport  region  is 
divided  into  a number  of  "sectors,"  each  with  its  own  controller.  For 
Instance,  the  New  York  area  contains  101  separate  sectors.  These 
sectors  are  divided  into  12  different  sector  functions,  namely.  Low 
Altitude  Transitional  (LT) , Low  Altitude  Enroute  (LE) , High  Altitude  En- 
route  (HI),  Clearance  Delivery  (CD),  Ground  Control  (GN),  Local  Control 
(LC) , Local  and  Ground  Control  (LG),  Approach  Control,  Non-Radar  (AP) , 
Radar  Departure  Control  (DP) , Radar  Arrival  Control  (AR) , Radar  Arrival 
and  Departure  Control  (AD),  and  Radar  Advisory  Position  (RA) . The  sector 
functions  are  partitioned  into  three  distinct  groups:  enroute  sectors 
(LT,  LE,  HI),  tower  cab  sectors  (CD,  GN,  LC,  LG,  AP)  and  IFR  room  sectors 
(DP,  AR,  AD,  RA). 

Within  a sector,  and  using  only  a single  communication  channel, 
a number  of  communications  are  made  between  the  controller  and  each  air- 
craft pilot  within  that  sector's  jurisdiction.  A completed  conversation 
between  an  aircraft  in  a sector  and  the  controller  is  termed  a "communi- 
cations transaction"  or  "CT."  Each  CT  is  composed  of  separate  "trans- 
missions" or  "TR’s"  which  are  made  alternatively  by  pilot  and  controller. 
Each  TR  in  turn  contains  one  or  more  "message  elements"  or  "ME's'.  " 
volume  I contains  extensive  dictionaries,  listing  the  ME's  comprising 
individual  TR's  along  with  statistics  characterizing  their  time  duration 
and  frequency. 


1.2.1  Simulation  Input  Variables 

The  first  input  variable  required  by  the  simulation  model  is  the 
aircraft  arrival  stream  entering  a sector.  It  was  found  that  a Poisson 
process  described  the  pattern  of  aircraft  arrival  to  a sector  adequately 
for  most  sectors  in  the  New  York  area.  A Poisson  process  is  the  one  in 
which  the  number  of  occurrences  of  an  event  (e.g.  the  aircraft  arrival) 
in  a unit  of  time  can  be  represented  by  the  probability  model 

p(j)  = e®  . /j  ! (1.1) 

where  p(j)  is  the  probability  of  j occurrences,  S is  the  arrival  rate 
and  j=0,  1,  2,  ...  is  the  number  of  aircraft  arrivals  obseryed  in  each 
time  unit. 


For  each  of  the  aircraft  assigned  to  the  control  sector,  the 
number  of  communications  transactions  (CT's)  associated  with  this  air- 
craft becomes  the  second  input  variable.  From  the  analyses  presented  in 
Volume  II,  chapter  4,  a negative  binomial  distribution  was  shown  adequate 


in  describing 
is 

P(r) 


this  variable 


k+r-1  k 

r 

k-1  ^ 


, In  mathematical 
(1-p)^,  r=0,l,2,.. 


form,  this  distribution 
(0<p<l),  (1.2) 


where  p and  k are  parameters  and  p(r)  is  the  probability  of  r+m  CT's 
per  aircraft  where  m=2  for  the  enroute  sectors  and  m=l  for  others.  A 
master  table  for  the  values  of  p and  k estimated  from  the  New  York  data 
for  each  of  the  sector  types  is  given  in  table  1.1.  Table  1.1  also 
indicates  the  test  statistic  for  a lack  of  fit  test.  The  significant 
lack  of  fit  observed  on  the  LC  and  AR  sectors,  indicating  the  inappropri- 
ateness of  the  negative  binomial  model,  is  due  to  the  nonhomogeneous 
performance  of  individual  sectors  within  the  function.  Further  dis- 
cussion on  this  issue  will  be  given  later  in  this  volume. 


The  third  variable  requiring  specification  is  the  time  length 
between  two  consecutive  CT's  for  the  same  aircraft,  termed  intercommuni- 
cation gap  lengths.  The  actual  procedure  for  generating  this  variable 
for  the  simulations  is  complicated.  It  is  determined  using  two  param- 
eters 3 and  relating  the  intercommunication  gap  length  to  the 
number  o?  CT's  per  aircraft,  A master  table  for  the  values  of  a and 
a^  is  given  in  table  1.2  for  each  of  the  sector  types. 


An  Important  input  variable  for  a simulation  is  the  length  of 
CT.  This  variable  is  a function  of  the  number  of  TR's  per  CT  and  trans- 
mission length  (TR  length).  Unfortunately,  the  number  of  TR  per  CT  is 
too  irregular  to  be  fitted  by  any  simple  mathematical  model,  and  there- 
fore histograms  compiled  from  historical  data  are  used  to  provide  this 
variable  in  the  simulations.  For  the  second  variable,  TR  length,  the 
gamma  distribution  proved  adequate  for  representation.  This  distribution 
is 


f(t) 


1 a-1  -t/> 

t e 

X“r(ct) 


0<t<«> 

u>0  , X>0 


(1.3) 


where  t is  the  TR  length,  and  m and  X are  parameters.  In  order  to  deter- 
mine the  values  of  nt  and  X for  each  of  the  sector  types,  two  simultaneous 
master  equations,  relating  the  values  of  a and  X to  the  aircraft  arrival 
rate,  were  compiled.  The  following  two  equations  apply  to  all  sectors 
specified:  (The  CD,  AP  and  RA  sector  functions  are  exceptions  due  to  their 
small  numbers  and  unique  behavior). 


(i) 

n • 

X 

= 3. 

,70  - 0.0158  e 

(General) 

= 2. 

.97  - 0.0042  6 

(LC) 

(ID 

(n  - 

1) 

X = 

= 2.0 

(General) 

(-1  - 

1) 

A = 

= 1.7 

(ON) 

(«  - 

1. 

13) 

X = 1.72 

(DP) 

(a  - 

1. 

5) 

X = 1.0 

(AD) 
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Table 

1.1 

Fit  of  Negative  Binomial  Distribution 
to  Number  of  Communications  Transactions 
per  Aircraft 

Sec  tor 
Function 

No. 

Obs. 

t 

m 

_E 

_k 

Degrees  (jf 
Freedom 

.L.. 

LT 

514 

2 

0.41 

2.36 

14 

21.07 

LE 

498 

2 

0.46 

2.30 

12 

13.89 

HI 

379 

2 

0.58 

4.03 

11 

4.14 

CD 

200 

I 

0.94 

15.87 

4 

8.30 

GN 

905 

1 

0.41 

0.88 

10 

16.15 

LC 

1075 

1 

0.71 

6.83 

11 

125.49* 

LG 

128 

1 

0.69 

6.83 

8 

10.78 

AP 

49 

1 

0.63 

2.63 

5 

7.02 

DP 

252 

1 

0.88 

30  45 

10 

7.27 

AD 

88 

1 

0.58 

4.85 

9 

8.28 

AR 

225 

1 

0.43 

4.54 

16 

33.46* 

RA 

95 

1 

0.46 

3.01 

10 

7.72 

Significant  lack  of  fit  at  the  0.01  level. 


To  be  comparable  to  a variable  from  the  negative  binomial  distribution, 
subtract  m from  the  actual  number  of  CT/AC. 
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Table  1.2  Average  Values  of  Regression  Estimates  a_,  a 
for  Cap  Lengths  in  Each  Sector  Function 


Sector 

Function 

o 
1 eg 

^1 

Sector 

Function 

io_^ 

LT 

(12)* 

4.903 

-0.109 

LG  (5) 

4.919 

-0.028 

LE 

(11) 

5.579 

-0.164 

AP  (4) 

4.671 

0.069 

HI 

( 9) 

5.184 

-0.083 

DP  (6) 
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-0.024 

IIA  (3) 

4.902 

-0.048 

Value  in  the  parenthesis  indicates  the  number  of  sectors  in  the 
function  specified. 
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Once  again  0 is  ttie  expected  number  of  aircraft  arrivals  per  hour.  A 
pair  of  parameter  values  a and  A can  be  obtained,  given  0,  by  solving 
thest!  two  equations.  Altogether,  seven  Input  parameters,  0,  p,  k,  a , 
a , u and  A are  required  to  Initiate  a simulation,  plus  the  appropriate 
histogram  for  the  number  of  TR/CT. 


1.2.2  Simulation  Response  Variables 

In  addition  to  these  Input  variables,  three  major  responses  are 
generated  from  each  of  the  simulation  exercises.  They  are  (1)  channel 
utilization;  (li)  aircraft  loading;  and  (111)  communications  queuing. 
Channel  utilization  (CU)  is  defined  as  the  proportion  of  time  the  single 
communication  channel  is  utilized.  Usually  In  our  analysis,  the  channel 
utilization  is  computed  for  each  60-second  interval.  The  time  series 
thus  compiled  is  then  used  for  further  analysis.  The  second  response, 
aircraft  loading,  is  defined  as  the  number  of  aircraft  in  a sector  at  a 
given  time  point.  As  in  the  case  of  channel  utilization,  this  response 
is  also  averaged  for  each  60-seconds  to  form  a time  series.  The 
third  response,  not  available  in  the  historical  records,  is  of  special 
interest.  There  are  two  aspects  to  this  response.  One  is  the  number  of 
aircraft  waiting  to  converse  with  the  controller  (or  conversely,  the 
number  of  aircraft  the  controller  has  to  converse  with  but  cannot  because 
he  is  occupied  by  other  aircraft)  at  a given  time  point.  The  other  aspect 
is  the  time  length  an  aircraft  waits  to  converse  with  the  controller  (the 
converse  also  applies)  for  each  CT  attempted.  These  two  aspects,  the 
number  of  waiting  aircraft  and  the  queue  time,  provide  information 
supplementary  to  the  other  in  the  analysis  of  queuing  status.  In  the 
simulation  60-second  averages  of  the  number  of  aircraft  waiting  to 
converse  are  compiled  as  a time  series. 


1.3  Validat ion  of  the  Simulation  Model 

Validation  of  the  simulation  model  was  performed  by  comparing 
the  generated  responses,  namely,  the  aircraft  loading  and  the  channel 
utilization,  compiled  as  time  series,  with  their  corresponding  historical 
time  series,  A test  procedure  designed  to  detect  the  difference  of  two 
time  series  with  respect  to  their  structural  parameters  (the  mean,  auto- 
correlation and  error  variance)  was  employed  to  do  the  comparison  more 
rigorously.  Two  test  statistics  were  used  in  this  comparison  procedure. 

One  is  the  test  statistic  G(iO,',)  distributed  as  a variable  and  useful 
for  comparing  two  time  series  witli  respect  to  their  estimated  auto- 
regressive parameters  and  error  variances.  The  other  is  a test  statistic 
distributed  as  Student’s  t appropriate  for  testing  the  difference  in  means 
of  two  autocorrelated  time  series. 

Results  of  validation  of  the  simulation  exercise  on  the  indivldiuil 
sector  basis  (the  seven  initiating  parameters  given  are  unique  to  the 
sector  being  simulated),  using  typical  sectors  from  each  of  the  eight 
types  of  control  function  considered,  were  sliown  in  cliapter  3 of  volume  ij; 
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for  the  responses:  aircraft  loading  series  and  the  channel  utili- 
zation series  separately.  The  resulting  statistics  indicate  that  the 
simulation  model  successfully  reconstructed  the  dynamic  responses  of 
concern.  Another  simulation  exercise  on  the  HI  sector  type,  in  which 
all  the  seven  initiating  parameters  were  obtained  ^rom  master  tables 
and  equations  descriptive  of  the  functions,  was  also  examined  and  shown 
satisfactory.  A continuing  study  on  the  validation  of  sector  type 
simulations  for  other  sector  functions  is  presented  in  chapter  2 of 
this  volume. 


1.4  An  Application  to  the  Estimation  of  Sector  Communication  Capacity 

The  simulation  model  was  exercised  for  sector  475  HI,  a repre- 
sentative enroute  sector,  in  an  experiment  in  changing  traffic  density. 

The  aircraft  arrival  rate  per  hour  was  increased  gradually  from  a level 
equal  to  the  historical  record  to,  finally,  a level  constrained  by  the 
celling  of  communication  channel  capacity.  Graphs  of  the  generated 
response  time  series,  the  aircraft  loading  series,  the  channel  utilization 
series  and,  most  important,  the  queue  lengths  series,  were  presented,  in 
volume  III,  chapter  4,  for  each  of  the  traffic  densities  tried.  The 
onset  of  the  operation  instability,  or  an  explosive  status,  was  clearly 
seen  by  the  sudden  leap  in  the  level  of  the  queue  length  series.  The 
aircraft  loading  series  revealed  similar  behavior.  In  addition  to  this 
eye  Inspection  of  the  response  series,  a test  statistic  useful  for  de- 
tecting nonstationarity  of  a time  series  model  was  used  to  Indicate  the 
operational  Instability  (unstable  or  overloaded  performance)  observed 
from  the  aircraft  loading  series.  A continuing  discussion  on  the  simu- 
lation experiments  for  communications  capacities  will  be  given  in  chapter  3 
of  this  volume. 


1.5  Network  Analyses 

In  order  to  construct  a simulation  model  that  is  applicable  to 
a whole  airport  region,  it  is  necessary  to  examine  the  patterns  of  trans- 
fer of  aircraft  from  one  sector  to  another  in  addition  to  the  within 
sector  characteristics.  The  historical  movements  of  aircraft  between 
sectors  were  followed  and  the  numbers  of  aircraft  transferring  between 
all  pairs  of  sectors  determined.  These  observed  transitions  were  then 
summarized  by  sector  function  and  are  shown  in  a transition  matrix  in 
table  1.3. 


More  detailed  information  on  the  patterns  of  movement  of  aircraft 
between  separate  sectors  was  presented  in  a number  of  schematic  diagrams 
obtained  from  further  analyses  of  the  data.  Figure  1.1  shows  an  example 
of  one  of  the  schematic  diagrams  of  transitions  between  various  sectors. 
This  work  is  an  extension  of  the  studies  on  the  networks  of  routes 
between  sectors  initially  reported  upon  in  chapter  7 of  volume  II. 
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An  important  network  problem  presented  by  the  study  of  sector 
transitions  is  the  determination  of  system  capacity,  given  the  structure 
of  individual  sectors.  An  initial  formulation  of  the  ATC  system  in  net- 
work terminology  was  discussed  in  chapter  5,  section  5.3,  of  the  third 
volume,  and  maximum  flow  algorithms  were  applied  to  a sample  network. 

In  the  present  volume  a more  complete  network  model  that  includes  the 
network  arrival  streams,  the  formulation  of  the  network  structure,  and 
the  construction  of  transition  patterns  is  presented  in  chapter  4. 

1.6  Previev/  of  the  Following  Chapters 


The  present  volume  continues  the  effort  presented  in  tlie  earlier 
three  volumes  listed  at  the  beginning  of  this  chapter,  and  summarized  in 
the  previous  sections.  One  objective  is  to  generalize  the  framework  of 
the  simulation  model  for  adaption  to  various  sector  types  and  geographical 
characteristics.  Further  validation  on  the  model  for  sector  type  simu- 
lations has  been  accomplished.  Applications  have  been  extended  to  studies 
on:  (i)  sector  type  communication  capacities;  (li)  effects  of  a tone- 
burst  scheme  for  aircraft  identification,  proposed  by  FCC;  and  finally 
(ill)  the  influence  of  hypothesized  changes  in  the  number  of  TR/CT.  An 
extensive  analysis  of  N.Y.  ARTCC  network  structure  — the  theoretical 
formulation,  empirical  modeling,  and  tests  — is  also  discussed  in  this 
volume  as  an  initial  step  toward  network  simulation  modeling.  This  net- 
work framework  can  hopefully  be  useful  to  the  study  of  ATC  petformance 
in  ocean  areas. 

In  chapter  2 of  this  volume  we  first  review  the  sector  type 
simulations  defined  in  the  third  volume  and  perform  the  simulation 
exercises  for  nine  out  of  the  12  sector  types  in  N.Y.  area.  Results 
indicate  that  the  model  can  work  adequately  for  all  of  these  sector  types 
given  a slight  modification  on  one  or  two  of  the  input  parameters,  for 
some  of  the  sector  types.  The  simulation  model  is  further  extended  and 
adapted  for  distinct  geographic  characteristics.  Using  the  Houston 
enroute  sector  groups  as  an  example,  a comparison  of  various  key  input 
parameters  for  New  York  versus  Houston  data  has  been  made  to  determine 
which  of  the  input  variables  should  be  adjusted.  No  basic  changes  in  the 
structure  of  the  simulation  model  proved  necessary.  Simulations  on 
Houston  sector  types  are  also  validated. 

Several  applications  of  the  sector  type  simulation  model  are 
demonstrated  and  described  in  chapter  3.  The  subjects  included  are: 

(i)  the  estimation  of  sector  type  communication  capacities  for  the  New 
York  sectors;  (il)  the  effects  of  a proposed  tone-burst  scheme  designed 
to  enforce  the  .air  traffic  control  performance;  and  (ili)  the  influence 
of  a change  in  the  expected  number  of  TR/CT.  The  first  subject  is  an 
extension  of  the  work  presented  in  chapter  4 of  the  third  volume,  while  the 
second  subject  is  an  interesting  realistic  exercise  of  the  simulation 
model.  In  the  latter  case,  a spline  fitting  method  is  used  to  generate  a 
complete  average  response  surface  for  the  major  variables  studied.  For 
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the  last  subject  mentioned,  an  exponential  weighting  function  Is  used  to 
"Improve"  the  Input  variable  available  as  historical  histograms. 

Chapter  4 provides  a more  complete  network  model  discussing  the 
mathematical  model  for  aircraft  arrival  streams  to  the  network  as  well  as 
the  second  order  Markov  model  for  the  transition  pattern  of  the  aircraft 
transferring  from  one  sector  to  another.  Extensive  network  data  analysis 
Is  given  to  examine  the  adequacy  of  the  postulated  model. 

Miscellaneous  related  studies  are  presented  In  chapter  5. 

Subjects  Included  are:  (1)  sensitivity  analysis  of  the  expected  channel 

utilization;  (11)  the  construction  of  operational  Indices  for  ATC  per- 
formance; (111)  transfer  function  models  for  aircraft  flows  at  adjacent 
sectors;  (iv)  modeling  long-term  time  series  for  the  Inputs  and  responses 
of  an  ATC  system;  and  (v)  queuing  time  forecasting  In  the  ATC  communication 
system. 

The  last  chapter  summarizes  the  conclusions  and  recommendations 
made  In  previous  chapters  and  sections. 
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CHAPTER  2 


VALIDATION  OF  A GENERAL  SIMULATION  MODEL  FOR  ATC  COMMUNICATIONS 


2.1  Introduction 

In  simulation  modeling,  we  note  that  a model,  useful  in  studying 
responses  for  a complicated  system,  can  work  appropriately  only  when  an 
adequate  set  of  input  variables  is  provided.  All  simulation  models  are 
constrained  to  be  simplifications  of  reality.  Thus  they  require  input 
rich  in  information  extracted  from  both  the  structure  and  history  of  the 
system.  In  the  ATC  simulation  model  constructed  to  date,  the  simulation 
system  starts  with  seven  Initiating  parameters  plus  a histogram  (for  the 
distribution  of  //TR/CT).  No  further  inputs  are  used,  or  required,  to 
exercise  the  present  model.  Subsequent  close  study  has  shown  that  the 
responses  derived  from  the  simulation  appear  to  depend  only  on  those 
components  which  can  be  specified  by  the  seven  input  parameters.  To 
have  Included  in  the  ATC  simulation  program  all  the  details  descriptive 
of  an  individual  ATC  sector  is  logically  untenable  and  computationally 
too  cumbersome. 

The  ATC  simulation  procedure  consists  of  two  stages:  (i)  structure 
and  data  analysis;  and  (ii)  verification  and  exercise  of  the  (GPSS)  com- 
puter program.  These  two  stages  are  equally  important.  The  ATC  structure 
and  data  analysis  have  been  described  in  full  detail  in  volumes  I and  II. 

The  computer  model  was  explained  and  illustrated  in  detail  in  the  third 
volume,  and  summarized  in  chapter  1 of  this  volume. 

In  these  previous  studies,  the  structure  of  the  simulation  model, 
and  the  corresponding  derivation  and  use  of  the  Initiating  parameters,  was 
established  and  verified  using  information  from  individual  sectors.  The 
next  step  ot  the  simulation  development  was  then  to  construct  a useful 
"general"  sector  representative  of  a function.  First  of  all,  this 
required  the  input  parameters  extracted  from  historical  data  of  all  the 
sectors  performing  a particular  function.  These  are  now  available  as 
"master  tables"  and  "master  equations."  From  the  discussion  given  in  the 
beginning  paragraph  of  this  chapter,  a "general"  simulation  is  therefore 
possible.  However,  validations  are  required  to  examine  the  usefulness  of 
the  "general"  model. 

The  validation  of  the  simulation  model  for  individual  sectors  was 
straightforward,  since  the  Individual  sector  data  was  in  hand.  However, 
the  validation  of  the  simulation  model  for  sector  function  encounters  the 
problem  of  non-homogeneity  of  performance  amongst  the  individual  sectors 
performing  that  function.  In  a few  cases,  responses  provided  by  the 
sector  function  model  did  not  agree  well  with  the  corresponding  responses 
for  the  individual  sectors.  After  much  study  it  was  discovered  that  modest 
changes  in  only  one  or  two  of  the  function  simulation  parameters  were 
required  to  bring  the  function  and  sector  responses  into  compliance.  Thus, 
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although  the  simulations  provided  by  the  function  parameters  may  not 
always  provide  response  time  series  sufficiently  similar  to  those  avail- 
able from  the  individual  sectors,  the  function  simulations  remain 
valuable.  Further  function  simulations  proved  very  useful  in  extension 
to  geographical  areas  other  than  New  York. 

In  section  2.2  of  this  chapter ^the  data  collected  from  Houston 
ARTCC  (consisting  exclusively  of  the  LT,  LE  and  HI  enroute  sectors)  are 
analyzed  and  compared  against  the  results  with  the  New  York  enroute 
sectors.  A great  similarity  between  the  regions  exists.  The  only 
exception  in  the  performance  of  these  two  distinct  geographical  areas 
concerns  the  variable  f»CT/AC.  A straight  extension  of  the  simulation 
model  to  the  Houston  sector  functions  is  made  using  entire  "master  kit" 
borrowed  from  N.Y.  ARTCC,  except  for  the  parameters  p and  k corresponding 
to  the  Houston  variable  #CT/AC.  Results  of  the  validation  of  this  per- 
formance presented  in  section  2.3  indicate  the  adequacy  of  the  general 
simulation  model  for  both  the  N.Y.  and  Houston  areas.  The  Inference 
that  the  simulation  model  is  appropriate  to  geographical  areas  other 
than  New  York  is  thus  sustained. 


2 o 2 Comparison  of  ATC  Communications  Performance  in  New  York  versus 

Houston  Enroute  Sector  Functions 


In  this  section  ATC  communications  data  collected  from  32  New 
York  enroute  sectors  and  23  Houston  enroute  sectors  are  compared.  Results 
of  the  analyses  indicate  that  the  only  major  difference  between  New  York 
and  Houston  ATC  communications  is  with  the  distribution  of  the  number  of 
communications  transactions  per  aircraft  (#CT/AC). 


2.2.1  The  Data  and  Their  Comparisons 

The  data  of  New  York  sectors  were  collected  and  described  in 
Volume  II  for  a 2-hour  period  (20/00/00  - 22/00/00  GMT  on  April  30, 

1969.  The  Houston  data  were  also  recorded  for  2-hour  periods,  and 
furthermore,  repeated  for  the  same  sector  at  different  times  in  a day 
and/or  the  same  time  at  different  dates.  The  sampling  periods  spread 
between  February  17,  1971, and  March  5,  1971.  Altogether  109  distinct 
sets  of  2-hour  data  are  available.  Detailed  description  of  the  data 
structure  is  given  in  table  2.1. 

Four  important  variables:  (1)  transmission  lengths;  (11)  trans- 
action lengths;  (ill)  Intertransaction  gap  lengths  and  (Iv)  number  of 
CT's  per  aircraft  are  analyzed  and  compared  for  Houston  versus  New  York 
enroute  sectors.  Gamma  distributions  were  fitted  to  the  transmission 
length  data  and  their  parameters  estimated.  The  scatter  diagrams  of 
these  parameter  estimates  are  shown  in  figures  2.1  and  2.2  for  New  York 
and  Houston  enroute  sectors,  respectively.  A striking  similarity  between 
these  two  diagrams  is  observed.  Their  locations,  dispersions  and  the 
relationships  between  the  estimates  n and  t for  both  regions  are  almost 
identical.  Similar  comparisons  for  transaction  lengths,  as  well  as  the 
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intertransaction  gap  lengths  are  given  In  figures  2.3  - 2.6,  The 
similarities  observed  are  again  remarkable. 

However,  the  distributions  of  the  number  of  communication  trans- 
actions per  aircraft  appear  to  be  different  for  the  LT  and  HI  functions  in 
these  two  distinct  geographical  areas. 

2.2.2  Analyses  of  the  Number  of  CT's  per  Aircraft 

Negative  binomial  distributions  were  fitted  to  the  Houston  data 
for  the  number  of  CT's  per  aircraft.  Details  of  the  parameter  estimates 
for  each  of  the  109  sets  of  records  are  displayed  in  table  2.2.  Results 
from  fitting  tlie  distribution  to  the  data  compiled  by  sector  functions 
are  shown  in  table  2.3.  Histograms  of  the  number  of  CT's  per  aircraft 
are  shown  for  the  Houston  data  in  contrast  to  those  of  the  New  York,  in 
■figure  2.7.  A further  list  of  the  estimates  of  p and  k and  the  expected 
values  is  given  in  table  2.4  for  both  the  Houston  and  New  York  data. 

In  these  tables  and  figures,  the  nonhomogeneous  behavior  observed 
from  the  Houston  HI  sectors  in  the  run  codes  194H  - 217H  versus  the 
remainder  (data  represented  bv  the  last  10  run  codes  were  collected  from 
two  controller  sectors  at  different  periods  of  time)  suggested  the  sep- 
aration of  the  Houston  HI  sector  function  into  two  subgroups.  The  shapes 
of  the  distribution  functions  and  the  estimates  of  p and  k,  particularly 
the  mean  value  of  the  distributions,  for  these  two  subgroups,  are 
apparently  distinct,  A careful  inspection  of  table  2.4  indicates  that 
except  for  the  LE  function,  the  p and  k values  for  N.Y.  sectors  are 
different  from  those  of  the  Houston  sectors.  A rigorous  statistical  test 
for  the  comparison  of  two  negative  binomial  variables  is  complicated  and 
not  presented  here.  For  the  purpose  in  hand  we  consider  the  distributions 
to  be  different  and  correspondingly  adjust  the  specification  of  parameter 
values  in  the  simulations  for  Houston  sector  functions  as  shown  in  the 
next  section. 


2.3  Validation  of  Simulations  for  New  York  and  Houston  Sector Functions 

A general  ATC  simulation  for  New  York  HI  sector  function  has  been 
demonstrated  in  chapter  3 of  volume  III,  In  the  simulation  mentioned, 
input  parameters  were  generated  from  various  master  tables  and  equations 
compiled  for  the  N.Y.  area  . To  further  explore  the  general  simulation 
for  other  sector  functions,  techniques  developed  for  the  comparison  of 
two  time  series,  with  respect  to  their  stochastic  structures,  are  again 
applied  to  detect  possible  differences  in  essential  structure  between  tiie 
generated  time  series  responses  and  their  historical  counterparts.  Since 
the  historical  responses  can  only  be  observed  for  individual  sectors,  not 
for  a sector  function  as  a whole,  comparisons  of  time  series  generated 
from  a function  simulation  versus  historical  time  series  of  individual 
sectors  within  the  function  were  carried  out  in  the  following  fashion.  A 
function  simulation  using  parameters  determined  by  master  tables  and 
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equations  was  first  executed  and  validated  against  responses  observed 
from  a single  representative  sector  within  a function.  Significant  dis- 
crepancies, if  any,  between  the  simulated  and  historical  responses,  were 
then  studied,  and  adjustment  made  on  the  value  of  the  input  function 
simulation  parameters.  Some  Individual  sectors,  exclusively  in  the 
terminal  sector  functions,  necessitated  adjustments  on  either  the 
#TR/CT  or  #CT/AC  input  variables.  With  these  occasional  adjustments, 
the  simulations  for  a sector  function  produced  responses  satisfactorily 
similar  to  the  historical  time  series  of  an  individual  sector.  Results 
of  the  final  validations  of  the  simulations  of  nine  sector  functions  of 
N.Y.  area  are  displayed  in  tables  2.5  and  2.6. 

The  master  tables  and  equations  compiled  from  the  New  York 
data  were  then  used  to  perform  function  simulations  for  the  Houston 
enroute  sector  functions.  The  only  changes  required  in  the  master  tables 
were  the  values  of  p and  k,  for  the  distributions  of  the  number  of  CT's 
per  aircraft.  The  new  values  of  p and  k were  estimated  from  the  Houston 
data.  The  simulation  model  gave  a surprisingly  good  performance  as  seen 
from  the  results  of  validation  tests  shown  in  tables  2.7  and  2.8.  This 
suggests  the  general  applicability  of  the  simulation  procedure  for 
sector  functions  in  areas  other  than  New  York, 

Complete  statistical  tables  and  figures  generated  from  the  New 
York  and  Houston  simulation  programs  for  the  validations  of  function 
simulations  were  compiled  along  with  the  historical  time  series  of  the 
individual  sectors  having  been  compared.  This  material  is  now  available 
in  a NAFEC-Prlnceton  University  Interim  Report  (PU-36). 


2.4  ^urana ry  and  Concluding  Remarks 

A previous  simulation  model  adequate  for  mimicking  the  dynamic 
behavior  of  an  individual  ATC  sector  has  been  generalized  and  applied 
to  the  function  simulations  for  both  New  York  and  Houston  areas. 

Studies  on  the  verifications  indicated  that  the  "general"  simulation 
model  is  rich  in  Information  and  convenient  for  use  in  extensions  to 
varied  studies.  Only  one  or  two  sensitive  variables,  mostly  likely  the 
distributions  of  #TR/CT  or  #CT/AC,  need  to  be  adjusted  in  order  to 
extend  the  simulations  to  other  geographical  areas.  Much  effort  is 
saved  by  concentrating  on  tuning  one  or  two  input  parameters,  instead 
of  reconstructing  a new  model,  for  the  simulations  of  ATC  communications 
in  different  geographical  areas,  as  well  as  for  distinct  types  of  simu- 
lation studies. 

Future  experience  with  more  extensive  data  collected  from  various 
environments,  such  as  the  ATC  systems  in  the  ocean  areas,  will  certainly 
help  develop  the  present  model  into  an  even  more  mature  form. 
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The  Structure  of  Houston  Data  (cont.) 
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LT : Low  altitude  transitional 


bSTimates  ot  the  Parameters  ot  Gamma  Distributions 
the  Tronsmission  Lengths  (New  York  Data) 


Figure  2.2  Estimates  of  the  Parameters  of  Gamma  Distributions  Fitted  to  the  Transmission  Lengths 
( Houston  Data ) 
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Figure  2.5  Mean  vs  Standard  Deviation  of  Inter  transaction  Gap 
(New  York  Data) 
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Table 

2.2  Results  of  Fitting  Shifted  Negative  Binomial  Distributions 
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(Houston  Individual  Sectors) 
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Table  2.2 


Pitting  Shifted  Negative  Binomial  Distributions  to 
he  Number  of  Communications  Transactions  per  Aircraft  (cont.) 
(Houston  Individual  Sectors) 
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HOU/N 

12 

0.60 

2.92 

0 

3 

2.36 

148H 

LT 

HOU/N 

17 

0.83 

13.78 

0 

5 

7.62 

149H 

LT 

HOU/N 

24 

0.56 

2.18 

0 

4 

3.61 

150H 

LT 

HOU/SW 

19 

0.38 

1.62 

0 

6 

6.78 

151H 

LT 

HOU/SW 

29 

0.93 

29.94 

0 

5 

2,02 

152H 

LT 

HOU/SW 

22 

0.65 

4.88 

0 

5 

1.76 

153H 

LT 

HOU/SW 

33 

0.53 

2.73 

0 

6 

7.47 

154H 

LT 

HOU/SW 

26 

0.75 

5.59 

0 

4 

2.08 

155H 

LT 

HOU/SW 

28 

0.77 

6.00 

0 

4 

2.31 

15f)H 

LT 

HOU/SW 

157H 

LT 

HOU/SW 

19 

0.46 

2.26 

0 

5 

2.51 

158H 

LT 

HOU/SW 

27 

0.82 

12.40 

1 

5 

6.72 

159H 

LT 

HOU/NW 

38 

0.70 

10.80 

4 

8 

10.74 

160H 

LT 

HOU/NW 

28 

0.61 

3.86 

0 

6 

3.44 

161H 

LT 

HOU/NW 

23 

0.40 

2.39 

0 

8 

5.97 

162H 

LT 

HOU/NW 

25 

0.55 

3.43 

0 

6 

8.78 

163H 

LT 

HOU/NW 

18 

0.61 

2.23 

0 

3 

3.84 

164H 

LT 

HOU/NW 

24 

0.46 

2.73 

0 

7 

4.76 

165H 

LT 

HOU/NW 

166H 

LT 

HOU/NW 

17 

0.48 

1.88 

0 

4 

7.16 

167H 

LT 

HOU/NW 

14 

0.85 

16.71 

0 

4 

1.33 

170H 

LT 

SKF 

28 

0.52 

2.53 

0 

5 

9.50 

171H 

LT 

SKF 

24 

0.60 

3.35 

0 

5 

1.73 
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Table  2.2  Results  of  Fitting  Shifted  Negative  Binomial  Distributions  to 
the  Number  of  Communications  Transactions  per  Aircraft  (ront.) 
(Houston  Individual  Sectors) 


Sector  ID 

No. 

Obs. 

A 

P 

k 

No. 

Adj 

d.f. 

2 

X 

172H 

LT 

SKF 

27 

0.48 

2.60 

0 

6 

5.81 

173H 

LT 

SAT 

17 

0.34 

1.23 

0 

5 

6.84 

174H 

LT 

SAT 

18 

0.75 

10.14 

0 

5 

4.50 

175H 

LT 

SAT 

19 

0.68 

3.97 

C 

4 

2.46 

176H 

LT 

SAT 

20 

0.67 

6.09 

0 

6 

2.94 

177H 

LT 

SAT 

30 

0.75 

10.75 

3 

7 

11.22 

179H 

LT 

8CS 

18 

0.79 

9.75 

0 

5 

4.25 

180H 

LT 

8CS 

12 

0.75 

8.64 

0 

4 

2.12 

181H 

LT 

8CS 

22 

0.35 

1.86 

0 

7 

5.62 

183H 

LT 

MTA 

21 

0.68 

3.84 

0 

4 

2.17 

184H 

LT 

MTA 

37 

0.43 

1.57 

0 

6 

4.46 

185H 

LT 

MTA 

18 

0.54 

2.52 

0 

4 

11.07 

187H 

LT 

MTA 

12 

0.37 

1.45 

0 

4 

5.59 

188H 

LT 

AUS 

11 

0.46 

2.50 

0 

4 

2.38 

189H 

LT 

AUS 

13 

0.55 

3.29 

0 

4 

4.48 

190H 

LT 

AUS 

23 

0.76 

4.49 

0 

3 

3.63 

191H 

LT 

AUS 

19 

0.58 

2.51 

0 

4 

4.02 

192H 

LT 

AUS 

39 

0.85 

16.30 

4 

6 

5.60 

194H 

HI 

MOB 

34 

0.70 

3.54 

0 

4 

16. 54’ 

195H 

HI 

MOB 

23 

0.48 

1.58 

0 

4 

4.95 

196H 

HI 

MOB 

27 

0.43 

1.02 

0 

4 

2.99 

197H 

HI 

MOB 

34 

0.45 

1.08 

0 

4 

7.20 

198H 

HI 

MSY/HI 

25 

0.41 

0.58 

0 

3 

3.60 

199H 

HI 

MSY/HI 

36 

0.57 

1.76 

0 

4 

2.30 

202H 

HI 

MSY/HI 

40 

0.53 

1.53 

0 

4 

3.71 

203H 

HI 

LCH/HI 

13 

0.22 

0.54 

0 

3 

6.94 

206H 

HI 

LCH/HI 

32 

0.84 

19.86 

2 

6 

4.32 

207H 

HI 

AEX/HI 

23 

0.24 

1.11 

0 

7 

8.30 

208H 

HI 

AEX/HI 

33 

0.47 

1.45 

0 

5 

2.25 

209H 

HI 

AEX/HI 

20 

0.67 

2.98 

0 

3 

1.47 

210H 

HI 

AEX/HI 

29 

0.46 

1.46 

0 

5 

2.91 

212H 

HI 

AEX/HI 

33 

0.48 

1.46 

0 

5 

10.43 
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Table  2.2 


Results  of  Fitting  Shifted  Negative  Binomial  Distributions  to 
the  Number  of  Communications  Transactions  per  Aircraft  (cont.) 
(Houston  Individual  Sectors) 


Sector  ID 

No. 

Obs. 

A 

P 

A 

k 

No. 

Adj. 

d.f. 

2 

X 

213H 

HI 

HOU/HI 

15 

0.61 

1.38 

0 

2 

4.45 

215H 

HI 

HOU/HI 

20 

0.61 

2.31 

0 

4 

9.72 

216H 

HI 

HOU/HI 

31 

0.80 

4.73 

0 

3 

5.58 

217H 

HI 

HOU/HI 

27 

0.43 

1.31 

0 

5 

9.19 

218H 

HI 

AUS/HI 

54 

0.93 

29.25 

4 

5 

2.00 

219H 

HI 

AUS/HI 

18 

0.25 

1.21 

0 

7 

6.30 

220H 

HI 

AUS/HI 

29 

0.66 

4.22 

0 

5 

12.76 

221H 

HI 

AUS/HI 

28 

0.57 

3.07 

0 

5 

6.39 

222H 

HI 

AUS/HI 

21 

0.18 

0.55 

0 

6 

24. 9l' 

223H 

HI 

SAT/HI 

224H 

HI 

SAT/HI 

17 

0.37 

3.51 

0 

9 

14.90 

225H 

HI 

SAT/HI 

19 

0.29 

1.59 

0 

7 

8.06 

226H 

HI 

SAT/HI 

33 

0.31 

1.95 

0 

10 

16.71 

227H 

HI 

SAT/HI 

20 

0.69 

6.11 

0 

5 

5.03 

Table  2.3  Results  of  Fitting  Shifted  Negative  Binomial  Distributions 
to  the  Number  of  CT/AC  (Houston  Enroute  Functions) 


Sector  Function 

No. 

Obs. 

A 

P 

k 

No. 

Adj. 

d.f. 

2 

X 

LE 

206 

0.42 

2.01 

0 

11 

10.42 

LI 

1412 

0.56 

2.98 

0 

12 

9.30 

HI 

704 

0,36 

1.14 

0 

13 

33.46’ 

HI  Group  1 
(194H-217H) 

482 

0.51 

1.59 

0 

9 

8.63 

HI  Group  2 

222 

0.31 

1.41 

0 

13 

26.37 

(218H-227H) 


* Significant  at  the  0.01  level 


2-17 


Table  2.4  Comparison  of  the  Estimates  of  p and  k 

for  Houston  and  New  York  Enroute  FunctlopR 


Houston  New  York 


A 

A 

A 

A 

A 

P 

k 

E(X) 

P 

k 

E(X) 

LE 

0,42 

2.01 

4.7757 

0.46 

2.30 

4.7000 

LT 

0.56 

2.98 

4.3414 

0.41 

2.36 

5.3961 

HI 

(group 

1) 

0.51 

1.59 

3.5276 

0.58 

4.03 

4.9183 

HI 

(group 

2) 

0.31 

1.41 

5.1384 

*E(X) 

= k(  f -1) 

+ m 
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I 
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Table  2.5  Comoarison  of  Historical  versus  Simulated  Aircraft 

Loading  Series  Generated  from  Sector  Function  Simulations 
(New  York  Metroplex) 


Sector 

Function 

Series 

A 

02 

^2 

a 

t^ 

LT 

Historical 

(453LT) 

4.4988 

1.1131 

-0.2015 

0.5884 

-0.7027 

1.9051 

Simulated 

3.7909 

1.1379 

-0.2363 

0.4313 

LE 

Historical 

(464LE) 

5.4113 

1.1533 

-0.2511 

0.6567 

0.3245 

0.6908 

Simulated 

5.7026 

1.1416 

-0.3038 

0.7000 

HI 

Historical 

(475HI) 

4.4104 

1.2874 

-0.4048 

0.5655 

0.7049 

1.0887 

Simulated 

4.9354 

1.1647 

-0.3319 

0.6849 

GN 

Historical 

(504GN) 

2.4156 

0.9155 

-0.0698 

0.5360 

0.3682 

1.6614 

Simulated 

2.6585 

1.1470 

-0.2833 

0.5362 

LC 

Historical 

(510LC) 

6.3220 

0.9737 

-0.1314 

0.5271 

-0.2922 

4.3406 

Simulated 

6.0167 

1.1618 

-0.2510 

0.8592 

LG 

Historical 

(524LG) 

2.9971 

1.2741 

-0.3681 

0.4132 

-0.0460 

0.2286 

Simulated 

2.9472 

1.2740 

-0.3446 

0.4607 

DP 

Historical 

(534DP) 

1.8776 

0.9404 

-0.2310 

0.3680 

1.8743 

5.3551 

Simulated 

2.8147 

1.1990 

-0.3514 

0.5868 

AD 

Historical 

(543AD) 

1.3789 

1.2273 

-0.3855 

0.2247 

1.1687 

3.4071 

Simul.ited 

2.0835 

1.1577 

-0.2600 

0.3558 

AR 

Historical 

(553AR) 

3.3706 

1.0482 

-0.0891 

0.4028 

-0.1142 

1.8212 

Simulated 

3.1905 

1.2261 

-0.3093 

0.3791 

a 

t 

.005<200)  = 

2.6 

*’l  2 

2 ~ 
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Table  2.6  Comparison  of  Historical  versus  Simulated  Channel  Utilization 


Series 

Generated  from  Sector 
(New  York  Metroplex) 

Function 

S imulations 

Sector 

function 

Series 

P 

<(>1 

4)2 

-^2 

a 

t 

a 

G(0,1)‘’ 

LT 

Historical  (453LT) 

0.4939 

0.1303 

0.2064 

0.0510 

-0.7299 

2.4383 

Simulated 

0.4534 

0.3496 

0.1231 

0.0689 

LE 

Historical  (464LE) 

0.4806 

0.3763 

0.1823 

0.0542 

1.3919 

0.4348 

Simulated 

0.5837 

0.4587 

0.1259 

0.0638 

HI 

Historical  (47 SHI) 

0.5740 

0.3502 

0.1485 

0.0491 

-0.8195 

0.9569 

Simulated 

0.5177 

0.4932 

0.1022 

0.0593 

GN 

Historical  (504GN) 

0.4964 

0.5094 

0.0667 

0.0742 

-1.6846 

1.4124 

Simulated 

0.3691 

0.2920 

0.1570 

0.0795 

LC 

Historical  (510LC) 

0.5595 

0.2194 

0.0407 

0.0470 

-2.2771 

2.0329 

Simulated 

0.4463 

0.4568 

0.0586 

0.0484 

LG 

Historical  (524LG) 

0.3096 

0.2989 

0.2286 

0.0520 

-0.1488 

0.6368 

Simulated 

0.2992 

0.4112 

0.1625 

0.0616 

DP 

Historical  (534DP) 

0.4379 

O.25I3 

0.2124 

0.0441 

0.6749 

5.6940 

Simulated 

0.4846 

0.5330 

0.0438 

0.0741 

AD 

Historical  (S43AD) 

0.5760 

0.3171 

0.1153 

0.0840 

-2.9322* 

1 7839 

Simulated 

0.3366 

0.4808 

0.1480 

0.0726 

AR 

Historical  (553AR) 

0.5521 

0.3268 

0.1879 

0.0574 

-0.1614 

0.9875 

Simulated 

0.5396 

0.5003 

0.1434 

0.0588 

a 

b 1 2 
2 ^.01 

(3)  = 5.6 
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Table  2.7 


Comparison  of  Historical  versus  Simulated  Aircraft 
Loading  Series  Generated  from  Sector  Function  Simulations 
(Houston  ARTCC) 


Sector 

Function 

Series 

u 

92 

-2 

a 

a 

t 

0(^,1)'’ 

LE 

Historical 

(LE107H) 

4.3157 

1.1472 

-0.1994 

0.2471 

-1.3858 

4.7444 

Simulated 

3.0292 

1.2237 

-0.3926 

0.3359 

LT 

Historical 

(LT145H) 

3.2067 

1.2395 

-0.3017 

0.2959 

-1.3546 

2.6192 

Simulated 

2.0489 

1.1751 

-0.3557 

0.3374 

HI 

Historical 

(HI226H) 

3.0424 

1.2056 

-0.2268 

0.2782 

-0.1192 

4.2715 

Simulated 

2.7642 

1.2418 

-0.3599 

0.3723 

"^005(200)  = 

2.6 

1 

2 ^ 

5.6 
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Table  2.8 


Comparison  of  Historical  versus  Simulated  Channel  Utilization 
Series  Generated  from  Sector  Function  Simulations 


(Houston  ARTCC) 


Sector 

Function 

Series 

M 

^2 

a 

t^ 

G()2,l)' 

LE 

Historical  (LE107H) 

.2642 

.1861 

,0901 

.0529 

.5767 

2.0493 

Simulated 

.2934 

.4425 

.0102 

.0601 

LT 

Historical  (LT145H) 

.3066 

.2926 

.1344 

.0543 

-.4998 

0.8263 

Simulated 

.2799 

.4150 

-.0080 

.0603 

HI 

Historical  (HI226H) 

.3046 

.3924 

.2049 

.504 

.1511 

0.3541 

Simulated 

.3151 

.3176 

.2039 

.0578 

a 


*^.005 


(200) 


2.6 


b 


1 

2 


2 

01 


fr^=s 


6 
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CHAPTER  3 


r 


[ 

I 

I FURTHER  APPLICATIONS  OF  THE  SIMULATION  MODEL 


3.1  Introduction 


In  the  evolution  of  the  ATC  simulation  model,  several  consider- 
ations have  guided  the  technical  developments: 

(i)  the  structure  of  the  simulation  procedure  should  be  general, 
but  adaptive  to  specific  applications  as  well, 

(li)  the  formulation  of  the  computer  program  should  be  simple, 
but  allow  additional  details  to  be  added. 

Experience  gained  from  many  exercises,  and  the  feedback  from  comparisons 
against  historical  data  have  resulted  in  many  modifications  of  the  simu- 
lation model.  The  simulation  has  been  extended  from  the  one  appropriate 
to  individual  sectors  to  a general  simulation  for  sector  functions  without 
requiring  extra  computational  effort  Applications,  either  for  hypo- 
thesized changes  of  the  system,  or  for  the  study  of  events  with  practical 
significance  should  continue  to  contribute  to  additional  generalizations 
and  understanding  of  the  simulation  model. 

In  this  chapter  three  different  applications  of  the  simulation 
model  are  presented  to  Illustrate  the  usefulness  of  the  model.  In  section 

3.2  a continuing  study  of  estimating  communications  capacities  for  st?ctor 
functions  is  described  and  results  analyzed  in  conjunction  with  some  con- 
ventional queuing  models.  Section  3.3  is  devoted  to  the  exploration  of 
the  effects  of  a proposed  FCC  regulation  concerning  tone  burst  identifi- 
cation of  privately  owned  transmitters,  for  ATC  alr/'ground  communications. 
Various  response  surfaces  are  constructed  on  the  two  causal  factors:  the 
tone  duration  and  traffic  density,  from  the  results  of  simulation  experi- 
mentation and  spline  extropolatlons.  In  section  3.4  a hypothesized  change 
in  the  distribution  of  //TR/CT  using  a smooth  weighting  function  is 
Investigated  and  results  analyzed. 

3.2  imat  1 on  of  Communlcati ons  Capacities  for  New  York  Sector 

Funct  ions 

As  earlier  demonstrated  in  the  report  volume  III,  chapter  4, 
techniques  useful  for  detecting  nonstatlonarity  of  the  response  time 
series  can  be  used  to  measure  capacities  of  communication  channels.  In 
order  to  gain  more  insights  about  the  mechanism  which  determines  the 
communication  capacity  and  to  explore  the  effects  of  regulating  Input 
variables  (e.g.  the  transaction  lengths,  communication  arrival  streams, 
etc.),  some  well-known  results  from  the  conventional  queuing  theory  are 
discussed  in  this  section  and  Information  particularly  important  for  the 
design  of  simulation  experiments  gathered. 
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In  the  following  discussions  we  first  compare  the  ATC  communi- 
cation queuing  with  a theoretical  model  of  Poisson  arrival  streams  and 
general  service  time  distributions.  Effects  of  the  nonhomogeneity  of 
communication  arrivals  upon  queuing  are  then  discussed.  Influences  of  the 
variance  reduction  of  communications  transaction  lengths  are  studied  and 
results  presented.  Finally,  a simulation  experiment  for  estimating  sector 
function  communications  capacities  is  described  and  results  tabulated  for 
practical  use. 


» 

3.2.1  A Theoretical  Model  for  ATC  Communications  Queuing 

In  chapter  2 of  report  \A)lume  III,  an  exponential  model  was 
analyzed  and  verified  for  intercommunication  gap  lengttis.  It  was  demon- 
strated that  the  parameter  values  descriptive  of  the  gap  length  associated 
with  a specific  aircraft  were  a function  of  the  nuiuber  of  transactions  made 
by  tiie  aircraft.  By  further  assuming  independence  between  Individual  gap 
lengths,  a theoretical  model  for  the  communication  arrivals  to  an  ATC 
communication  channel  can  be  formulated.  To  explain: 

(1)  For  a single  aircraft,  arrivals  to  the  communications  channel 
follow  a homogeneous  Poisson  process.  The  value  of  the  arrival  intensity 
(i.e.  the  value  of  the  Poisson  parameter  A.  for  the  i*-^  aircraft)  is 
determined  by  the  number  of  communication  transactions  associated  with  this 
alrcraf  t . 


(ii)  When  there  are  many  aircraft  handled  by  a controller  at  one 
time,  and  within  a moderate  short  time  length  where  no  aircraft  enters  or 
leaves  the  sector,  the  combined  communication  arrival  stream  is  still 
homogeneous  Poisson.  The  intensity  of  this  combined  process  is  equal  to 
the  sum  of  the  intensities  of  all  aircraft  presently  in  the  sector  (i.e. 

x*=?;x.). 

(ill)  Wlienever  an  aircraft  enters  or  leaves  the  sector,  tlie 
Intensity  of  the  combined  communication  arrival  stream  X changes  depending 
on  the  number  of  aircraft  staying  in  the  sector  and  their  associated 
communication  Intensities.  From  a simple  simulation  experiment  it  was 
found  that  X follows  a second  order  autoregressive  model,  when  jjata  based 
on  per-minute  averages  were  examined.  In  fact,  the  values  of  X follow 
the  so-called  immigration-death  process.  (For  reference,  see  Cox  and 
Miller:  The  Theory  of  Stochastic  Processes,  1966,  John  Wiley  Co.) 

A schematic  presentation  of  the  process  described  above  is  given 
in  figure  3.1.  From  the  time  plot  for  X , the  communication  arrival 
intensity  for  the  combined  process,  one  can  see  differences  between  tlie 
arrival  stream  to  an  ATC  communication  channel  and  those  provided  by  the 
conventional  assumption  of  a homogeneous  Poisson  process  which  has  a 
constant  X throughout  the  entire  time  period. 

As  has  been  discussed  in  earlier  cliapters,  the  communication  trans- 
action lengths  are  gamma  variables  (for  the  transmission  lengths)  com- 
pounded with  a distribution  for  the  number  of  transmissions  per  communl cat  ion 
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transact  Ion.  In  tlie  literature  of -queuing  theory,  the  compounded 
distribution  is  called  the  service  time  distribution.  Combining  a non- 
homogeneous  Poisson  arrival  stream  with  the  distribution  of  transaction 
lengths,  a complete  queuing  model  for  the  ATC  communications  is  specified. 

Although  analytic  solutions  to  such  models  are  difficult  to 
attain,  this  formulation  is  helpful  in  providing  insights  about  the 
queuing  properties  of  communication  transactions.  Particularly,  solutions 
to  a model  consisting  of  a homogeneous  Poisson  arrival  and  a general 
service  time  distribution  can  be  considered  as  providing  lower  bounds  to 
the  expected  queuing  times  and  queue  sizes.  These  lower  bounds  are  use- 
ful for  selecting  ranges  for  simulation  experiments. 


3.2.2  Communications  Capacity  and  Queuing 

It  can  be  shown  that  the  expected  channel  utilization,  expressed 
by  number  of  seconds  per  hour,  is 

E(CU)  = 0[E(TR)  ] fE(//TR/CT)  1 [E(//CT/AC)  1 

= 0[  b -b ,0] [t] [k(-  -l)+ml  (3.1) 

o 1 p 

where 

E = Expectation  Operator 

0 = the  number  of  aircraft  arrivals  per  hour; 

b ,b,  : the  parameters  used  to  relate  E(TR)  with  O; 
o 1 

T : the  expected  number  of  transmissions  per  CT; 
k,p,m  : parameters  in  the  distribution  of  #CT/AC; 


A straightforward  determination  of  the  sector  communication  capacity  is 
to  set  E(CU)  equal  to  3600  seconds  x 0.85,  the  physical  limitation  for  the 
period  of  1 hour,  and  then  to  solve  the  Equation  (3.1)  for  0.  In  other 
words,  the  value  of  0 which  satisfies  the  equation: 

0(b  -b,0)  = 3600  (0.85)*t~^-  fk(--l)+  m]"^  (3.2) 

o 1 P 

is  the  maximum  number  of  aircraft  arrivals  tlieoretlca lly  allowed. 

We  now  turn  to  a discussion  of  queuing  associated  with  Equation 
(3.2).  Let  us  hypothesize  that  one  is  able  to  stabllze  tlie  communicat ion 
arrivals  to  an  ATC  communication  channel  such  that  the  intensity  X stays 
constant  over  a fairly  long  period  of  time  The  results  of  conventional 
queuing  theory  can  now  be  used  to  explore  the  nature  of  the  communication 
capacity. 


It  is  \;e' ’-known  In  the  literature  of  queuing 
expected  queuing  time  for  an  M/C,/]  model  (the  system 


tlieory,  that  the 
consists  of  a 
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Poisson  arrival  stream,  a general  service  time  distribution,  and  a single 
server) , in  a long  run  equilibrium,  is 

E(Q)  = expected  queuing  time 


2 2 

1 ■ 

- X*  ^ 2(l-p) 


expected  service  time  (1.3) 


where 


X ; parameter  of  the  homogeneous  (after  stabilization, 
in  our  case)  Poisson  arrival  process; 

p : a traffic  intensity  (=  mean  service  tirae/mean  inter- 
arrival time); 

C,  : the  coefficient  of  variation  of  the  service  time 
distribution  (=  standard  deviation/mean). 


One  interesting  fact  is  that  p is  equal  to  E(CU)/3600.  This  can  be  seen 
from  the  following  algebraic  manipulations.  Let  us  write 

- . . _ mean  service  time 

' ^ mean  communication  Interarrival  time 


(li)  mean  service  time  = mean  communication  transaction 

length  = E(T1^E(//TR/CT)  = (b  -b,0)*T 

o i 

(ill)  mean  communication  interarrjval  time  = “ 

36OO/[0*E((<CT/AC)]  = 3bOO.0~^  [k(i  - l)+m]"^ 

r 

Substituting  (ii)  and  (iii)  into  (1),  one  finds 


p = 0. (b  -b,0>T*[k(-  -l)+ml/3600 

^ o 1 p 

= E(CU)/3600 

Q. E.D. 


Furthermore,  since 

X*  = 0[k(-  -l)-Hn]/3600 
P 

the  p can  be  expressed  as 


= X 


(3.4) 


(3.5) 


>-5 


P 


(b  -b  0)*T 

o 1 


(3.6) 


From  Equation  (3.4)  one  may  see  that,  to  find  the  value  of  0 which  f>,ives 
p = 0.85,  one  has  to  refer,  once  again  to  Eqj^ation(3. 2) • Now  let  us  denote 
the  solution  value  of  0 in  Equation^!, 2)by  0,  The  complete  formula  for 
E(Q),  when  p = 0.85,  is  thus 


E(Q) 


'^2  ^1.^2, 

^ P 

P + 

2(1-P) 


- (b^-bj^^)T 


where 

II 

S[k(-  -1)4-  Dl 

P 

1/3600 

(=0.85 

p = 

^ (b 

0 1 

(3.7) 


An  illustration  of  the  calculation  of  E(Q)  for  HI  sectors  is  given  below. 

From  previous  reports  volume  IT  and  III,  the  parameter  estimates 
from  an  HI  sector  function  are  as  follows: 


b = 3.70 
o 

b^  = 0.0158 
T = 3.56 
k = 4.03 


p = 0.58 
m = 2 

a = 0.8054 

b 


Substituting  these  values  into  Equation  (3.2),  one  obtains 


^•'3.70  - 0.0158b)  = 3600  (0.85)  (3.56)"^  [4.03(0. 724)+2] 
= 3600  (0.85)  (0.281)  (0.203) 

= 174.85 

The  above  equation  can  be  rewritten  as 


Thus , 


- 234.180  + 11066.46  = 0 
0 = 65o675  ^66 


The  second  root  of  0 is  Ignored  because  the  value  of  E(CU)  is  higher  than 
3600x0.85.  The  figure  of  66  is  thus  the  theoretical  maximum  number  of 
aircraft  allowed  to  enter  an  average  HI  sector,  per  hour,  without  exploding 
the  system.  This  does  not  mean  that  one  may  actually  load  this  many  air- 
craft into  the  channel  at  the  specified  rate.  Another  important  response, 
the  communications  queuing  times,  which  is  also  closely  related  to  the 
aircraft  arrival  rate,  should  be  simultaneously  considered. 


It  must  be  emphasized  that  the  expected  queuing  time  for 
communication  channel,  which  has  a nonhomogeneous  Poisson  arrival 


an  ATE 
St  ream. 
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is  still  unknown  to  us.  A lower  bound  of  this  expected  queuing  time  can 
be  calculated  from  Equation  (3,7).  Again,  by  substituting  all  estimates 
of  the  parameters  of  concern  into  Equation  (3.7),  we  obtain 


0 [k(-  -l)+m]/3600 
P 


= 66  [4.03(0.724)+2]/3600 


Thus , 


= (0.09)  [0.85  + -- 


(0.85)~(l+0.65). 


2(1-0.85) 


[3. 70-(0. 0158)66] (3.56) 

= 44.10  (seconds) 

This  lower  bound  for  E(Q)  indicates  that  when  an  HI  sector  channel  is 
loaded  at  66  aircraft  per  hour,  the  average  queuing  time  per  communi- 
cat ion  is  at  least  44.10  seconds.  Comparing  this  result  with  a simu- 
lation result  previously  reported  in  volume  III,  chapter  4,  we  found 
that  the  actual  average  queuing  time  was  equal  to  56.82  seconds,  where 
50  aircraft  were  loaded  per  hour.  Clearly  the  expected  queuing  time 
derived  from  a hypothesized  homogeneous  Poisson  process,  in  particular, 
the  M/G/1  model,  is  Indeed  too  small  but  it  is  still  valid  if  viewed  as 
a lower  bound. 

3.2.3  Effects  of  Regulating  (Communications  Transaction  Lengths 

One  interesting  problem  related  to  the  study  of  communication 
capacity  concerns  how  much  Improvement  in  queuing  can  be  anticipated 
if  one  is  able  to  reduce  the  variability  of  the  communication  trans- 
action lengths.  Although  a precise  solution  to  this  problem  is  still 
unavailable,  the  effects  of  regulating  transaction  lengths  can  be 
elucidated  by  studying  the  shifts  of  the  lower  bound  of  the  expected 
queuing  time,  E(Q),  at  different  values  of  the  coefficient  of 
variation  of  the  transaction  lengths.  Using  the  same  example  as  in  the 
preceding  section,  with  p = 0.85,  we  compute  the  E(Q)  at  various  values 
of  ranging  from  0 to  1 . The  results  are  displayed  in  figure  3.2. 

From  this  plot  one  may  see  the  lower  bound  of  44.10  seconds  has  dropped 
to  the  level  of  26.74  seconds  when  is  reduced  from  its  present  value 
of  0.8  to  0.0.  This  is  a reduction  of  3d  per  cent.  In  general,  regulating 
the  transaction  lengths  will  help  reduce  queuing  times,  particularly  under 
heavy  traffic  situations. 


3.2.4  Results  of  a Simulation  Experiment 


ATC  communications  simulations  for  New  York  and  Houston  ARTEC 
sector  functions  were  validated  and  results  reported  in  chapter  2. 


ing  time  (seconds) 


0 .25  .50  .75  1.0  Cb 

Coefficient  of  variation  of  transaction  lengths 


Figure  3.2  Lower  Bound  of  the  Expected  Queuing 
Time  vs  Cb  (HI  Sectors) 


Making  use  of  this  verified  model,  we  proceeded  to  estimate  the  communi- 
cations capacities  by  varying  the  expected  number  of  aircraft  arrivals 
per  hour  in  a sequence  of  simulation  exerclseso  The  capacity  measure- 
ments were  made  by  simultaneously  considering  the  following  five  variables: 

(1)  Average  aircraft  loading,  n^; 

(li)  Average  channel  utilization,  c^; 

(ill)  Average  queuing  time  per  communication  transaction,  Q^; 

(iv)  The  degree  of  nonstationarlty  measured  from  the  time  series  of 
n^,  denoted  by  Pr(f;>l); 

(v)  The  probability  of  n^  being  larger  than  a critical  value  k, 
denoted  by  Pr(n^>k). 

Capacity  estimates  through  experimentation  with  the  simulation  model  were 
obtained  for  eight  different  sector  functions,  namely,  LT,  LE,  HI,  LC,  LG, 
DP,  AD,  and  AR.  Details  of  these  experiments  are  given  below.  A special 
study  on  the  sector  function  GN,  using  an  analytical  formula,  is  presented 
later  in  this  subsection. 

For  all  simulation  runs  actually  performed,  the  three  rules  listed 
below  were  followed: 

(i)  The  simulation  is  run  for  a period  of  1-hour  (simulation  time) 
before  any  record  is  taken;  this  is  to  eliminate  the  effects  of 
initial  conditions. 

(ii)  The  pseudo  random  numbers  used  are  kept  identical  whenever 
possible  for  the  simulations  of  the  same  sector  function. 

(iii)  For  each  of  the  computer  runs  of  2-hour  simulation  time  the 

three  major  responses:  the  aircraft  loading  (n  ),  the  channel 
utilization  (c^),  and  the  queue  lengths  (q  ) are  recorded  as 
time  series,  on  the  basis  of  60-second  averages. 

Our  simulation  experiments  for  exploring  the  communications 
capacities  were  performed  in  an  Iterative  fashion  as  described  below. 

An  initial  traffic  density  was  selected,  usually  being  equal  to  or 
slightly  higher  than  the  value  observed  from  the  historical  data.  The 
experimentation  with  the  simulation  model  is  then  carried  out  bv 
increasing  the  expected  number  of  aircraft  arrivals  per  hour  until  the 
average  channel  utilization  generated  reached  0.85. 

Summary  statistics  of  the  responses  of  concern  are  furnished 
in  table  3.1  (where  the  "bar"  on  the  top  of  the  variable  means  "averaged 
over  the  2-hmir  data").  The  nonstatlonarltv  detection  measure  Pr(f,>l), 
where  E = + $2*  which  has  been  described  In  the  report  volume  III, 

(Chapter  4,  was  applied  to  the  simulated  aircraft  loading  time  series, 
n^.  The  numerical  results  of  this  measure  under  various  experimental 
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Table  3.1  Summary  Results  of  the  Simulation  Experiments 
for  Capacity  Study  (New  York  Sector  Functions) 


• Function 

Traffic  Density 
()5fac/hr) 

Expected  Actual 

n 

t 

"t 

Qt 

(sec.) 

Pr(C>l) 

LT 

40 

40 

6.628 

.689 

17.247 

.0031 

45 

44 

7.401 

.775 

23.273 

.0038 

55 

49 

10.204 

.844 

46.112 

.0262* 

LE 

35 

39 

5.852 

.591 

13.357 

.0001 

45 

44 

7.669 

.636 

24.660 

.0057 

55 

58 

9.362 

.748 

30.638 

.0092 

70 

68 

16.458 

.857 

81.866 

.0696* 

HI 

45 

52 

8.216 

.728 

21.017 

.0000 

50 

57 

9.555 

.800 

28.629 

.0032 

60 

68 

16.640 

.869 

89.104 

.3749* 

LC 

85 

74 

11.652 

.760 

23.953 

.0018 

100 

93 

26.157 

.877 

120.406 

.0298* 

110 

101 

45.195 

.877 

265.149 

.3001* 

LG 

80 

77 

10.995 

.730 

24.531 

.0095 

90 

93 

15.079 

.786 

46.835 

.0357* 

100 

108 

20.097 

.822 

72.632 

.4555* 

DP 

50 

44 

6.137 

.686 

19.205 

.0034 

55 

50 

6.895 

.745 

21.293 

.0032 

60 

56 

8.059 

.807 

28.792 

.0015 

70 

66 

12.573 

.875 

62.698 

.0590* 
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Table  3.1  Summary  Results  of  the  Simulation  Experiments  for 
Capacity  Study  (New  York  Sector  Functions)  (cont.) 


Sector  Function 

Traffic 

(//ac/hr) 

Expected 

Density 

Actual 

n 

t 

c 

t 

^t 

(sec. ) 

Pr(C>r 

AD 

40 

34 

4.089 

.603 

13.532 

.0017 

45 

42 

4.694 

.665 

15.514 

.0000 

50 

43 

5.865 

.714 

25.307 

.0040 

60 

56 

7.095 

.773 

33.039 

.0007 

70 

67 

11.781 

.873 

76.798 

.0796 

AR 

35 

34 

6.443 

.780 

27.855 

.0037 

40 

42 

9.079 

.873 

49.853 

.0153' 

Indicating  critical  conditions 
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conditions  are  displayed  in  the  last  column  of  table  3.1.  In  this  study 
we  use  Pr(C>l)  = 0.01  as  a critical  point  which  seems  sensible,  although 
somewhat  arbitrary,  for  the  detection  of  the  onset  of  an  unstable 
operation.  From  the  results  displayed  in  table  3,1,  capacities  of  each 
control  function  in  N.Y.  metroplex  were  decided  and  listed  in  table  3.2. 
Additional  measurements  of  the  operational  instability,  expressed  by  the 
probability  Pr(n^>k),  projecting  from  the  fitted  AR(2)  process,  where  k 
is  a predetermined  value,  are  listed  in  tables  3.3-3.10  for  each  of  the 
eight  functions  studied.  Tables  containing  detailed  statistics  computed 
from  the  simulated  data,  as  well  as  the  time  series  plots  of  the  major 
responses  are  compiled  and  available  as  a separate  Interim  report  (PU-37). 


3.2.5  A special  Note  oti  the  Capacity  of  the  ON  Function 

Because  of  the  low  average  value  on  the  number  of  CT's  per  air- 
craft (=2,27)  and  the  relatively  high  speed  of  decrease  on  the  expected 
transmission  length  (expected  TR  length  ^ 3.70  - 0.0158* (expected  number 

of  aircraft  arrivals  per  hour)),  the  GN  function  has  an  extraordinarily 
large  communications  capacity.  Due  to  the  limitation  of  the  storage  space 
in  a computer,  a simulation  experiment  is  difficult  to  carry  out  for  this 
control  function.  Nevertheless,  an  analytic  approximation  to  the  ATC 
channel  model  is  available  for  the  estimation  of  the  capacity.  Formulas 
useful  for  this  purpose  are  sketched  below.  Related  numerical  results  are 
displayed  in  figures  3.3  and  3.4.  From  these  graphs  0 = 180  aircraft 
arrivals/hour  seem  to  be  a moderate  estimate  of  the  capacity  of  the  (!N 
function. 


Analysis  of  the  GN  function  capacity: 

E(#CT/AC)  = k(-  -1)  + m = 0.88(2.439  -1)  + 1 = 2.2663 
P 

E((tTR/CT)=  3.4695 

E(TR)  = 3.70  - 0.01580,  0^  = 0.99  - (0-  50)  0.005 

Additional  Conditions: 

(i)  TR  length  on  average  > 2.0  j 

(li)  0.7  < C,^  < 0.99 

— b - 


p = 0 E(TR)E(/iTR/CT)E(i'/CT/AC)/3600 
X = 2.2663  0/3600 


E(Q) 


X 2(l-p)  J 


- E(TR 


) 3,4695 


Both  p and  E(Q)  are  plotted  on  graphs  for  various  traffic  densities 
(#ac/hr) . 
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Table  3.2 


Estimates  of  the  Communications  Capacities 


for  New  York  Sector  Functions 


Sector  Function 

Expected  Capacity 
(//  aircraft/hour) 

+ 

Conservative  Capacity 
(//  aircraft/hour) 

LT 

50 

41 

LE 

55 

46 

HI 

50 

41 

GN 

180* 

163 

LC 

90 

78 

LG 

85 

73 

DP 

60 

50 

AD 

60 

50 

AR 

37 

29 

Estimated  from  an  analytic  solution  , 

Conservative  Capacity  = Expected  Capacity  - 1 . 282 (Expected 

Capacity)^' 

= 90%  confidence  range 


L. 
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0 50  100  150  200 

Number  of  aircraft  arrivals  per  hour  (ac/hr) 

Figure  3.4  Expected  Queuing  Time  as  a Function  of 
ac/hr  (GN  Sectors) 
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F,£fects  of  Tone  Disturbance  upon  ATC  Communications 

This  section  describes  the  results  of  simulation  experiments 
designed  to  show  the  effect  of  the  proposed  FCC  regulation  concerning  tone 
burst  identification  of  privately  owned  transmitters  on  air  traffic  control 
(ATC)  air/ground  communications.  The  nine  functions  subjected  to  investi- 
gation are:  HI,  LE,  LT,  GN,  LC,  LG,  DP,  AD  and  AR  selected  from  the  N.Y. 

area.  The  results  for  a particular  function  demonstrate  the  effect  that 
the  proposed  regulation,  if  implemented,  will  have  upon  communications  in 
a typical  sector  within  that  function.  The  complete  results  for  each 
function  are  contained  in  20  pages  of  tables  and  graphs,  thereby  giving  a 
total  of  180  pages  of  output.  In  this  report  only  a complete  set  of 
material  for  HI  function  is  included  for  illustiation.  Instead,  the 
summary  statistics  and  spline  extrapolated  response  surfaces  for  each  of 
the  nine  sector  functions  are  presented.  A complete  collection  of  the 
material  related  to  this  study  is  available  as  a NAPEC-Pr inceton  Univer- 
sity Interim  Report  (PU-39). 


3.3.1  Simulation  Input  Variables  and  the  Average  Responses 

In  this  study,  a simulation  experiment  is  characterized  by  two 
basic  input  parameters,  namely,  tone  duration  and  expected  aircraft 
arrival  rate  per  hour.  In  the  present  study,  simulation  experiments  were 
conducted  at  tone  durations  of  0,  .5,  1,  and  1.5  seconds.  A tone  burst 
appeared  at  the  beginning  and  end  of  each  pilot-initiated  transmission. 
Since  only  a negligible  number  of  transmissions  last  for  more  than  30 
seconds,  no  attempt  was  made  to  incorporate  tone  bursts  at  intervals  of 
30  seconds  within  the  main  body  of  a transmission  as  proposed  by  FCC. 

Hence,  the  effect  of  the  tone  burst  was  to  increase  the  length  of  pilot- 
initiated  transmissions  by  0,  1,  2,  or  3 seconds.  Average  aircraft  arrival 
rates  per  hour  for  each  of  the  typical  sectors  involved  in  table  2.5, 
estimated  from  data  collected  in  1969,  were  used  in  each  simulation  con- 
ducted for  this  study.  In  addition,  all  the  remaining  parameters  were 
obtained  from  master  equations  and  tables.  Thus  a total  of  four  2-hour 
experiments  were  performed  for  each  of  the  nine  sector  functions.  The 
results  of  each  experiment  were  originally  displayed  on  four  pages  of  out- 
put consisting  of  one  summary  page  followed  by  three  graphs  exhibiting  the 
dynamic  behavior  of  communications  during  the  2-hour  observation  period. 
However,  only  a grand  summary  of  all  simulations  is  presented  for  each  of 
the  sector  functions  as  shown  in  tables  3.11  - 3.19,  (A  sample  of  the  raw 
tables  and  graphs  is  given  at  the  end  of  this  chapter.) 

From  a careful  inspection  of  the  tables,  one  is  not  intrigued  by 
the  fact  that  the  average  number  of  aircraft  in  sector  n increases  with 
increasing  tone  duration  at  constant  aircraft  arrival  rafe.  In  practice, 
it  might  be  expected  that  pilots  and  controllers  will  alter  the  flow  of 
information  to  compensate  for  the  additional  burden  imposed  upon  the 
channel  by  long  tone  hurst  lengths  rather  than  resort  to  such  measures  as 
speed  control,  stacking,  etc.  On  the  other  hand,  the  simulation  program 
was  developed  from  data  collected  in  a normal  busy  working  environment 
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rather  than  a crisis  situation.  As  a result,  it  does  not  presently 
provide  for  compensation  via  alteration  of  information  flow,  and  so  the 
time  interval  of  an  aircraft's  being  in  the  sector  Increases  with 
increasing  tone  burst  duration.  Thus,  the  simulation  results  must  be  inter- 
preted under  the  assumption  that  th2  amount  ot  aTC  information  flow  between 
an  aircraft  and  ground  is  Independent  of  tone  burst  length.  Succinctly 
stated,  system  degradation  due  to  Increasine  tone  durations  must  be 
reflected  by  the  three  response  variables,  number  of  aircraft  in  sector, 
channel  state,  and  aircraft  queuing  state,  or  a reduction  in  the  quality 
and/or  quantity  of  ATC  air/ground  information  flow.  In  this  study,  the 
degradation  is  displayed  by  changes  in  the  three  response  variables. 

However,  this  is  not  to  say  that  degradation  in  terms  of  information  flow 
cannot  be  incorporated  into  the  simulation  software.  Obviously,  once  a 
mechanism  of  compensation  by  alteration  of  information  flow  is  specified, 
the  possibility  of  incorporation  arises.  In  fact,  the  development  of  the 
present  ATCC  simulation  program  was  undertaken  for  the  purpose  of  studying 
the  effect  that  luch  variations  from  normal  operating  conditions  have  on 
communications. 


3.3.2  Estimation  of  Average  Response  at  Intermediate  Tone 
Durations 

Since  the  GPSS  program  can  execute  simulations  for  oniy  Integer 
number  of  seconds  of  transmission  length,  the  four  trials  for  each  of  the 
sector  function  simulation  experiments,  at  the  tone  durations  of  0.0,  0.5, 
1,0  and  1.5  seconds,  are  all  one  ran  perform  within  this  programming 
restriction.  Nevertheless,  for  tone  durations  falling  bet\/een  these  four 
points,  a polynomial  spline  can  be  used  to  interpolate  the  four  average 
response  points  obtained  from  the  simulation  exercises.  At  a traffic 
density  specified  (for  instance,  for  the  N.Y.  HI  sector  function  the 
expected  traffic  density  Is  specified  as  33  aircraft  arrivals  per  hourl , a 
cubic  polynomial  curve  which  relates  the  average  response  (y)  with  the 
tone  duration  (x)  can  be  written  as 

2 3 

y = c^  + Cj^x  + + c^x  (3,3) 

where  c c^^,  c^  and  c are  constant  parameters  to  be  determined  from  a 
spline  fitting  procedure.  Since  now  we  have  four  observations  on  v and 
X and  the  same  number  of  unknown  parameters,  the  spline  equation  can  be 
completely  specified  by  solving  four  simultaneous  equations  each  express- 
ing a value  of  y by  its  corresponding  value  of  x,  in  the  polynomial 
function.  Values  other  than  the  multiples  of  0. 5-second  tone  duration  can 
then  substitute  the  variable  x in  Equation  (3.8),  and  the  values  of  y 
determined . 

On  the  other  hand,  for  a given  tone  duration,  but  varied  traffic 
density,  anotlier  spline  can  be  specified,  estimated  and  combined  with  the 
one  just  described.  Making  use  of  the  values  displayed  in  Cable  1.1,  in 
which  different  traffic  densities  were  imposed  upon  the  sector  function 
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simulations  with  a zero  tone  duration,  a spline  with  three  parameters 
as  expressed  below  is  adequate  for  relating  the  response  to  the  expected 
traffic  density; 

y = dj^z  + d^z^  + d^z^  (3.9) 

where  y is  the  average  response  results  from  a simulation  run  and  z is  the 
expected  traffic  density  (number  of  aircraft  arrivals/hour).  Again,  using 
three  pairs  of  moderate  values  of  y and  z,  from  table  3.1,  one  can 
completely  determine  the  parameters  in  the  polynomial  function  considered. 
For  values  of  z other  than  those  actually  tried  the  corresponding  values 
of  y can  be  interpolated,, 

Furthermore,  a convenient  method  useful  for  combining  the  two 
splines  discussed  above  to  form  a complete  surface  of  y for  various 
combinations  of  x and  z is  described  below.  Assume  that  the  values  of  y 
at  different  values  of  x are  proportionally  invariant  at  every  value  of  z, 
i.e.  the  ratio 


V 

c 


0 


.22 
+ — X 

"o 


, 3 3 

H-  X 

"o 


(3.10) 


remains  constant  for  all  values  of  z.  Let  y^  denote  the  value  of  y when 
the  tone  duration  Is  equal  to  zero,  i.e.  x = 0,  and  the  traffic  density  is 
equal  to  z^,  (y  = d^^z  +d2Z'^r)+d-^z^Q).  Then  the  value  jf  y for  a given  tone 
duration  x and  the  traffic  density  z^  can  be  computed  in  the  following 
way : 

c,  ^2  2 *^13 

y = ^0  c ^ ^ ^ ^ ) (3.11) 

^ 0 ‘'o  ^0 


Through  the  use  of  this  procedure  a complete  average  response  sur- 
face was  constructed  and  tabulated  for  each  of  the  three  responses:  aircraft 
loading,  channel  utilization  and  the  queuing  time  per  communication  trans- 
action. Totally  nine  sets  of  tiiree-tables  were  compiled  for  each  of  the 
nine  functions  investigated  (see  tables  3.20  - 3,46). 


For  practical  use  of  these  tables  two  special  notes  are  made 
below.  First,  in  occasional  cases  the  polynomial  functions  determined 
were  found  having  negative  slope  in  certain  range  of  the  values  of  the 
driving  variable.  To  remedy  this  difficulty  a horizontal  straight  line 
was  used  to  replace  the  portion  of  valley-like  curve.  Second,  the  average 
channel  utilization  rate  over  the  2-hour  observation  period  was  assumed 
upper  bounded  by  fractions  in  the  neighborhood  of  0.85,  due  to  the  common 
pauses  that  occur  during  the  ebb  and  flow  of  human  conversations.  For 
convenience,  all  values  which  were  originally  greater  than  0.85  were 
reduced  to  0.85,  thereby  to  indicate  the  existence  of  an  unstable  channel 
status  under  unduly  heavy  demands  for  conversations.  The  tabulated  values 
of  aircraft  loading  and  queuing  time  corresponding  to  a 0,85  channel  rate 
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are  in  fact  providing  the  very  lower  bounds  to  the  responses  expected 
from  a real  situation  represented  by  the  combination  of  tone  duration 
and  traffic  density. 

Several  extra  simulations  were  executed  to  provide  cFiecks  on  the 
accuracy  of  the  spline  projected  response  surface.  Results  are  summarised 
as  follows: 


Case  (1): 


Case  (ii) : 


Case  (ill) : 


New  York  HI  function 

0 = 20  aircraft  arrivals  per  hour,  tone  duration  = 1.0 
second 


Simulated  responses: 

(Ht.  C^.  Q^) 

= (3.88, 

0.48, 

7.31) 

Projected  responses: 

(Ht.  c^,  Q,) 

= (3.85, 

0.39, 

8.01) 

New  York  LT  Function 

0 = 20  aircraft  arrivals  per  hour; 

tone  duration 

= 1.0 

second 

Simulated  responses: 

(nt,  c^,  Q^) 

= (3.39, 

0.53, 

11.45) 

Projected  responses: 

(Ht,  C^,  Q^) 

= (2.50, 

0.56, 

13.24) 

New  York  AR  function 

0 = 10  aircraft  arrivals  per  hour; 

tone  duration 

= 1.0 

second 

Simulated  responses: 

(n^.  Ct,  Q^) 

= (1.25, 

0.30, 

3.83) 

Projected  responses: 

(iij.,  C^,  Qj.) 

= (1.90, 

0.32, 

9.78) 

Impressions  obtained  from  an  inspection  of  these  results  appear  to  be  In 
favor  of  an  acceptance  of  the  spline  projections  as  convenient , and  low 
cost,  approximations  to  the  responses  desired. 


1.3.3  Concluding  Remarks  on  the  Tone  Burst  Effect 

From  an  inspection  of  the  tables  provided,  it  is  found  that  tlie 
range  of  percentage  increases  in  average  channel  utilization  for  all  nine 
sector  functions  is  11  to  18  for  a tone  burst  length  of  .5  second  and  35 
to  54  for  1.5-second  bursts.  The  corresponding  ranges  of  percentage 
increases  in  average  queuing  time  are  28  to  80  and  138  to  475.  This  leads 
us  to  believe  that  implementation  of  the  proposed  FCC  regulation  may 
severely  impact  the  ability  of  the  agency  to  provide  adequate  air  traffic 
control  service  unless  some  way  is  found  to  prevent  tone  identification 
procedures  from  extending  the  duration  of  pilot-initiated  transmissions. 

A sample  set  of  tables  and  graphs  for  HI  sector  function,  from 
which  Cables  3.11  - 3.19  were  constructed  is  appended  to  tliis  chapter  (see 
tables3.47  - 3.50  and  figures  3.5  - 3.16). 
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3.4 


Effects  of  Changes  m Expected  Number  of  TR/CT  upon  Channel 
Responses 

From  studies  performed  to  date, It  was  found  that  among  the  major 
inputs  to  the  simulation  model  the  number  of  TR/CT  is  the  variable  most 
difficult  to  describe  using  conventional  statistical  models.  Historical 
histograms  were  therefore  suggested  to  be  used  as  inputs  to  the  simulation 
programs.  Furthermore,  nonhoraogeneity  of  this  variable  was  found  to  exist 
among  individual  sectors  within  a function.  All  these  facts  Indicate  the 
flexibility  or  adaptability  of  this  variable  to  varied  environments.  Ue 
suspect  that,  for  improving  the  ATC  performance  in  the  future,  the  number 
of  TR/CT  could  be  the  first  candidate  subjected  to  modifications. 

In  this  section,  without  specifying  detailed  changes  in  the 
comraunicational  technology,  we  hypothesize  a transfer  function,  which  trans- 
fers an  old  distribution  into  a new  distribution,  to  express  the  improve- 
ment due  to  a change  in  the  ATC  method  and/or  facilities.  The  effects  on 
various  responses  are  then  studied  through  simulation  experiments. 


3.4,1  Specification  and  Estimation  of  the  Transfer  Function 


Let  the  old  distribution  of  #TR/CT  be  denoted  by  p(j),  while  the 
new  one  by  p*(j);  in  both  cases  j = 1,2,3,  ...  . The  transfer  function 
w(x),  which  maps  p(j)  into  p*(j)  and  reduces  the  expected  value  to  a pro- 
portion, K,  of  the  original  level,  is  the  function  satisfying  the  following 
conditions : 


I P*(j) 

j = l 


I w(j)p(j)  = 1; 

i=i 


(il) 


) jp*(-1) 

J=1 


y j w(j)p(j) 

j=i 


X y j p(.i) 

j=i 


(3.12) 


The  first  condition  is  to  assure  that  p*(j)  is  a probability  function  and 
the  second  is  to  "shift"  the  expected  level  to  a new  point. 


A transfer  function  useful  for  /*TR/CT  is  the  one  which  can  be 
written  as: 


w(x) 


-x/b 


0 < X < "° 


(3.13) 


where  a and  b are  constants  to  be  determined  by  the  two  conditions  given 
in  Equation  (3.12).  This  function  means  that,  in  order  to  reduce  the 
average  level  of  the  variable  to  a given  point,  the  frequencies  at  the 
lower  range  of  the  variable  are  to  be  Increased  while  those  at  the  higher 
range  decreased.  The  shifts  of  the  weights  are  of  an  exponential  form. 
This  appears  to  be  an  efficient  way  of  Improving  the  responses. 
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3.4,2  Results  of  Simulation  Experiments 


To  Investigate  the  effects  of  this  shift  in  the  expected  value  of 
#TR/CT  upon  the  responses  of  concern,  a simulation  experiment  for  N.Y.  HI 
function  was  carried  out  using  distributions  determined  by  Equations  (3.13) 
and  (3.12).  For  a given  distribution  of  #TR/CT,  as  shown  in  figure  3.17, 
the  function  w(x)  and  the  new  distribution  p*(j)  are  shown  in  figures  3.18 
- 3.23  for  K = 0.9,  0.7,  and  0.5.  The  resulting  distributions  of  p*(j) 
were  then  used  as  Inputs  to  the  simulation  exercises,  at  a traffic  density 
equal  to  60  aircraft  arrivals  per  hour.  The  computer  outputs  for  the  four 
trials  are  provided  in  tables  3.51  - 3.54  and  figures  3,24  - 3.35  and 
appended  to  this  chapter.  A summary  table  (table  3.55)  is  given  for  the 
review  of  the  average  responses.  Significant  Improvements  in  the  responses 
are  observed  for  each  of  the  decrements  on  the  value  of  K.  The  explosive 
status  in  the  case  of  K = 1.0  has  been  transferred  to  a normal  operation  due 
to  a 30  per  cent  reduction  on  the  y*TR/CT  (for  R=0.7).  This  result  will  be 
useful  for  the  future  design  for  system  improvement. 

3. 5 Summary  and  Concluding  Remarks 

In  this  chapter  three  applications  of  the  simulation  model  have 
been  discussed;  namely,  the  estimation  of  sector  function  communications 
capacities,  the  evaluation  of  the  effects  of  tone  disturbance,  and  an 
illustration  of  Influence  of  changes  in  expected  number  of  TR/CT.  All 
three  applications  provide  valuable  information  for  the  future  design  for 
system  Improvement.  These  examples  also  suggest  the  general  appl icabi 1 1 ty 
of  the  simulation  model  to  various  studies. 
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w(x)  p{j) 


Figure  3.17  Original  Probability  Function  of  TR/CT 
(New  York  HI  Function) 


Figure  3.18  The  Weighting  Function  (k  = 0.9) 


Table  3.55 


Comparisons  of  Major  Responses  at  Various  Rates  of 
Reduction  on  #TR/CT 

Chi  Function,  N.Y,  ARTCC;  Expected  Traffic  Denslty=60A/HR) 


* 


«=  1. 

K = 0.9 

K = 0.7 

11 

o 

^t 

16.640 

13.447 

10.260 

9.343 

^t 

0.869 

0.824 

0.669 

0.475 

Qt 

89.104 

46.511 

12.222 

4.670 

(sec. ) 

Pr(5>l) 

0.3749 

0.090 

0.003 

0.000 

* ^ 

K = (New  average  //TR/CT)  / (Old  average  //TR/CT) 
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T&ble  3.3 


Pr(nj.>k)  at  Various  Traffic  Densities 


(LT  Sector  Function,  N.Y.  ARTCD 


k 

40A/HR 

45A/HR 

55A/HR 

10 

0.1219 

0.1359 

0.5225 

11 

.0653 

.0641 

.4121 

12 

.0316 

.0259 

.3082 

13 

.0138 

.0090 

.2178 

14 

.0054 

.0026 

.1450 

15 

.0019 

.0007 

.0906 

16 

.0006 

.0001 

.0531 

17 

.0002 

.0000 

.0291 

18 

.0000 

.0000 

.0149 

19 

.0000 

.0000 

.0071 

20 

.0000 

.0000 

.0032 

21 

.0000 

.0000 

.0013 

22 

.0000 

.0000 

.0005 

23 

.0000 

.0000 

.0002 

24 

.0000 

.0000 

.0001 

25 

.0000 

.0000 

.0000 

h 

1.2470 

1.1598 

1.2819 

h 

-0.3352 

-0.2603 

-0.3349 

-'2 

a 

0.9489 

0.7989 

0.8894 

y 

6.6280 

7.4010 

10.2040 

A/HR  = expected  number  of  aircraft  arrivals  per  hour 
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Table  3.4 


Pr(n^>k)  at  Various  Traffic  Densities 
(LE  Sector  li'unction,  N.Y.  ARTCC) 


k 

35A/HR 

45A/HR 

55A/HR 

70A/HR 

10 

0.0107 

0.2131 

0.4203 

0.8786 

11 

.0021 

.1278 

.3026 

.8382 

12 

.0003 

.0697 

.2026 

.7899 

13 

.0000 

.0344 

.1255 

.7341 

14 

.0000 

.0154 

.0717 

.6716 

15 

.0000 

.0062 

.0376 

.6039 

16 

.0000 

.0022 

.0181 

.5329 

17 

.0 

.0007 

.0080 

.4608 

18 

oO 

.0002 

.0032 

.3900 

19 

.0 

.0001 

.0012 

.3227 

20 

.0 

.0000 

,0004 

.2607 

21 

.0 

.0 

.0001 

.2055 

22 

.0 

.0 

.0 

.1579 

23 

.0 

.0 

.0 

.1182 

24 

.0 

.0 

.0 

.0861 

25 

.0 

.0 

.0 

.0611 

26 

.0 

.0 

.0 

.0421 

27 

.0 

.0 

.0 

.0282 

28 

.0 

.0 

.0 

.0184 

29 

.0 

.0 

.0 

.0116 

30 

.0 

.0 

.0 

.0071 

31 

.0 

.0 

.0 

.0043 

32 

.0 

.0 

.0 

.0025 

33 

.0 

.0 

.0 

.0014 

34 

oO 

.0 

.0 

.0008 

35 

.0 

.0 

.0 

.0004 

36 

.0 

.0 

.0 

.0002 

37 

oO 

.0 

.0 

.0001 

38 

oO 

oO 

.0 

.0 

A 

<(>1 

1.1417 

1.1516 

1.1946 

1.0966 

CM 

< -e 

-0.3038 

-0.2154 

-0.2637 

-0.1283 

a 

0.6882 

0.8377 

0.9945 

1.6658 

A 

y 

5.8520 

7.6690 

9.3620 

16.4580 
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Table  3.5  P (n  >k)  at  Various  Traffic  Densities 
r t 

(HI  Sector  Function,  N.Y.  ARTCC) 


k 

45A/HR 

50 A/ HR 

60A/HR 

10 

0.1934 

0.4380 

0.6733 

11 

0.0884 

0.3059 

0.6485 

12 

0.0332 

0.1953 

0.6232 

13 

0.0101 

0.1132 

0.5972 

14 

0.0025 

0.0592 

0.5709 

15 

0.0005 

0.0279 

0.5442 

16 

0.0001 

0.0118 

0.5173 

17 

0.0000 

0.0045 

0.4904 

18 

0.0000 

0.0015 

0.4634 

19 

0.0000 

0.0005 

0.4367 

20 

0.0000 

0.0001 

0.4102 

21 

0.0000 

0.0000 

0.3842 

22 

0.0000 

0.0000 

0.3586 

23 

0.0000 

0.0000 

0.3337 

24 

0.0 

0.0 

0.3095 

25 

0.0 

0.0 

0.2861 

26 

0.0 

0.0 

0.2635 

27 

0.0 

0.0 

0.2419 

28 

0.0 

0.0 

0.2213 

29 

0.0 

0.0 

0.2018 

30 

0.0 

0.0 

0.1833 

31 

0.0 

0.0 

0.1659 

32 

0.0 

0.0 

0.1496 

33 

0.0 

0.0 

0.1344 

34 

0.0 

0.0 

0.1204 

35 

0.0 

0.0 

0.1073 

36 

0.0 

0.0 

0.0954 

37 

0.0 

0.0 

0.0844 

38 

0.0 

0.0 

0.0744 

39 

0.0 

0.0 

0.0654 

40 

0.0 

0.0 

0.0572 

1.1522 

1.1942 

1.3581 

A 

1^2 

-0.3326 

-0.2906 

-0.3626 

a 

0.9532 

1.0674 

1.2398 

b 

8.2160 

9.5S50 

16.6400 
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Taale  3.6 


Pr(n^>k)  at  Various  Traffic  Densities 
(LC  Sector  Function,  N.Y.  Area) 


k 

85A/HR 

lOOA  /HR 

llOA/HR 

10 

0.6708 

0.9992 

0.9973 

11 

0.5693 

0.9985 

0.9965 

12 

0.4629 

0.9973 

0.9956 

13 

0.3591 

0.9951 

0.9945 

14 

0.2648 

0.9915 

0.9931 

15 

0.1850 

0.9857 

0.9915 

16 

0.1222 

0.9769 

0.9894 

17 

0.0761 

0.9638 

0.9870 

18 

0.0446 

0.9453 

0.9841 

19 

0.0246 

0.9199 

0.9807 

20 

0.0127 

0.8866 

0.9767 

21 

0.0062 

0.8443 

0.9720 

22 

0.0028 

0.7927 

0.9665 

23 

0.0012 

0.7323 

0.9602 

24 

0.0005 

0.6640 

0.9530 

25 

0.0002 

0.5898 

0.9447 

26 

0.0001 

0.5123 

0.9353 

27 

0.0000 

0.4343 

0.9247 

28 

0.0 

0.3588 

0.9128 

29 

0.0 

0.2885 

0.8996 

30 

0.0 

0.2254 

0.8850 

31 

0.0 

0.1710 

0.8b89 

32 

0.0 

0.1258 

0.8514 

33 

0.0 

0.0897 

0.8323 

34 

0.0 

0.0619 

0.8118 

35 

0.0 

0.0413 

0.7987 

36 

0.0 

0.0267 

0.7662 

37 

0.0 

0.0167 

0.7414 

38 

0.0 

0.0101 

0.7152 

39 

0.0 

0.0059 

0.6878 

40 

0.0 

0.0033 

0.6593 

A 

<I>1 

1.2824 

1.2551 

1.1921 

CM 

< -e- 

-0.3747 

-0.3035 

-0.2008 

a2 

a 

1.5563 

1.7169 

2.2231 

u 

11.6520 

26.1570 

45. 1950 
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Table  3.7 


Pr(n^>k)  at  Various  Trattic  Densities 
(LC  Sector  Function,  N.Y.  Area  ) 


k 

80A/HR 

90A/HR 

lOOA/HR 

10 

0.6048 

0.8072 

0.6560 

11 

0.4995 

0.7571 

0.6413 

12 

0.3943 

0.7007 

0.6263 

13 

0.2963 

0.6389 

0.6112 

14 

0.2113 

0.5733 

0.5958 

15 

0.1426 

0.5056 

0.5804 

16 

0.0908 

0.4377 

0.5648 

17 

0.0545 

0.3716 

0.5491 

18 

0.0308 

0.3092 

0.5333 

19 

0.0163 

0.2517 

0.5175 

20 

0.0081 

0.2005 

0.5016 

21 

0.0038 

0.1561 

0.4857 

22 

0.0017 

0.1187 

0.4699 

23 

0.0007 

0.0882 

0.4541 

24 

0.0003 

0.0639 

0.4384 

25 

0.0001 

0.0452 

0.4228 

26 

0.0 

0.0311 

0.4073 

27 

0.0 

0.0209 

0.3919 

28 

0.0 

0.0137 

0.3767 

29 

0.0 

0.0087 

0.3617 

30 

0.0 

0.0054 

0.3469 

31 

0.0 

0.0033 

0.3324 

32 

0.0 

0.0019 

0.3181 

33 

0.0 

0.0011 

0.3040 

34 

0.0 

0.0006 

0.2903 

35 

0.0 

0.0003 

0.2768 

36 

0.0 

0.0001 

0.2636 

37 

0.0 

0.0001 

0.2508 

38 

OoO 

0.0 

0.2383 

39 

0.0 

0.0 

0.2262 

40 

0.0 

0.0 

0.2144 

A 

1.2028 

1.2828 

1.1654 

<t>2 

-0.2747 

-0.3242 

-0. 1671 

^ X 

0 

1.4232 

1.8864 

1.7836 

u 

10.9950 

15.0790 

20.0970 

L. 


Table  3.8  Pr(n^>k)  at  Various  Traffic  Densities 
(DP  Sector  Function,  N.Y.  Area) 


k 

50A/HR 

55A/HR 

60A/HR 

70A/HR 

10 

0.0491 

0.0551 

0.2082 

0.7094 

11 

0.0187 

0.0173 

0.1091 

0.6321 

12 

0.0060 

0.0043 

0.0494 

0.5490 

13 

0.0017 

0.0008 

0.0193 

0.4636 

14 

0.0004 

0.0001 

0.0064 

0.3800 

15 

0.0001 

0.0000 

0.0018 

0.3016 

16 

0.0000 

0.0000 

0.0004 

0.2314 

17 

0.0000 

0.0000 

0.0001 

0.1715 

18 

0.0000 

0.0 

0.0 

0.1225 

19 

0.0 

0.0 

0.0 

0.0843 

20 

0.0 

0.0 

0.0 

0.0558 

21 

0.0 

0.0 

0.0 

0.0355 

22 

0.0 

0.0 

0.0 

0.0217 

23 

0.0 

0.0 

0.0 

0.0127 

24 

0.0 

0.0 

0.0 

0.0072 

25 

0.0 

0.0 

0.0 

0.0039 

26 

0.0 

0.0 

0.0 

0.0020 

27 

0.0 

0.0 

0.0 

0.0010 

28 

0,0 

0,0 

0,0 

0,0005 

29 

0.0 

0.0 

0.0 

0.0002 

30 

0.0 

0.0 

0.0 

0.0001 

31 

0.0 

0.0 

0.0 

0.0 

A 


4>i 

1.1326 

1.0023 

1.1631 

1.2635 

A 

(j)2 

-0.2413 

-0.1155 

-0.2796 

-0.3009 

^2 

0 

0.8606 

0.7180 

0.9139 

1.1208 

A 

M 

6.1370 

6.8950 

8.0590 

12.5730 
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Table  3.9 


Pr(rij^>k)  at  Various  Traffic  Densities 
fAD  Sector  Function.  N.Y.  Area  ^ 


k 

40A/HR 

45A/HR 

50A/HR 

60A/HR 

70A/HR 

10 

0.0001 

0.0000 

0.0374 

0.1380 

0.6211 

€1 

0.0000 

0.0000 

0.0135 

0.0715 

0.5538 

12 

0.0 

0.0 

0.0041 

0.0329 

0.4850 

.13 

0.0 

0.0 

0.0011 

0.0134 

0.4165 

14 

0.0 

0.0 

0.0002 

0.0048 

0.3506 

15 

0.0 

0.0 

0.0 

0.0015 

0.2889 

16 

0.0 

0.0 

0.0 

0.0004 

0.2328 

17 

0.0 

0.0 

0.0 

0.0001 

0.1833 

18 

0.0 

0.0 

0,0 

0.0 

0.1410 

19 

0.0 

0.0 

0.0 

0.0 

0.1059 

20 

0.0 

0.0 

0.0 

0.0 

0.0775 

21 

OoO 

0.0 

0.0 

0.0 

0.0554 

22 

0.0 

0.0 

0.0 

0.0 

0.0385 

23 

0.0 

0.0 

0.0 

0.0 

0.0261 

24 

0.0 

0.0 

0.0 

0.0 

0.0173 

25 

0.0 

0.0 

0.0 

0.0 

0.0111 

26 

0.0 

0.0 

0.0 

0.0 

0.0070 

27 

0.0 

0.0 

0.0 

0.0 

0.0042 

28 

0.0 

0.0 

0.0 

0.0 

0,0025 

29 

0.0 

0.0 

0.0 

0.0 

0.0014 

30 

OoO 

0.0 

0.0 

0.0 

0.0008 

31 

0.0 

0.0 

0.0 

0.0 

0.0004 

32 

0.0 

0.0 

0.0 

0.0 

0.0002 

33 

0,0 

0.0 

0.0 

0.0 

0.0001 

34 

0.0 

0.0 

0.0 

0.0 

0.0001 

35 

0,0 

0.0 

0.0 

0.0 

0.0 

A 

<0l 

1.0631 

1.1136 

1.2060 

1.2472 

1.3505 

A 

^2 

-0.2108 

-0.3297 

-0.3123 

-0.3715 

-0.3769 

-Z 

a 

0.5605 

0.4908 

0.7561 

1.0602 

1.0876 

A 

b 

4.0890 

4.6940 

5.8650 

7.0950 

11.7810 
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Table  3.10  Pr(n^>k)  at  Various  Traffic  Densities 
CAR  Sector  Function.  N.Y.  Area 


k 

35A/HR 

40A/HR 

10 

0o0991 

0.3952 

11 

0.0A96 

0.2896 

12 

0.0222 

0.1996 

13 

0.0088 

0.1289 

U 

0o0031 

0.0777 

15 

0.0010 

0.0437 

16 

0.0003 

0.0229 

17 

0.0001 

0.0111 

18 

0.0000 

0.0050 

19 

0.0000 

0.0021 

20 

0.0000 

0.0008 

21 

0.0000 

0.0003 

22 

0.0000 

0.0001 

23 

0.0000 

0.0000 

A 

<t>l 

1.324A 

1.3496 

A 

<t>2 

-0.4000 

-0.4073 

a2 

a 

0.6743 

0.8047 

A 

6.4430 

9.0790 
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Table  3.11 


Effects  of  Tone-Disturbance  upon  ATC  Connnunicatlons  Channels 
(N.Y,  ARTCC,  HI  Sector  Function) 

Traffic  Density  = 33  Aircraft  Arrivals/Hour 


Response 

Tone 

Duration  (seconds) 

0.0 

0.5 

1.0 

1.5 

n 

5.579 

5.982 

6.497 

6.933 

t 

c 

0.529 

0.611 

0.685 

0.755 

t 

Qt 

9.146 

16.448 

21.251 

29.576 

(sec) 

“ average  aircraft  loading  per  second 
c^  » average  channel  utilization 
“ average  queuing  time. 
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Table  3.12 


Effects  of  Tone-Disturbance  upon  ATC  Communications  Channels 
(N.Y.  ARTCC,  LE  Sector  Function) 


Traffic  Density  = 35  Aircraft  Arrivals/Hour 


Response  Tone  Duration  (seconds) 

0.0 0.5  1.0. 


“t 

5.852 

6.228 

7.313 

9.300 

^t 

0.591 

0.676 

0.766 

0.859 

sec) 

13.357 

21.016 

40.303 

76.829 

n 

t 

“ average  aircraft 

loading  per 

second 

c^  ■ average  channel  utilization 
~ average  queuing  time. 
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Table  3.13 


Effects  of  Tone-Disturbance  upon  ATC  Communications  Channels 
(N.Y.  ARTCC,  LT  Sector  Function) 


Traffic  Density  = 29  Aircraft  Arrivals/Hour 


Tone  Duration  (seconds) 


0.0 

0.5 

1.0 

1.5 

4.257 

4.692 

4.875 

5.320 

0.545 

0.616 

0.700 

0.770 

9.410 

15.683 

19.368 

28.148 

“ average  aircraft  loading  per  second 
c^  » average  channel  utilization 
■ average  queuing  time. 


Table  3.14 


Effects  of  Tone-Disturbance  upon  ATC  Connsunlcatlons  Channels 
(N.Y.  Area^GN  Sector  Function) 

Traffic  Density  = 43  Aircraft  Arrivals/Hour 


Response 

Tone  Duration  (seconds) 

0.0 

0.5 

1.0 

1.5 

°t 

2.374 

2,496 

2.608 

2.763 

0.289 

0.340 

0.388 

0.436 

(sec) 

3.830 

6.145 

10.091 

13,246 

n 

t 

” average  aircraft 

loading  per 

second 

c^  - average  channel  utilization 
■ average  queuing  time. 
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Table  3.15 


Effects  of  Tone-Disturbance  upon  ATC  Communications  Channels 
(N.Y.  Area  , LC  Sector  Function) 

Traffic  Density  • 42  Aircraft  Arrivals/Hour 


Response 


Tone  Duration  (seconds) 


0.0 


0.5 


1.0, 


1.5 


n 

t 

5.905 

6.275 

6 . 644 

7.291 

c 

t 

0.454 

0.537 

0.609 

0.700 

\ 

(sec. ) 

6.689 

8.540 

13.484 

22.656 

n ■ 

t 

average  aircraft 

loading  per 

second 

c ■ 

t 

average  channel 

utilization 

^t  ■ 

average  queuing 

time. 
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Table  3.1f> 


Effects  of  Tone-Disturbance  upon  ATC  Communications  Cliannels 
(N.Y.  Area  , LG  Sector  Function) 

Traffic  Density  = 25  Aircraft  Arrivals/Hour 

Response  Tone  Duration  (seconds) 

0.0 O^J T;0_ U5_ 

n 2.990  2.998  3.2AO  3.606 


0.306 

0.354 

0.405 

0.453 

7.232 

10.944 

20.381 

24.628 

(sec . ) 


n^  • average  aircraft  loading  per  second 

c^  ■ average  channel  utilization 

“ average  queuing  time. 
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Table  3.17 


Effects  of  Tone-Disturbance  upon  ATC  Communications  Channels 
(N.Y.  Area  , DP  Sector  Function) 


Traffic  Density  • 25  Aircraft  Arrivals/Hour 


Response  Tone  Duration  (seconds) 

0.0  0.5  1.0  1.5 


n 

t 

3.295 

3.565 

3.907 

4.191 

c 

t 

0.A69 

0.535 

0.591 

0.660 

\ 

(sec . ) 

10.288 

13.884 

19.099 

28.000 

n 

* average  aircraft 

loading  per 

second 

c^  « average  channel  utilization 
■ average  queuing  time. 
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Table  3.18 


Effects  of  Tone-Disturbance  upon  ATC  Communications  Channels 
(N.Y.  Area  , aD  Sector  Function) 

Traffic  Density  = 16  Aircraft  Arrivals/Hour 


Response 


Tone  Duration  (seconds) 


0.0 

0.5 

1.0 

1.5 

n 

t 

1.592 

1.602 

1.741 

1.813 

0.240 

0.271 

0.304 

0.337 

^t 

(sec.  ) 

4.779 

6.398 

9.981 

11.390 

n - 

t 

average  aircraft 

loading  per 

second 

^t  " 

average  channel 

utilization 

^t  - 

average  queuing 

time. 
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icy  tn 


Table  3.19 


Effects  of  Tone-Disturbance  upon  ATC  Communications  Channels 
(N.Y.  Area  , AR  Sector  Function) 

Traffic  Density  = 22  Aircraft  Arrivals/Hour 

Tone  Duration  (seconds) 

0.0 0^5 UO. 1-5 

2.892  3.151  3.469  3.780 

0.527  0.585  0.644  0.712 

11.232  15.087  21.512  27.277 


Response 


, Ic.) 


ii^  *•  average  aircraft  loading  per  second 
«■  average  channel  utilization 
» average  queuing  time. 
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Table  3.47  Sttmmary  Statistics  of  Simulations  for  the  Study 


NOHBEP  0?  ZE33  ENTRIES  (NON-HAITINl)  = 1S3 

PERCENT  0?  ZERO  ENTRIES  = 44. 

A?BPA3R  NiJNSEo  OF  AIPCRAFT  TK  OfJFfJF  = .434 

PAJCTFa*  rjHSEP  CF  AIRCRAFT  T*»  OUBhE  = 


Function 


Figure  3.8  Simulated  Aircraft  Loading  Time  Series  for  N.Y.  HI  Function 
33  A/H  and  0.5  Second  Tone  IXiration 


Figure  3.10  Simulated  Queue  Length  Time  Series  for  N.Y.  HI  Function 


Disturbance  (Duration  =1.0  Seconds) 
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Figure  3.12  Simulated  Channel  Utilization  Time  Series  for  N.Y.  HI  Function 


Figure  3.13  Simulated  Queue  I.ength  Time  Series  for  N.Y.  HI  Function 


Table  3.50  Summary  Statistics  of  Simulations  for  the  Study 
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Figure  3.24  Simulated  Aircraft  I.oadlng  Time  Series  for  N.Y.  HI  Function  at  60  A/H  and  K=l. 


Figure  3.25  Simulated  Channel  Utilization  Time  Series  for  \.Y.  HI  Function  at  60  A/H  and  K= 
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Figure  3.28  Simulated  Channel  Utilization  Time  Series  for  N.Y.  HI  Function  at  60  A/H  and  K=0.9 
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Table  3.53  Suimnary  Statistics  of  Simulations  for  the  Changes 
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Figure  3.30  Simulated  Aircraft  l.oadlng  Time  Series  for  N.Y.  HI  Function  at  60  A/H  and  K=0.  7 


Figure  3.31  Simulated  Channel  Tltlllgatlon  Time  Series  for  N.Y  HI  Function  at  60  A/H  and 
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Figure  3.33  Simulated  Aircraft  Loading  Time  Series  for  N.Y.  HI  Function  at  60  A/H  and 


Figure  3.34  SinmlafeH  Channel  Iltlllzatlon  ’’‘Ime  Series  for 


CHAPTER  4 


NETWORK  MODELING  - NETWORK  STRUCTURE,  AIRCRAFT  ARRIVAL 
STREAMS,  AND  SECTOR  SEQUENCES 


4.1  Introduction 


4.1.1  Motivation 

The  ATC  voice  communications  study  presented  in  previous  Princeton 
reports  has  concentrated  on  the  analysis  and  simulation  of  air-ground 
communications  models  for  a single  sector  or  sector  function.  The  simulated 
sector  is  excited  by  a stream  of  traffic  which  is  generated  without  regard 
to  sector  flows.  Valuable  insights  regarding  various  responses  such  as 
queuing  times,  channel  utilization,  etc.,  have  been  gained  by  exercising  the 
model  in  many  different  ways,  but  little  consideration  has  yet  been  given 
to  the  Important  transfer  phenomena  of  some  significant  events,  such  as 
traffic  congestion,  in  a sequence  of  adjacent  sectors.  Indeed,  a systematic 
study  of  a group  of  sectors  as  a network  remains  to  be  done. 

This  chapter  presents  the  results  of  research  conducted  into  various 
aspects  of  a network  composed  from  a group  of  ATC  enroute  sectors,  based  upon 
the  New  York  ARTCC  data  file.  This  chapter  could  form  the  baais  for  a more 
comprehensive  study  aimed  at  providing  valuable  information  regarding  the 
interrelationships  of  the  ATC  activities  among  the  sectors  in  a network  which 
might  lead  to  the  suggestions  of  methods  useful  for  the  Improvement  of  the 
communications  efficiency  of  the  network  as  a whole. 

The  selection  of  the  enroute  sectors  was  suggested  by  a number  of 
considerations.  As  a group,  these  sectors  appear  to  be  the  most  homogeneous, 
in  terras  of  the  probability  distributions  used  to  describe  the  major 
communication  variables,  and  in  the  type  of  traffic  pattern  handled.  The 
traffic  flows  among  these  sectors  are  the  most  diffuse  when  compared  with  the 
sector  groups,  due  to  the  large  number  of  paths  passing  through  each  sector. 
Further,  enroute  sectors  are  felt  to  be  analogous  to  oceanic  sectors  where 
network  analysis  and  simulations  may  have  applications.  The  ultimate 
goal  of  this  study  is  to  Integrate  the  single-sector  simulation  into  a 
general  network  model  which  will  be  useful  for  the  purposesjust  mentioned. 

4.1.2  Scope  of  this  Report 

The  three  facets  of  the  multiple-sector  network  on  which  this  first 
report  will  concentrate  relate  to  the  mechanics  of  flow  among  the  enroute 
sectors.  The  first  section  develops  a formulation  for  a general  ATC  network 
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using  terminology  from  graph  theory  which  has  commonly  been  applied  to 
ground  transportation  systems.  The  next  section  discusses  the  specifications 
of  traffic  streams,  considering  sources  and  sinks  for  traffic  in  the  network 
and  characterizing  movement  of  aircraft  through  sequences  of  sectors.  The 
final  section  applies  these  techniques  to  a network  of  enroute  sectors  from 
the  New  York  data  file. 


4.2  Network  Structure 


4.2.1  Attributes  of  the  Network 

A general  ATC  network  consists  of  a set  of  sectors  , k = l , 2, . . . ,m) , 
each  controlled  by  a single  control  position  with  its  assignea  radio  fre- 
quency. Aircraft  move  among  the  various  sectors  as  permitted  by  geographical 
proximity  or  functional  relationship.  Adopting  the  notation  used  by  Ford  and 
Fulkerson  (1962),  Hu  (1969),  Potts  and  Oliver  (1972),  and  others,  the  network 
is  defined  by  G=[N;A]  where  N is  the  set  of  sectors  or  "nodes"  and  A is  a 
set  of  "arcs"  representing  feasible  transition  paths  between  adjacent  sectors. 

It  is  convenient  to  define  r super-source  s and  super-sink  t to 
represent  the  points  of  entry  and  exit,  respectively,  for  all  traffic  in  rlie 
network.  Then  if  (X,Y)  is  the  set  of  all  arcs  from  nodes  in  set  X to  nodes 
in  set  Y, 

A = (s,  N)  U (N,  N)  U (N,  t) 
where  (J  is  the  symbol  for  set  union. 


2 For  a network  containing  m sectors,  A could  contain  as  many  as 

(m  +2m)  directed  arcs,  each  representing  flow  in  one  direction  between  a pair 
of  nodes.  In  most  networks,  however,  the  number  of  arcs  is  much  less  than 
the  maximum  possible  since  many  pairs  of  sectors  are  not  physically  or 
functionally  adjacent.  To  specify  which  arcs  are  present  in  a network,  a 
node-node  incidence  matrix  H of  dimension  (m  x m)  may  be  defined  where  element 


d 


ij 


1,  if  flow  is  possible  from  node  i to  node  j ; 
0,  otherwise  . 


For  convenience,  it  is  assumed  that  all  2m  arcs  in  (s,  N) __and..(-Nv"T)''exTst , 
which  does  not  restrict  the  generality  pf,.i-he-ftrrmulation  since  the  proba- 
bility of  using  infeasible  arCs'Tor  the  capacity  of  such  arcs)  can  be  set 
to  zero.  ■ 

4.2.2  Defining  Flow  on  the  Network 

Flow  over  the  network  is  defined  on  the  arcs  in  A.  Let  f be  a 
real-valued  function  on  A with  dimensions  of  aircraft  per  unit  time,  hence- 
forth called  the  "traffic  flow  rate."  Then  if  the  total  rate  of  flow  througli 
the  network  is  F, 
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f (s,N) 


(4.1) 


F = r f(s,N,^) 

Equivalently, 

F = E f(N  ,t) 
Nk^N 


(4.2) 


Actually,  (4.1)  and  (4.2)  are  special  cases  of  Klrchhoff’s  flow  conservation 
law  which  states  that  In  equilibrium  the  total  flow  directed  Into  a node 
must  equal  the  total  flow  directed  outward.  Thus 

F,  If  x=s 

I f(x,Nk)  - „ f(Nk,x)  = 0,  If  x^s,t  (4.3) 

NktN  k -F,  If  x=t. 


In  determining  flow  rates  In  an  actual  network  from  data  collected 
over  a finite  sample  period,  intuitive  estimators  may  give  rise  to  estimates 
which  do  not  satisfy  (4,3)  and  for  which  (4,1)  and  (4.2)  are  not  equivalent. 

A common  reason  for  this  is  unequal  numbers  of  aircraft  In  each  node  (sector) 
at  the  beginning  and  end  of  the  sample  period.  In  order  to  arrive  at  estimates 
which  satisfy  (4,3),  It  is  necessary  to  isolate  separate  streams  of  traffic 
and  to  estimate  the  rate  of  flow  in  each  stream  individually,  or  to  character- 
ize the  traffic  flow  pattern  more  economical Iv  while  generating  consistent 
estimates.  Sections  4.3  and  4.4  discuss  this  problem  In  more  detail. 


4.3  Network  Arrivals 


4.3.1  Arrival  Sources 

The  arrivals  of  aircraft  to  the  network  0 are  considered  to  he 
generated  at  the  super-source,  s.- ■ From- •s',' 'entry  to  the  network  occurs  over 
one  of  the,^m. arcs' 'fh' (s , N),  In  reality,  the  flow  through  s is  a super- 
p'osftion  of  traffic  streams  from  many  sources,  such  as  airports  in  the 
region  or  other  ATC  centers.  Similarly,  the  flow  through  the  super-sink  t 
Is  a superposition  of  traffic  streams  hound  for  many  destinations. 

Denote  the  set  of  sources  feeding  traffic  to  the  s>iper-source  by 
fe^,  1=1,2,  ...,  m } and  the  set  of  sinks  from  the  super-sink  by 
(tj,  j=l,2,  ...,  m'^}.  Then  all  flow  through  the  network  originates  at  some 
s.,  passes  through  the  super-source  s Into  the  network,  then  exits  tlirough 
t to  some  destination  t.. Figure  4.1  Is  a schematic  diagram  of  such  a network 
with  four  sources  and  fiur  sinks, 

4.3.2  Determining  Flow  Rates  on  (s,  N) 

Assume  that  the  rate  of  traffic  generated  at  each  source  Sj^  is 
and  that  the  rate  of  traffic  attracted  to  each  sink  t^  is  ;i  ^ . Then  an 
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Figure  4.1  ATC  Network  Diagram 


application  of  the  flow  conservation  equation  requires  that 


F 


m 

s 


i=l 


A. 

1 


t 

J=1 


(4.4) 


where  F is  defined  to  be  the  rate  of  flow  over  the  entire  network,  (or 
equivalently  the  rate  of  flow  from  s to  t).  From  (4.4),  it  is  evident 
that  the  rate  of  flow  over  the  arcs  in  (s,  N)  must  equal  both  the  sum 
of  the  generation  rates  A^  and  the  sum  of  the  attraction  rates  pj. 

Given  A^  and  it  is  necessary  to  determine  each  f(s,  Nj^)  such 
that  the  above  conservation  equations  are  satisfied.  To  formulate  the 
problem  let 

p = Pr  (an  aircraft  enters  the  network  over 

arc  (s,  N )) 
k 

where  Pr  represents  the  probability  and  let 

Pj^l^  = Pr  (an  aircraft  enters  the  network  over 

arc  (s,  Nj^)  given  it  was  generated  at 
source  1) . 


Then 


f(s,  N^) 


ni 

s 


(4.5) 


If  the  selection  of  entry  arc  from  s does  not  depend 
which  an  aircraft  was  generated,  then  p,  I . = P,  and 

K.  I 1 cC 

m 

V = Pk 

X=1 

or  equivalently 

f(s,  Nj^)  = F. 

The  simple  form  of  (4.7)  will  not  be  appropriate  for  most  net- 
works representing  physical  systems  since  the  point  of  entry  to  the  network 
will  usually  depend  on  the  source  of  the  arrival.  For  an  ATC  network,  man> 
Pkli  equal  to  zero  for  positive  pj^  simply  from  geographical  con- 

siderations. For  example,  a departure  from  an  airport  in  the  eastern  part 
of  the  region  can  hardly  enter  the  network  at  any  sector  in  the  western 
part.  Thus  (4.5)  is  normally  required  for  specification  of  flow  rates  over 
the  arcs  in  (s,  N) , and  all  must  be  determined. 


On  the  source  at 

(4.6) 


4-5 


A consideration  which  further  complicates  the  situation  is  that 
the  Aj^  and  p . are  often  functions  of  each  other.  In  an  analysis  of  the 
network,  the  specification  of  sink  attraction  rates  may  be  more  appro- 
priate to  the  problem  at  hand  than  specification  of  source  generation 
rates.  It  is  even  possible,  as  will  be  shown  in  a later  section,  that 
some  combination  of  the  A^'s  and  U:'s  will  be  given,  requiring  that  the 
others  be  determined.  A discussion  of  techniques  for  solving  this  prob- 
lem must  be  postponed  until  the  specification  of  sector  transition  proba- 
bilities is  given  in  section  4.4. 

4.3.3  Estimation  of  Arrival  Rate  Parameters 

Assume  that  a given  traffic  stream  entering  the  super-source  s 
has  flow  rate  A.  Let  {A(t),  t > 0}  be  the  corresponding  counting  process 
for  the  number  of  aircraft  in  that  stream  passing  a given  point,  such  as 
the  point  of  entry  to  s.  Such  a point  may  have  a physical  meaning  or 
be  simply  an  abstract  reference.  Then  the  structure  of  {A(t)}  will  be 
Important  in  determining  an  appropriate  estimate  for  A.  The  most  fre- 
quent assumption,  which  has  been  applied  successfully  in  the  past  for 
streams  of  aircraft,  is  that  of  a Poisson  process.  Parzen  (1962) 
defines  a Poisson  process  as  a counting  process  such  that 

(i)  {A(t),  t>0}  has  stationary  independent  increments;  (4.S) 

(ii)  for  any  Interval  (t]^,t2),  A(t2)-A(t2)  has  a Poisson 

distribution  with  mean  A(t2-t,J  where  A is  the  rate  of 
occurrence  per  unit  time  or  events  being  counted. 


The  Poisson  process  arises  naturally  in  circumstances  where  events 
occur  independently  and  has  several  interesting  properties.  Among  the 
most  commonly  used  is  that  the  length  of  time  between  k events  has  a 
gamma  distribution  with  parameters  (k,A).  In  particular,  the  time  be- 
tween consecutive  arrivals  (the  "interarrival  time")  follows  the  exponen- 
tial distribution  with  probability  density  function 


f(x) 


Ae 


-Ax 


X > 0 


(4.^) 


A second  Interesting  property  is  that  the  superposition  of  n independent 
Poisson  processes  with  rate  parameters  A.,  i = l,2,...,n,  is  also  a Poisson 
process  with  rate  parameter 

n 

'''  = I . (4.10) 

i=l 


Consider  the 
traffic  generated  at 
parameter  A^  and  all 
the  stream  through  s 


stream  of  traffic  through  the  super-source  s.  If 
each  of  the  sources  s^  is  a Poisson  process  with  rate 
the  traffic  streams  are  generated  independc’nt Iv, then 
is  also  Poisson  with  rate  parameter 


4-u 


(A. 11) 


A 

s 


I >• 

1=1  " • 


Further,  the  joint  distribution  of  arrivals  per  unit  time  at  s from  each 

of  the  sources  s.  follows  the  multinomial  distribution 
1 


X ! 

N - S 

[111 

^1 

^2 

^2 

'1 

m 

s 

> X.-v  1 , , 

2 m X,  J x^ ! . . .X  ! 

s 12  m 

A 

• S'* 

A 

^ s ' 

where  x 

s 


m 

s 

y X . is  the  total 
1=1  ^ 


= 0,1,2, . . . 


(4.12) 


number  of  arrivals  at  s In  one  unit  of  time. 


Assume  now  that  traffic  In  the  network  Is  observed  during  a finite 
period  and  that  estimates  are  desired  for  each  of  the  arc  flow  rates  f (s,  Nj^) . 
Then  from  14.5),  we  require  estimates  of  the  source  generation  rate  para- 
meters {A^,  1 =l,2,...,m  } and  of  the  conditional  probabilities  p.  ' . • 

For  convenience  the  A.'s^can  be  represented  by  a vector  of  parameters  A 

m xl 

and  the  conditional  probabilities  by  a transition  matrix  P . All 


parameters  are  assumed  to  be  constant  over  time. 


m xm 
s 


Suppose  that  over  the  duration  of  the  sample  period  we  observe  R 
arrivals  through  s.  If  the  sample  period  consists  of  T discrete  intervals 
of  time,  then  for  each  arrival  with  entry  time  t , the  data  can  be 
represented  by  (x^(t),  t=l,2,...,T)  where  ^ 

s.  if  t<tr  and  the  r*"^  aircraft 
was  generated  at  s^^ 

if  and  the  r*"*^  aircraft 

entered  the  network  over  arc 

(s,H^). 

VJe  will  show  that  the  likelihood  function  of  the  parameter  set  0 = (A,  P) 
given  a sample  X can  be  factored  into  two  terms,  one  involving  ^ only  and 
the  other  involving  P only.  The  proof  follows  that  of  Duncan  and  Lin  (1972) 
for  Markov  chains  with  stochastic  entry. 


For  each  observed  arrival,  no  information  is  lost  if  x (t)  is 
replaced  by  the  triplet  (s.,  Nj^,  t ).  Then  the  sampling  distribution  can 
be  written  as  ^ 


f0(s<)  = ^(s,  N,  t)  = f^(N|s,t)  fg(s,t) 
= fp(N|s,t)  f^(s,t)  . 


(4.11) 
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Thus  the  likelihood  function  of  the  parameters  given  x can  be  written 


L (0)  = L , . (P)  • L (A' 

X N|S,t  . s,t  ~ 


(4.14) 


If  the  transition  probabilities  Pj^i^  are  constant  over  time,  then  (4.14) 
reduces  to  ' 


(4.15). 


The  consequence  of  (4.15)  is  that  the  overall  likelihood  function 
for  6 Is  maximized  when  maximum  likehood  estimates  are  obtained  separately 
for  P and  A.  Further,  the  estimates  of  P need  involve  only  the  number  of 
aircraft  entering  each  sector  from  each  source,  and  the  estimates  of  A 
require  only  knowledge  of  the  time  and  soiarce  of  each  entry  to  s. 


Suppose  that  in  a sample  period  wh  observed  a source-node  entry 

count  matrix  C with  elements  defined  by 
m xm 
s 

c.,  = # of  aircraft  generated  at  source  s.  which 

lit  , N r 

entered  sector  . 

k 

Then  the  total  number  of  aircraft  generated  at  source  s.  is  given  by 

m 


k=l 


‘^ik 


and  the  number  of  aircraft  entering  the  network  over  arc  (s,  N^^)  is  given  by 

s 

".k=  ^ik- 
1=1 

Since  the  arrival  streams  from  the  Individual  sources  are  assumed  to  be 
independent  Poisson  processes,  the  maximum  likelihood  estimates  of  the 
source  generation  rate  parameters  are  given  by 

c . 


'^i  T 


(4.16) 

where  T is  the  number  of  time  units  in  the  sample  period,  a well-known  result. 


Now  since  the  conditional  probabilities  P|^  | ^ have  been  assumed  to  be 

constant  over  time,  then  given  c.  , the  individual  transition  counts  c^^  have 
a multinomial  distribution.  Thus  the  maximum  likelihood  estimates  are  simply 


'ik 


^k  i 


c . 

1 . 


(4.171 


The  estimate  of  tlie  overall  network  flow  rate  F can  be  obtained  using 
(4.4)  and  (4.16),  and  is 


4-8 


which  could  be  obtained  directly  by  realizing  that  the  flow  through  the 
super-source  s is  Poisson  with  rate  given  by  (a.H),  The  estimates  of  the 
flow  rates  over  the  arcs  in  (s,  N) , using  (4.5)  and  ("4. 17),  are 


^ nig  Q 

t(s,  N ) = I . (4.1U) 

i=l  ^ 

Since  the  i are  constant  over  time,  these  traffic  streams  are  also  Poisson. 
The  estimate  of  the  flow  rate  through  s is  given  by 

m 

A s 

A = y A.  = F . (4.20) 

s . , 1 

1=1 

4.3.4  Summary 

The  results  of  this  section  show  that  under  two  basic  assump- 
tions, the  arrival  rates  at  each  sector  of  a network  can  be  determined  from 
the  arrival  rates  at  individual  sources.  The  maximum  likelihood  estimates 
of  traffic  flow  rates  are  easily  obtained  from  observing  entries  to  the 
network  during  a sample  period.  The  basic  assumptions  required  are: 

(i)  that  the  generation  of  aircraft  at  the  sources  are  inde- 
pendent Poisson  processes; 

(ii)  that  the  probability  of  an  arrival  from  a given  source  en- 
tering a specific  sector  is  constant  over  time. 

When  simulating  the  system,  the  number  of  aircraft  generated  at  s 
per  unit  time  will  have  the  multinomial  distribution  given  by  (4.12)  with 
whatever  generation  rates  for  the  sources  are  desired.  Thus  only  a 
single  generator  at  s is  required,  not  one  for  each  source.  Using  the  esti- 
mates of  Pk|i  from  (4.17),  aircraft  can  then  be  assigned  to  the  proper 
arcs  in  (s,  N)  in  a random  fashion.  The  movement  of  aircraft  after  they 
enter  the  initial  sector  is  the  subject  of  the  next  section. 


4.4  Sector  Sequences 
4.4.1  Definition 

As  aircraft  move  through  an  ATC  network,  they  pass  from  sector  to 
sector  (or  from  node  to  node)  in  sequences  determined  by  their  origin  and 
destination.  In  a network  of  many  sectors,  the  number  of  possible  sequences 
is  enormous,  making  the  specification  of  the  relative  frequencies  of  all 
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sequences  prohibitive.  In  order  to  reduce  the  complexity  of  the  problem 
while  still  maintaining  the  general  patterns  of  network  flow,  an  approach 
based  on  Markov  chains  will  be  presented. 

To  state  the  problem  formally,  let  {s^(k),  k = 0,  1,  2,  ...}  be  the 
sequence  of  sectors  through  which  the  n*"^  aircraft  in  some  index  set  passes, 
where 

i,  if  the  n^^  aircraft  is  a departure 
from  the  i^^  source 

r,  if  the  k^^  sector  entered  by  the 
n aircraft  is  sector  r 

1 < k < m 
— — n 

j,  if  the  n*"^  aircraft  is  an  arrival  at 

, .th  . , 
the  j sink 

k > m 

n 

where  m^  is  the  number  of  network  sectors  in  the  sequence  for  the  n*"*^  aircraft. 
Then  the  m nodes,  mg  sources,  and  m^  sinks  can  be  considered  to  be 
states  of  a Markov  chain  of  unknown  order  from  which  {s  (•)}  is  one 
possible  observation. 

It  will  be  recognized  that  the  nodes  and  sources  are  transient 
states,  while  the  sinks  are  absorbing.  Further,  the  nodes  form  a communi- 
cating class,  accessible  from  the  sources,  but  the  sources  (non-return  states) 
are  not  accessible  from  any  states  in  the  chain. 

The  characterization  of  the  sequences  will  thus  involve  state  trans- 
ition matrices  of  very  special  form.  While  all  sequences  begin  in  one  of 
the  source  states,  the  probability  of  ever  returning  to  those  states  must 
be  zero.  In  describing  the  sequences,  we  state  first  the  initial  dis- 
tribution of  s (0)  and  then  discuss  the  state  transition  probabilities, 
n 


s (0) 
n 


s (k) 
n 


Sn(k) 


4.4.2  Initial  Distribution  of  s (0) 


The 


int lal 


distribution 


of 


s 

n 


n 

(0)  has  parameter  set 

^™s  1 


where 


4-m 


(A.21) 


W 


Pr(Sn(0)  = i)  = ^ 

m s 

s 

and  \ X.  = X . Thus  the  relative  generation  rates  at  the  sources  deter- 
i=l 

mine  the  nrobabilitv  distribution  ^or  ^ natural  way.  Further, 

the  properties  of  the  Poisson  model  for  source  generation  streams  ensure 
that  successive  observations  of  s (0)  for  n=l,2,...,  will  be  independent. 


Transition  Probabilities 

If  the  sector  sequences  {Sj^(’)}  can  be  regarded  as  realiza- 
tions of  a Markov  chain  of  order  q,  then  the  distribution  of  &n(l<-) 
depends  on  the  past  history  of  the  sequence  only  through  s (k-1) , 

s (k-2) , and  s (k-q).  More  explicitly,  let 

n n 

Pj  = Pr(s^(k)  = j) 

= Pr(s^(k)  = jls^(k-l)  = i) 

and  p^^^  = Pr(s^(k)  = jls^(k-l)  = i,  s^(k-2)  = h)  . 

% 

Then  If  the  sequences  are  zero-order  Markov  chains. 


(4.22) 

(4.23) 

(4.24) 


For 

For 


The 


first  order 


second-order 


extension  to 


’’i  ■ ’’ij  ■ »hi., 
Markov  chains, 

"i  * hi  ■ fhij 

chains, 

'’i  ^.i  ^ ’’hi.1 

higher  orders  is 


direct. 


(4.25) 

(4.26) 

(4.27) 


In  studying  sequences  of  sectors,  it  is  therefore  necessary  to 
determine  both  the  order  of  the  chain  and  all  relevant  transition  probabili- 
ties. This  is  most  easily  handled  by  defining  a series  of  transition 
matrices  p^  p*'  , ...  where  p'^  , the  "q-step  transition  matrix" 

has  element 


Pij  = Pr(s^^(k)  = j|s^(k-q)  = i) 

q = 1,  2,  . 

For  a zero-order  Markov  chain,  (4.25)  can  be  expressed  by 
P^''^  = [p  1 p = p . 

Ij  1.1  1 


For  first-order  Markov  chains,  a consequeni'e  of  (4.26)  is  that 
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p<"> . 


(4.28) 


(4.30) 


i 

'd 


For  second-order  and  higher  Markov  chains,  no  such  simple  representation 
is  possible. 

For  an  ATC  network,  it  is  instructive  to  determine  what  elements 
of  the  one-step  transition  matrix  p^  '^are  non-zero.  Figure  U ,2  shows  a 
sample  network  with  two  sources,  two  sinks,  and  five  interior  nodes.  The 
location  of  all  non-zero  elements  in  the  corresponding  one-step  transition 
matrix  is  shown  in  figure  4. 3.  The  matrix  has  been  partitioned  into 
nine  sections  for  ease  of  interpretation.  The  following  observations  are 
important ; 

(1)  Since  all  flow  from  the  sources  and  to  the  sinks  is 
uni-directional,  all  elements  in  sections  (i) , (iv) , 

(vii) , (viii) , and  (ix)  are  zero.  The  same  is  true 
for  all  higher-step  transition  matrices. 

(2)  Since  all  flow  in  the  network  passes  through  at 
least  one  interior  sector,  all  elements  in  section 
(iii)  are  zero.  This  is  not  true  in  higher  step 
transition  matrices.  In  fact,  for  large  q,  the  only 
significant  non-zero  elements  will  be  located  in 
section  (iii) . 

(3)  Since  transitions  fnim  any  state  to  itself  are 
not  allowed,  all  elements  on  the  main  diagonal 
are  zero. 

(4)  Since  the  interior  nodes  of  the  network  allow  J 

flow  in  both  directions  between  adjacent  pairs,  j 

the  non-zero  elements  in  section  (v)  are  sym-  j 

metrically  located  about  the  main  diagonal.  ! 

This  would  not  be  true  in  general  for  all 

ATC  networks. 

The  one-step  transition  matrix  is  closely  related  to  the  node- 
node incidence  matrix  D defined  in  section  4 . 2 . 1 . The  number  of  non- 
zero elements  in  the  matrix  is  simply  the  number  of  arcs  in  the  network 
and  the  location  of  the  non-zero  elements  in  section  (v)  is  equivalent  to 
the  location  of  I's  in  D. 

Higher  step  transition  matrices  determine  the  accessibility  of  nodes 
in  the  network  from  each  other.  Since  all  flow  originates  at  a source  and 
ends  at  a sink,  the  only  elements  of  the  matrix  which  do  not  converge  to  j 

zero  as  q becomes  large  are  those  in  section  (iii).  Let  the  limiting  values 
of  the  elements  in  section  (iii)  be  given  by 

(4.31)  i 

i = 1,  2 m^  j 

j — 2,  ...,  m • 

t I 
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Then  the  sink  attraction  rates 
by 

m 

bj=  l' 

i=l 


are  related  to  the  source  generation  rates 


A.  e.  . . (4.  ’2) 

1 ij 


To  prove  (4.32),  sum  both  sides  over  all  sinks,  giving 

m mm 

t t s 


I 

j = l 


I 

1=1 


i=l 


A. 
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e . . . 

ij 


Exchanging  summations  on  the  right. 


in 

r 

j=i 

m 

= r 
1=1 

m 

^i 

i-l 

e . . . 

m 

t 

Now  y e. . = 1 for 

all 

sources  i.  so 

that 

j=l 

m 

t 

I 

P 

m 

s 

= I 

A.  . 

j=i 

j 

i=l 

1 

But  from  (4.4),  the  sum  of  sink  attraction  rates  and  the  sum  of  source 
generation  rates  are  both  equal  to  the  total  flow  rate  over  the  network, 
so  that  (4.32)  must  be  true. 


The  importance  of  (4.32)  is  that  if  the  e^j  are  known,  then 
knowledge  of  the  sink  attraction  rates  is  equivalent  to  knowledge  of 
the  source  generation  rates.  Thus  the  problem  mentioned  in  section  4.  (.2 
of  knowing  some  p's  and  some  A's  (together  with  the  e^j)  will  give  an 
inconsistent  solution  to  (4.32)  except  in  very  special  circumstances. 
Luckily,  one  such  circumstance  which  arises  naturally  for  an  ATC  network 
is  where  all  p's  and  A's  refer  to  distinct  streams  of  aircraft.  A 
practical  example  of  this  is  discussed  in  section  4.5.6. 


4.4.4  Estimation  of  Transition  Probabilities 


The  estimation  of  transition  probabilities  in  Markov  cliains  has 
been  studied  in  detail  by  Anderson  and  Goodman  (1957)  and  otlier  autliors. 
Their  results  are  applicable  to  the  analysis  of  sector  sequences  witli  onlv 
minor  modifications. 


Suppose  there  are  available  records  of  c independent  sector 

sequences,  each  assumed  to  be  a random  observation  from  a q*’*'  order  Markov 

chain  with  M states.  Let  n^^  be  the  number  of  Limes  an  aircraft  enters 

state  k from  state  i,  n...  be  the  number  of  times  an  aircraft  enters  state 

1 _1  k 
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k from  state  j having  entered  j from  state  1,  and  so  forth.  Equivalently, 


n=l  m=0 


(4.33) 


where  is  an  indicator  function  such  that 


Ik(x)  = 


1,  if  X = k ; 


0,  otherwise. 


Similarly, 


m 

c n 

"jk  = I,  ^jk  (^(">-1)* 
11=1  m=l 


(4.34) 


where  y)  is  a two-dimensional  indicator  function  which  is  equal 

to  one  ,r  x=]  and  y=k  and  zero  otherwise.  Generalizing  in  a similar  manner. 


n=l  m=2 


(4.35) 


Assuming  stationary  transition  probabilities,  Anderson  and  Goodman 
show  that  if  the  Markov  chain  is  first  order,  then  the  nj|.  form  a set  of 
sufficient  statistics  for  the  state  transition  probabilities.  For  a 
second  order  Markov  chain,  the  are  sufficient.  The  results  can  be 

generalized  to  higher-order  chains. 


The  maximum  likelihood  estimates  of  the  transition  probabilities 
depend  on  the  order  of  the  Markov  chain.  For  a second-order  chain,  s^(0) 
and  s^(l)  are  assumed  to  be  nonrandom,  while  Sj^(k),  k _>2,  are  assumed 
to  be  random  variables.  Then  the  maximum  likelihood  estimates  of  (4.24), 
(4.23)  and  (4.22)are  given  by 


ft  = 

‘^ijk  M 


1 " • ■ I 
^ 11 


(4.36) 


’,jk  ’ M M " 

y y n 

^ i i P, 
i=l  H I 


(4.37) 


M M 

1=1  ,1=1 
M M M 


):  y.  y n 

1=1  j = l H=1 


(4.38) 


J 


Notice  that  since  8^,(0)  and  s (1)  are  assumed  to  be  nonrandom,  the 
estimates  of  the  transition  probabilities  involve  summations  over  n 
rather  than  using  n.,  and  n,  directly  (the  results  are  not  equivaleh 

J K K 


It  should  be  remarked  that  a second  order  Markov  chain  of  M states 
can  be  reformulated  as  a first-order  Markov  chain  of  M^  states  by  defining 
a state  in  the  new  cliain  to  be  a transition  between  two  particular 
states  in  the  original  chain.  This  is  equivalent  in  network  terms  to 
considering  the  arcs  of  the  network  to  be  states  rather  than  the  nodes. 

The  increase  in  number  of  states  is  in  practice  much  less  than  a factor 
of  M since  many  possible  arcs  do  not  exist.  Further,  the  one-srep 
transition  matrix  of  the  new  chain  has  many  zeros.  This  change  of 
formulation  is  illustrated  in  detail  in  section  4.5.4. 


4.4.5  Determining  the  Appropriate  Order  of  Harkov  Chain 

Using  (4.37)  and  (4.38),  a likelihood  ratio  test  was  derived 
by  Anderson  and  Goodman  (1957)  for  testing  the  null  hypothesis  that  a 
first-order  Markov  chain  is  sufficient  to  characterize  a set  of  obser- 
vations against  the  alternative  of  a second-order  chain.  The  test  cri- 
terion is 

M 

^ = • • 1 
1,J ,K=1 

Then  if  the  null  hypothesis  is  true,  the  quantity  -2  log^  has  an  asymptotic 
X^-distribution  with  M(M-l)'^  degrees  of  freedom.  .Similar  test.s^  were 
given  for  testing  the  null  hypothesis  that  a chain  is  of  (r-1)  order 
against  the  alternative  that  it  is  of  r^h  order.  If  some  of  the  are 

known  a priori  to  be  zero,  then  the  degrees  of  freedom  should  be  adjusted 
downward  appropriately. 

This  technique  can  be  applied  effectively  to  ATC  sector  sequences 
only  if  the  number  of  sectors  is  small  and  the  number  of  observed  sector 
sequences  quite  large.  For  other  situations,  a graphical  technique  based 
on  information  theory  given  by  Chatfield  (1973)  (which  can  be  related  to 
the  likelihood  ratio  test)  is  all  that  the  data  will  support.  The 
technique  involves  plotting  the  "conditional  uncertainty"  remaining  in  the 
data  as  successively  higher  order  Markov  chains  are  entertained. 

The  concept  of  conditional  uncertainty  is  defined  via  Shannon's 
measure  of  information.  Given  an  event  with  M mutually  exclusive  and 
exhaustive  outcomes,  each  with  prohab  il  i t v Pj^(k=l,  2,  ...,  M)  , tluni  Liie 
average  amount  of  information  provided  by  one  observation  is  given  by 

M 

n = E(-log  p)  = - ) p.log  p . (4.40) 

k=l  " " 


-1-^- 


n. 

ijk 


(4.31) 
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similarly,  the  average  Information  about  pairs  of  successive  events  provided 
by  a single  experiment  is  given  by 

M M 

H(pairs)  = - I I P..  log  p.,  . 
j=l  k=l 

The  difference  between  H(pairs)  and  H is 

H2  = H(pairs)  - H , 

termed  the  average  conditional  uncertainty  about  an  event  given  the 
preceding  event.  The  concepts  are  easily  extended  to  determining  H^j 
the  average  conditional  uncertainty  about  an  event  given  the  two  pre- 
vious events  and  so  forth. 

The  information  theory  technique  for  determining  the  order  of  a 
Markov  chain  involves  plotting  consecutively  the  average  conditional 
uncertainties  about  an  event  given  no  knowledge  of  the  previous  event, 
given  knowledge  only  of  the  previous  event,  knowing  the  two  previous  events, 
and  so  on.  The  resulting  sequence  {log  M,  H,  H2,  H^,  ...}  will  be  monotonically 
decreasing,  and  its  structure  will  help  to  indicate  the  order  of  Markov 
chain  appropriate. 

In  order  that  the  estimated  conditional  uncertainties  will  add 
together  in  the  appropriate  manner,  it  is  necessary  to  first  determine 
how  many  values  will  be  required,  and  then  to  determine  all  estimates  from 
sequences  of  the  same  length.  For  example,  if  the  chain  is  thought  to 
be  of  second  order,  the  first  conditional  uncertainty  expected  to  be 
insignificant  is  that  involving  two  previous  events,  H^.  Wishing  to  plot 
values  up  to  H4,  it  is  therefore  necessary  to  compute  and  estimate 

the  conditional  uncertainties  using  summations  involving  The 

analogy  with  the  maximum  likelihood  estimates  in  section  4, 4,4  is  direct, 
creating  the  first  three  events  in  the  sequence  as  nonrandom.  This  technique 
is  illustrated  in  more  detail  in  section  4.5.5. 


(4.41) 


(4.^2) 


4.4.6  Summary 

The  use  of  Markov  chain  models  to  characterize  sequences  of  sec- 
tors over  which  aircraft  pass  provides  a convenient  formulation  for  traffic 
through  ATC  networks.  It  eliminates  the  need  to  specify  the  relative 
frequency  of  all  possible  sequences  while  maintaining  the  basic  integrity 
of  streams  through  the  network.  The  estimation  of  transition  probabilities, 
together  with  estimates  of  traffic  generation  rates  of  the  sources  to  the 
network,  completely  describe  the  passage  of  traffic  through  the  network  on 
a sector-to-sector  basis.  Sector  transit  times  and  resultant  network  transit 
times  will  be  considered  In  a future  report. 
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A further  note  is  appropriate  about  the  likely  order  of  Markov 
chain  needed  to  describe  sector  sequences.  It  is  obvious  that  since 
aircraft  are  unlikely  to  return  to  sectors  from  which  they  have  just 
departed,  first  order  Markov  chain  models  will  not  be  appropriate.  Tlie 
second-order  model,  however,  would  be  sufficient  to  preserve  straight  line 
streams  through  individual  sectors  and  might  well  prove  adequate  to  des- 
cribe the  important  aspects  of  the  network  flow  patterns.  Higher  order 
models  would  probably  gain  little  in  terms  of  reality  while  sacrificing 
much  in  terms  of  simplicity  and  reliable  estimation.  Only  very  extensive 
data  samples  could  clearly  support  models  of  higher  than  second  order. 

It  thus  remains  to  apply  these  techniques  to  an  actual  ATC  network. 
The  following  section  makes  use  of  a data  sample  from  the  New  York 
ARTCC  which  has  been  analyzed  extensively  in  Hunter  et  al.  (1973)  and 
Hunter  et  al.  (1974).  The  applicability  of  the  methods  to  this  point 
will  be  readily  apparent. 


4 . 5 Analysis  of  the  New  York  Enroute  Network 

4.5.1  Network  Structure 

The  New  York  enroute  sectors  consist  of  32  control  positions: 

9 high  altitude  enroute  (HI)  sectors,  11  low  altitude  enroute  (LC) 
sectors,  and  12  low  altitude  transitional  (LT)  sectors  (see  table  4.1). 

These  sectors  cover  a geographical  area  consisting  of  the  entire  states  of 
New  Jersey  and  Delavrareand  parts  of  New  York,  Pennsylvania,  (Connecticut,  and 
Maryland  (see  figure  4.4).  The  relative  position  of  the  various  sectors  can 
be  seen  in  figures  4.5  and  4.6.  The  HI  sectors  control  traffic  at  and  above 
24,000  ft  while  the  LE  and  LT  sectors  handle  traffic  at  altlttides  of  less 
than  24,000  ft. 

In  describing  the  network  of  enroute  sectors,  there  are  m=32  nodes 
plus  a super-source  and  super-sink.  The  node-node  Incidence  matrix  T),  dis- 
played In  figure  4.7,  indicates  152  arcs  between  sectors  in  the  network. 

With  an  additional  32  arcs  in  (s,N)  and  32  arcs  In  (N,t),  there  are  212  .arcs 
In  total  (of  a possible  1088). 

There  are  12  alriJorts  in  tlie  region  covered  by  the  New  York 
enroute  sectors  wliich  contributed  at  least  five  aircraft  to  the  network  f 1 ow 
observed  during  a 2 -Hour  sample  period.  These  airports  are  listed  in 
table  4.2  with  the  number  of  arrivals  and  departures  observed  tiirough  tlie 
network.  These  airports  are  commonly  referred  to  by  th'‘ir  standard  three- 
letter  codes. 

Figure  4.8  illustrates  the  general  flow  of  aircraft  through  tlu‘  net- 
work. From  the  data  base,  it  was  possible  to  distinguish  four  classes  of 
traffic : 
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Table  4.1 

Identification  of  N.Y.  Enroute  Sectors 


Sector  //  Sector  Sector  Ident. 
Type Code 


451 

LT 

R 5 

ESR 

453 

LT 

R 7 

OOD 

454 

LT 

Rll 

SBJ 

455 

LT 

R12 

ARD 

456 

LT 

R15 

COL 

457 

LT 

R16 

COL 

458 

LT 

R17 

7ID 

459 

LT 

R23 

BDR 

460 

LT 

R24 

RVH 

461 

LT 

R25 

SAX 

462 

LT 

R27 

HPN 

463 

LT 

R28 

HPN 

464 

LE 

R 2 

ABE 

465 

LE 

R 3 

SEC 

466 

LE 

R 4 

HAR 

467 

LE 

R 8 

CYN 

468 

LE 

R 9 

SBY 

469 

LE 

R18 

7 HP 

470 

LE 

R20 

IPT 

471 

LE 

R21 

BGM 

472 

LE 

R26 

STW 

473 

LE 

R29 

HUO 

474 

LE 

R30 

IGN 

475 

HI 

R 1 

PSB 

476 

HI 

lOA 

CYN 

477 

HI 

lOB 

SIE 

478 

HI 

R13 

ARD 

479 

HI 

R14 

HAR 

480 

HI 

R19 

HUO 

481 

HI 

R22 

CMK 

482 

HI 

R35 

7HP 

483 

HI 

R36 

JFK 
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ARTCC  Sector  Control  Boundaries  — High  Altitude  (4/19/68) 


Figure  4.6  New  York  ARTCC  Sector  Control  Boundaries  — Low  Altitude  and  Low  Transitional  (4/19/68) 
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(i)  aircraft  departing  from  airports  within  the  region  covered 
by  the  New  York  enroute  sectors  and  bound  for  destinations 
beyond  the  network  boundaries; 

(ii)  aircraft  arriving  from  outside  the  region  and  bound  for  air- 
ports within  the  region; 

(iii)  aircraft  both  departing  from  and  arriving  at  airports  within 
the  region; 

(iv)  traffic  passing  through  the  enroute  sectors  which  were  neither 
arrivals  nor  departures  from  airports  in  the  region. 


Table  4.2 

Traffic  from  Airports  within  the  New  York  Region 


Airport 

Code 

Departures  ttirough 

Arrivals  through 

enroute  network 

enroutt!  network 

1. 

Newark 

EWR 

53 

40 

2. 

.John  F. 
Kennedy 

JFK 

52 

64 

3. 

LaGuardia 

EGA 

58 

56 

4. 

Philadelphia 

PHI. 

36 

40 

5. 

Atlantic  City 

ACY 

9 

5 

6. 

Wilmington 

TEG 

5 

9 

7. 

Wilkes-Barre 

AVP 

3 

2 

8. 

Binghamton 

BGM 

3 

5 

9. 

Harrisburg 

HAR 

12 

0 

10. 

Allentown 

ABE 

8 

3 

11. 

Elmira 

EEM 

7 

3 

12. 

Westchester 

HPN 

7 

11 

253 

238 

Table  4.3 

Traffic  Breakdown  by  Class 

No.  observed 
in  samp  1 e period 


(1) 

Source  in  N.Y. 

region 

237 

(11) 

Terminus  in  N.Y 

. region 

222 

(ill) 

Source  and  term 

inns  in  N.Y.  region 

16 

(Iv) 

Through  traffic 

90 
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A breakdown  of  the  traffic  observed  during  the  sample  period  is  given 
in  table  4.3.  Of  the  four  classes,  (i)  and  (il)  make  up  the  bulk  of  the 
traffic  observed.  Class  (ill)  appears  to  be  of  little  significance  in 
determining  total  flow. 

Appendix  4-A  describes  the  tabulation  of  traffic  in  the  pnrouto 
network  during  the  sample  period.  Each  aircraft  is  classified  as  a 
departure  from  an  airport  in  the  region,  an  arrival  to  an  airport, 
through  - traffic,  or  unidentifiable  because  its  sector  sequence  was 
Incomplete  at  the  end  of  the  sample  period.  Those  fb  aircraft 
which  both  arrived  at  and  departed  from  airports  in  the  region  are 
classified  as  departures  but  indicated  by  asterisks. 

The  count  of  all  transitions  between  the  sources,  sinks,  and  nodes 
of  tiie  network  is  presented  in  table  4.4.  The  form  of  the  table  is  timt 
of  a one-step  transition  matrix  as  discussed  in  section  4.4.3,  only 
transition  probabilities  have  been  replaced  by  the  counts  n^ . . For  a 
first-order  Markov  chain,  the  table  gives  all  the  information  needed  to 
estimate  transition  probabilities. 


4.5.2  Modeling  Network  Arrivals 

To  apply  the  techniques  of  section  4.3  to  the  N.Y.  enroute  net- 
work, it  is  necessary  to  show  that  the  streams  from  each  of  the  sources 
to  the  network  are  well-represented  by  independent  Poisson  processes  and 
that  the  conditional  probabilities  of  an  arrival  entering  over  any  of  the 
32  arcs  in  (s,N)  given  its  source  remains  constant  over  time. 

For  purposes  of  analysis,  it  is  convenient  to  consider  25  sources 
of  traffic.  Twelve  of  these  sources  represent  departures  from  airports 
in  the  region.  Another  12  represent  the  attraction  of  aircraft  to 
each  of  the  airports  in  the  region.  Tlie  final  source  represents  all 
through -traf  f ic . 

It  is  reasonable  to  assume  that  if  any  of  the  traffic  streams 
from  the  various  sources  deviate  greatly  from  the  structure  of  a Poisson 
process,  it  is  likely  to  be  those  of  the  major  airports  (FVfK,  IFK,  EGA, 

PHL)  where  scheduling  of  the  heavy  traffic  may  reduce  the  variance  of  the 
Interarrival  or  Interdeparture  times.  It  is  important  to  recognize  that 
the  structure  of  the  streams  is  relevant  upon  entrv  to  the  enroute  net- 
work, not  upon  takeoff  and  landing.  For  the  departure  streams,  tliis 
requires  looking  at  the  streams  as  thi>v  leave  the  radar  departure  cc'iit  ro  I 1 1 rs . 
The  arrival  streams  are  harder  to  identify  upon  entry  to  the  network  since 
they  originate  through  many  arcs  in  (s,N),  but  it  is  reasonable  to  ttssume 
that  if  they  appear  to  be  Poisson  upon  leaving  the  network  and  entering 
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Table  4.4  Counts  of  Departures,  Arrivals  and 
Transitions  in  New  York  Network 

'if-T  471  2 47  1 474  4'7S  47A  477  4'^fl  474  unQ  4«1  442  441 

0)02500  7 00000  7 0 0 5 1 

4 0 0 0 0 0 0 0 0 4 0 0 0 0 0 52 

Lli\  0 )0000000000000  5f 

’ ” 7 J 0 0 0 0 0 0 c 0 0 0 0 0 0 

.V.Y  0 J 0 0 0 0 0 0 0 0 0 7 0 '7  0 t- 

M;  070000000  0 00  )70  5 

A V 7 j 0 0 0 0 0 0 0 0 0 0 0 7 C 1 

V’.''  7 7 1 0 0 0 0 0 0 0 0 0 0 0 0 1 


Table  4.4  Counts  of  Departures,  Arrivals  and 
Transitions  in  New  York  Network  (cont.) 


EW9 

JFK 

L 0 A 

PHL 

ACY 

TIG 

AV  ? 

Ecn 

HAR 

ABE 

F.LM 

HPN 

000 

451 
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Table  4.4  Counts  of  Departures,  Arrivals  and 
Transitions  in  New  York  Network  (cont.) 
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the  radar  arrival  controllers,  they  are  likely  to  be  Poisson  at  their 
original  entry  to  the  network. 

Figures  4.9  - 4.16  are  histograms  of  the  empirical  Interarrival 
time  distributions  for  the  departure  and  arrival  streams  at  the  four 
major  airports.  If  the  underlying  streams  are  Poisson,  then  the  inter- 
arrival  and  interdeparture  time  distributions  should  be  exponential. 
Table  4.5  lists  various  statistics  calculated  from  the  empirical  distri- 
butions. In  judging  whether  the  distributions  are  exponential,  one 
should  remember  that  the  mean  and  standard  deviation  of  the  exponential 
distribution  are  equal. 


Since  we  observed  the  traffic  streams  over  a time  period  T of 
2 hours  at  intervals  of  1 second  each,  it  is  best  In  estimating  the 
flow  rate  parameters  to  use  (4.16),  the  number  of  events  observed  divided 
by  the  length  of  the  sample  period.  This  estimator  is  known  to  be  the 
minimum  variance  unbiased  estimator  for  a Poisson  distribution,  satisfy- 
ing the  Cramer-Rao  lower  bound.  Table  4.5  includes  these  estimates  for 
A with  the  appropriate  upper  and  lower  95  per  cent  confidence  limits  eiven 
by 


2A 

9. 


1 

2 


2 

^ 2n,a/2 


(4.43) 


^^u  2 ^ 2(n+l),  l-u/2 


(4.44) 


where  n is  the  number  of  events  observed  in  the  2 -hour  sample  (Johnson 
and  Kotz,  1969).  For  x>100,  the  normal  approximation  to  the  “ distri- 
bution was  used. 

In  testing  whether  the  arrival  and  departure  traffic  streams 
mi  gilt  be  reasonably  represented  by  Poisson  processes,  goodness-of-l  i t 
tests  were  performed  on  the  empirical  interarrival  and  interdeparture 
time  distributions  against  hypotliesi zed  exponentials.  The  first  test 
performed  was  a chi-square  goodness-of-f i t test,  based  on  the  classes 
represented  in  the  empirical  liistograms  (with  an  open  class  to  the 
right  beyond  the  fourth  or  fifth  interx’al  of  60  seconds).  Tlie 
test  statistic  calculated  was 

k (N  -F  )^ 

X?  = );  — (4.45) 

i=l  i 


where  N.  was  the  observed  number  of  Interarrival  times  falling  in  each  of 
tlie  k ctasses  and  F,  was  the  expected  number  given  tlial  the  true  distri- 
bution was  exponential  v;ith  parameter  A.  If  the  null  hyiiothesis  is  correct  , 
tlien  x^  will  liave  a chl-squart'  distribution  with  k-2  degrees  ol  freedom  it 
A is  estimated  from  the  grouped  data  (the  estimates  obtained  from  t lu> 
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actual  data  values  are  "super-efficient"  and  could  give  misleading 
results).  If  0=X  <X  <...<Xj^_j<“  are  the  dividing  points  of  the  k 
classes  (for  thls°exainple,  0.0,  59.5,  119.5,  ...  ),  Johnson  and  Kotz 
(1970)  quote  a result  from  Kuldorf  (1961)  that  the  maximum  likelihood 
estimator  for  X is  given  by  the  unique  solution  to 

k = 0 . (A. 46) 

- I V±-i 

1=2 

A 

Solving  (4.46)  numerically,  the  estimates  X can  be  used  to  calculate 
E and  to  obtain  the  solution  to  (4.45).  Table  4.5  Includes  the  values 
or  from  this  procedure.  Only  one  departure  stream  and  no  arrival 

streams  showed  significant  lack-of-fit  at  the  a=.05  significance  level. 

The  second  goodness-of-f it  test  performed  was  a one-sample 

Kolmogorov-Smlrnov  test.  This  test  compares  an  empirical  distribution 

S (X)  with  a hypothesized  distribution  F(X)  using  the  statistic 
n 

D = sup  {|S  (X)  - F(X)1}.  (4.47) 

n n 

X 

When  F(x)  is  exponential  but  the  parameter  X is  unknown,  then  Monte  Carlo 
studies  have  shown  that  the  standard  tables  for  the  critical  values  of  D 
are  extremely  conservative.  Lilllefors  (1969)  in  a Monte  Carlo  study 
suggested  that  values  obtained  from  standard  tables  such  as  Massey  (1957) 
for  a=.20  were  approximately  a=.05  critical  values  when  X was  estimated 
from  the  data.  For  sample  sizes  over  N=30,  he  notes  that  approximate  5 
per  cent  critical  values  are  given  by  Dp=1.06/»^,  10  per  cent  values  by 
,96//^  and.  1 per  cent  values  by  1.25/v/n. 

Comparing  these  critical  values  to  the  values  of  D calculated  for 
the  arrival  and  departure  streams,  no  arrival  streams  appear  to  be  sig- 
nificant at  a=.l0  or  lower,  but  one  departure  stream  is  significant  at 
a=.05  and  three  at  a=.10. 

The  various  tests  conducted  Indicate  that  the  Poisson  model  for 
aircraft  arriving  at  the  major  airports  is  adequate  but  that  the  model  is 
slightly  questionable  for  the  departure  streams.  The  small  size  of  the 
samples,  however,  makes  it  difficult  to  determine  an  alternative  represent- 
ation for  the  streams.  Further,  the  magnitude  of  the  scheduling  in  the 
departure  streams,  which  reduces  the  variance  of  interdeparture  times,  will 
have  only  a negligible  effect  on  the  streams  over  the  arcs  in  (s,N)  which 
handle  some  fraction  of  the  aircraft  from  several  sources.  These  streams 
are  considered  further  in  section  4.5.3. 

The  assumption  that  the  streams  from  the  various  sources  are 
Independent  is  more  difficult  to  test.  However,  there  is  no  convincing 


k-1  N^(Xj-X^_j) 
i=l^  UX^-X^_^) 
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reason  why  the  streams  should  be  dependent,  at  least  under  the  assumption 
of  constant  traffic  generation  and  flow  rates.  In  fact,  the  real  value 
of  the  model  lies  in  being  able  to  manipulate  the  various  source  generation 
rates  separately.  A more  serious  question  of  Independence  arises  In 
section  4.5.4.  ] 


4.5.3  Modeling  of  Traffic  Entering  Sectors 

The  streams  of  traffic  entering  the  network  over  each  of  the  arcs 
In  (s,  N)  are  under  the  assumptions  of  section  4.3,  Poisson  processes. 
The  maximum  likelihood  estimates  of  the  flow  rates  f(s,  N.)  from  (4.19) 
are  simply  the  total  number  of  entries  to  the  sectors  divided  by  the 
length  of  the  sample  period.  Table  4.6  gives  the  number  of  entries 
observed  in  a 2 -hour  sample  and  the  corresponding  flow  rate  estimates 
for  each  of  the  32  arcs  In  (s,  N). 

It  was  argued  in  section  4.5.2  that  the  effects  of  scheduling 
In  the  departure  streams  would  no  longer  be  apparent  in  the  streams 
entering  the  sectors  of  the  network.  In  order  to  test  this  conjecture, 
the  empirical  distributions  of  Interarrival  times  to  the  32  sectors  of 
the  network  were  tested  against  hypothesized  exponentials  with  the 
parameter  Estimates  in  table  4.6.  Using  Kolmogorov's  statistic  D with 
the  adjusted  critical  values  of  Lllliefors  (1969),  it  was  found  tRat  of 
28  sectors  with  greater  than  five  arrivals  during  the  sample  period,  2 
were  significant  at  the  a=.01  significance  level,  but  the  other  26  were 
within  the  90  per  cent  confidence  limits  for  the  null  hypothesis. 

Since  occasional  significant  values  will  occur  even  if  the  null  hypothe- 
sis is  correct,  the  overall  conjecture  that  the  streams  over  the  arcs  in 
(a,  N)  are  Poisson  appears  to  be  fairly  adequate  and  will  be  adopted. 


4.5.4  Source-Node  Transition  Probabilities 

Once  the  streams  of  traffic  through  the  network  have  been 
identified  and  the  flow  rates  estimated,  all  traffic  can  be  considered 
to  originate  at  the  super-source.  Under  the  assumptions  of  section  4.3, 
the  flow  through  s per  unit  time  will  have  a Poisson  distribution  with 
rate  parameter  equal  to  the  sum  of  the  individual  stream  parameters. 
Conditionally  on  k arrivals  at  s in  a given  interval,  the  distribution 
of  the  k aircraft  over  the  various  sources  will  have  the  multinomial 
distribution  of  (4.12). 

In  determining  where  arrivals  at  the  super-source  actually  enter 
the  network,  it  is  necessary  to  estimate  source-node  transition  proba- 
bilities PuU  defined  in  section  4.3.3.  It  will  be  recalled  that  the 
estimates  of  the  source-node  transition  probabilities  are  given  in  (4.17) 
to  be 

* ik 

Pk|l  “ » 

^ • 
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Table  4.6  Test  Statistics  for  Poisson  Arrivals  at 


Enroute 

Sectors 

Sector 

No.  of  arrivals 

f(s,Nj^) 

Interarrival  times 

X 

s 

451 

28 

14.0 

262.6 

249.6 

0 

453 

32 

16.0 

226.7 

217.6 

0 

454 

20 

10.0 

349.7 

219.4 

0 

455 

3 

1.5 

456 

0 

0.0 

457 

3 

1.5 

458 

32 

16.0 

225.8 

170.8 

0 

459 

28 

14.0 

254.6 

282.5 

0 

460 

18 

9.0 

402.3 

350.4 

0 

461 

33 

16.5 

220.7 

145.9 

0 

462 

28 

14.0 

258.3 

192.1 

0 

463 

12 

6.0 

474.5 

568.9 

0 

464 

11 

5.5 

541.4 

665.9 

0 

465 

12 

6.0 

538.5 

459.9 

0 

466 

36 

18.0 

276.0 

241.6 

0 

467 

34 

17.0 

206.8 

218.9 

0 

468 

27 

13.5 

267.8 

284.2 

0 

469 

10 

5.0 

730.6 

619.1 

0 

470 

20 

10.0 

336.5 

385.2 

0 

471 

22 

11.0 

300.9 

257.7 

0 

472 

27 

13.5 

253.2 

197.0 

0 

473 

8 

4.0 

915.4 

944.4 

0 

474 

21 

10.5 

306.1 

304.7 

0 

475 

31 

15.5 

224.9 

198.6 

0 

476 

24 

12.0 

282.0 

236.7 

0 

477 

13 

6.5 

554.5 

633.9 

0 

478 

22 

11.0 

301.0 

373.2 

0 

479 

26 

13.0 

251.4 

217.2 

0 

480 

17 

8.5 

428.4 

487.2 

0 

481 

28 

14.0 

251.4 

255.6 

0 

482 

8 

4.0 

638.3 

529.6 

0 

483 

13 

6.5 

531.7 

346.9 

0 

** 

stat.  slg,  at 


f(s,N  ) Is  In  aircraft/hour 
_ k 

X,  s given  In  seconds 
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n 


.1344 

•0844,* 

.3039 


.1554 

.0609 

.0963 

•1855,, 

.2580 

.1824 

.2029 

.1490 

.1525 

.1428 

.1332 

.1937 

.1422 

.1495 

.1152 

.2202 

.1662 

.0667 

.0860 

.1772 

.1737 

.1159 

.1262 

.0973 

.0789 

.2185 


a=.01 


0 


where  is  the  number  of  aircraft  generated  at  source  i during  the 
sample  period  and  c.,  is  the  number  of  aircraft  generated  at  s which 
entered  the  network  trough  sector  Nj^.  Table  4.7  lists  the  esHmates 
of  Pj^j.  obtained  from  the  number  of  aircraft  at  each  source  which 
enteria  each  sector  of  the  network. 

The  assumption  that  the  transition  probabilities  are  constant 
over  time  can  be  tested  in  part  by  determining  whether  the  selection  of 
sectors  entered  by  aircraft  in  a source  stream  are  independent  or  whether 
some  attempt  is  made  to  alternate  the  sectors  entered.  This  is  related 
to  the  question  of  independence  among  the  traffic  streams  entering  nodes 
in  the  network  and  requires  a more  detailed  look  at  individual  sources 
of  traffic. 


The  most  likely  areas  in  which  the  streams  would  fail  to  be 
independent  is  where  several  LT  and  T.F,  sectors  are  fed  by  the  same 
departure  controller  at  a major  terminal.  The  following  sectors  fall 
into  such  a category; 


From  EWR 
From  JFK 
From  LGA 
From  PHL 


454,461,472 
458, 463, others 
459,461,462 
451,467 , others • 


If  the  transition  probabilities  are  constant  over  time,  then  the 
probability  that  any  given  exit  will  go  to  a selected  enroute  sector 
would  not  depend  on  the  sector  to  which  the  previous  exit  was  directed. 

For  example,  consider  a departure  controller  feeding  three  enroute  sectors 
labelled  Nj^,  N2,  and  N_.  Now  let  p^  be  the  unconditional  probability  that 
a given  exit  from  the  aei 


eparture  sector  enters  the  j'"'*  enroute  sector,  and 


th 


let  p be  the  conditional  probability  ot  entering  the  J '■ 
given  that  the  last  exit  entered  the  1*-^  sector.  Then  if 
of  enroute  sectors  is  Independent, 


th 


enroute  sector 
the  selection 


Pjj  = Pj  o (4.48) 

The  problem  of  testing  whether  (4.48)  is  true  is  similar  to  that 
encountered  in  the  analysis  of  Markov  chains  in  determining  whether 
successive  events  are  Independent.  Using  the  notation  adopted  by  Chat- 
field  (1973),  let  n be  the  number  of  times  an  exit  to  the  i'-^  sector 
was  followed  by  an  exit  to  the  j^^  sector.  If  as  is  often  done  in  testing 
random  numbers,  the  process  is  made  circular  by  considering  the  last  exit 
to  be  followed  by  the  first,  then  estimates  of  Pj  are  given  by 

T.  n. . 

- i 

p = ^ (4.49) 

^ N 
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Table  4.7  Estimates  of  Source-Node  Transition  Probabilities 


where  the  total  number  of  exits  is 


N = 


I n 


Ij. 


Then  if  successive  exits  are  independent,  the  expected  number  of  times  an 
exit  to  sector  i will  be  followed  by  an  exit  to  sector  j is 


'ij 


- N PjPj 


(4.50) 


Using  (4.49)  the  estimates  of  e^^  are  given  by 


n s n.  . 

1 ‘J  i 

N 


(4.51) 


Anderson  and  Goodman  (1957)  have  shown  that  if  the  assumption  of  independ- 
ence is  correct,  then 


X 


2 = 


(4.52) 


is  asymptotically  distributed  as  with  (c-1)  degrees  of  freedom  where 
c is  the  number  of  enroute  sectors  fed  from  the  single  departure  sector. 


Table  4.8  gives  the  estimates  p,  determined  for  each  of  the 
departure  controllers  at  the  four  major  ■airports,  the  expected  values  from 

(4.50),  estimates  from  (4.51),  and  x^  values  from  (4.52).  A significant 
X^  value  was  obtained  for  LGA,  where  successive  departures  tended  to 
alternate  between  sectors  461  and  462.  No  such  interdependence  was  observed 
among  departure  streams  from  the  other  terminals. 


On  the  basis  of  this  analysis,  there  is  no  compelling  evidence 
that  the  streams  of  traffic  over  the  arcs  in  (s,  N)  are  dependent.  In 
particular,  the  assumption  of  constant  source-node  transition  probahillties 
appears  adequate,  at  least  at  the  level  of  analysis  for  which  the  model 
would  ij'e  appropriate. 


4.5.5  Characterization  of  Sector  Sequences 

The  sector  sequences  discussed  in  appendix  4-A  describe  the  movement 
of  aircraft  through  the  network.  A visual  picture  of  this  traffic  has 
previously  been  presented  in  figures  7.13  - 7.15  of  volume  II,  which  dis- 
played the  number  of  transitions  between  all  pairs  of  high  altitude  sectors, 
pairs  of  low  altitude  sectors,  and  pairs  of  high  and  low  sectors  respectively. 
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Table  4.8  Tests  for  Interdependence  In  Departure  Streams 


EWR 


Observed 


18  10  25  53 


1^ 


Expected 

^^"*^454  461  472 


454 

4 

5 

9 

18 

454 

6.11 

3.40 

8.49 

461 

4 

1 

5 

10 

461 

3.40 

1.89  j 

4.72 

472 

10 

4 

11 

25 

472  ] 

8.49 

4.72 

LI. 79 



1 

JFK 


2nd  Observed 
lst^\^  458  463  othey 


458 

— 

17 

6 

8 

463 

6 

4 

1 

other 

8 

1 

1 

31  11 


10 


LGA 


\^2nd  Observed 
1st  459  461  462 


459 

4 

6 

6 

461 

4 

3 

13 

462 

8 

11 

3 

16  20 


22 


PHL 


Observed 
451  467 


31 

11 

10 

52 


451 

12 

8 

2 

467 

8 

2 

1 

other 

2 

1 

0 

22 

11 

3 

\^">nd  Expected 

1^ 


463 

other 


\'’nd 

lst^\ 


16 

459 

20 

461 

22 

462 

58 

\ '’nti 
lst^\^ 

22 

451 

11 

467 

3 

other 

36 

Expected 
451  467 


458 

463 

other 

L8.48 

6.56 

5.96 

6.56 

2.33 

2.11 

5.96 

2.11 

1.92 

T 

Expected 
459  461 

462 

4.41 

5.52 

6.07 

5.52 

6.90 

7.59 

6.07 

7.59 

8.35 

other 


•v 

13.44 

6.72 

1.83 

6.72 

3.36 

0.92 

1.83 

0.92 

0.25 

* 


Significant 


at  a=.05 


A more  comprehensive  summary  of  observed  transitions  among  the  enroute 
sectors  Is  given  in  table  4.9,  a part  of  the  one-step  transition  matrix 
of  counts  presented  earlier. 

In  order  to  determine  what  order  of  Markov  chain  would  be 
appropriate  to  describe  the  sector  sequences,  the  information  theory 
approach  described  by  Chatfield  (1973)  was  selected. 

Figure  4.17  is  a plot  of  the  estimated  conditional  uncertainties 
in  the  sector  sequences,  using  all  observed  quadruplets  in  the  sample. 
From  the  data  in  appendix  4-A,  = 777  quadruplets  were  tabulated.  The 

following  formulas  were  used  to  calculate  the  conditional  uncertainties: 


H = log  44  (44  states  in  chain) 

o 


H 

H 

H 

H 


1 

2 

3 

4 


-1  Z 

log  N,  - N,  n.  log  n. 

4 4 j 1. • • 1. , 


N 


-1 


Z n.  log  n 


i • • • 


. - . j:.  n 

i • # • i > J i j • 


log 


n..  1 

ij.. 


-1  r 


„ n.  . - Z n,  log  n.  ,, 

i -i  i,  Ijk.  ijk. 


Z n^.  ^ log 


-Ir 


" "ijk. 
i.j.k 


Z n.  „ log  n,  „ 

ijk2.  ijki 

^ ? j » ^ ^ 


It  will  be  noticed  that  sharp  drops  occur  until  H , which  is  typical  of 
second-order  chains.  The  drop  from  H_  to  may  or  may  not  be  significant, 
but  it  does  not  appear  to  be  as  important  as  the  earlier  drops.  No 
rigorous  statistical  tests  such  as  discussed  in  section  4.5  were  performed 
because  of  the  large  number  of  states  in  the  chain  and  the  consequently 
small  number  of  counts  for  all  observed  pairs,  triplets,  and  higher 
sequences  during  the  2-hour  sample  period. 

On  the  basis  of  the  above  analysis,  it  was  decided  that  second 
order  Markov  chains  would  be  sufficient  to  describe  the  patterns  observed 
in  the  sector  sequences.  It  will  be  recalled  from  section  4.4.4  that 
second  order  Markov  chains  on  the  nodes  of  a network  can  be  expressed  by 
equivalent  first  order  chains  on  the  arcs.  Intuitively,  if  the  sector 
which  an  aircraft  enters  depends  on  the  previous  two  sectors  through  which 
it  passed,  an  equivalent  statement  is  that  the  arc  through  which  the  air- 
craft leaves  a particular  sector  depends  upon  the  arc  by  which  it  entered 
that  sector,  but  beyond  that  the  past  history  of  arcs  is  Irrelevant. 


Table  4.10  presents  the  first-order  transition  matrices  on  the 

arcs  of  the  network,  a separate  matrix  for  each  enroute  sector.  The 

maximum  likelihood  estimates  of  p. ..  in  (4.36)  can  all  be  derived  from  the 

table,  sin^e  the  element  in  the  ^ 1*^^  row  and  k*^^  column  of  the  matrix 

for  the  j sector  is  n. , and  the  total  of  the 

Ijk 
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Table  4.9  New  York  Enroute  Network  Transition  Counts 


Table  4.10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network 

ARC-ARC  TRANSITION  MATRIX  FOR  SECTOR  451 


000  453  455  464  466  479  ILG  AVP  TOT 
000  1 0 1 0 2 0 0 1 5 

453  C00200002 

454  15  2 0 0 0 0 0 0 17 

455  300000003 

456  1 0 0 0 0 0 0 0 1 

464  01  00^0001 

466  1 0 0 0 0 0 1 0 2 

PHL  0 0 06  5 10  f'0  21 

TOT  21  3 1 8 7 10  1 1 52 


ARC-ARC  TRANSITION  MATRIX  FOR  SECTOR  453 


000  451  455  456  457  466  467  479  PHL  ACY  ILG  TOT 

000  2 1 5 1 9 0 1 1 4 1 1 26 

451  1 000CC2C0003 

467  003031000007 

468  01000^1000001 

476  007050C0000  12 

PHL  100000000001 
ILG  200000000002 
TOT  6 2 15  1 17  1 3 1 4 1 1 52 


ARC-ARC  TRANSITION  MATRIX  FOR  SECTOR  454 


000  451  455  456  472  EWR  LGA  PHL  TOT 

0 
0 

1 1 

0 

1 18 
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Table  4.10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network  (cont.) 

AhC-ABC  TBANSITIOB  MATRIX  FOB  SECTOR  «55 


000  451  464  469  EMR  LGA  TCT 
000  1 0 1 0 1 0 3 

451  1 0 0 C 0 0 1 

453  0 0 0 0 15  0 15 

454  0 0 0 0 0 1 1 

458  0 0 0 1 0 0 1 

461  0 3 0 0 0 0 3 

464  0 0 0 C 13  2 15 

476  0 0 0 0 1 0 1 

TOT  2 3 1 1 30  3 40 


AEC-ABC  TRANSITICK  MATRIX  FOR  SECTOR  456 


000  451  476  JFK  TOT 


453 

1 

0 

0 

1 

0 

4 54 

0 

0 

4 

0 

4 

457 

0 

0 

0 

i 

1 

458 

0 

0 

0 

2 

2 

461 

0 

2 

3 

0 

c 

476 

0 

0 

C 

10 

10 

478 

0 

0 

0 

16 

16 

TOT 

1 

2 

7 

30 

40 

ARC-ARC  TRANSITION  MATRIX  FOR  SECTOR  457 


456  JFK  LGA  TCT 
OOC  0 0 1 1 

453  0 1 15  16 

464  1 1 0 2 

467  0 2 0 2 

PHL  0 1 0 1 

ACY  0101 
TOT  1 6 16  23 
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Table  4.10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network  (cont.) 

ABC-ABC  TBANSITIOi  MATRIX  FOR  SECTOR  458 


000 

455 

456 

467 

468 

469 

477 

482 

483 

TOT 

000 

0 

0 

0 

0 

0 

1 

Q 

0 

0 

1 

460 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

467 

0 

1 

0 

0 

0 

2 

0 

0 

0 

3 

466 

C 

0 

0 

0 

0 

1 

0 

0 

1 

2 

469 

3 

0 

0 

2 

1 

0 

1 

0 

0 

7 

482 

0 

0 

2 

0 

0 

0 

0 

0 

0 

2 

JFK 

1 

0 

0 

3 

2 

0 

17 

6 

0 

29 

TOT 

4 

1 

2 

c 

4 

4 

18 

6 

1 

45 

ARC-ARC  TRAMSITICK  MATRIX  PCF  SECTOR  459 


OOC 

460 

473 

474 

480 

TCT 

000 

0 

3 

1 

0 

0 

4 

460 

0 

0 

1 

0 

0 

1 

460 

1 

0 

0 

0 

0 

1 

JFK 

1 

1 

0 

1 

0 

3 

L3A 

7 

5 

0 

3 

0 

15 

HPN 

2 

0 

1 

0 

1 

4 

TOT 

n 

9 

3 

4 

1 

28 

ARC-ARC  TRANSITION  HATRIX  FOR  SECTOR  460 


000 

458 

459 

469 

JFK 

TOT 

COO 

2 

1 

2 

0 

1 1 

16 

459 

8 

0 

0 

1 

1 

10 

469 

0 

0 

0 

0 

3 

3 

481 

1 

0 

c 

0 

10 

11 

483 

0 

0 

0 

0 

2 

2 

TOT 

11 

1 

2 

1 

27 

42 
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Table  4.10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network  (cont.) 


iHC-ASC  T8AMSITI0M  HATHIX  FOR  SECTOR  461 


000 

000  0 
459  0 

462  1 

463  0 

464  0 
EHR  0 
LGA  0 
TOT  1 


454  455  456  462  463  471 


1 0 0 

1 0 0 

5 1 1 

4 C 0 

0 0 C 

0 0 0 

7 1 4 

10  2 5 


0 0 1 

0 0 0 

coo 

0 0 0 

0 0 0 

2 3 1 

0 0 0 

2 3 2 


472  473  478  E8R  TOT 

0 0 0 0 2 

0 0 0 0 1 

1 0 0 4 13 

0 0 0 0 4 

0 0 10  1 
3 1 0 0 10 

f 1 0 0 19 

10  2 1 4 50 


ARC-ARC  TRAISITIOM  HATBIX  FOR 
461  463  473  TOT 


461 

0 

0 

2 

2 

473 

2 

0 

0 

2 

481 

8 

0 

0 

8 

JFK 

0 

0 

3 

3 

L6A 

0 

21 

0 

21 

HPM 

3 

0 

0 

3 

TOT 

13 

21 

5 

39 

SECTOR 


462 


ARC-ARC  TRANSITION  MATRIX  FOR 


SECTOR  463 


000 

000  0 

461  0 

462  0 

472  0 

474  1 

481  C 
JFK  0 
TOT  1 


461  471  473  480  TOT 
0 0 10  1 
0 0 4 0 4 

0 4 2 14  20 

0 0 9 0 9 

0 0 C 0 1 

4 0 0 0 4 

0 1 0 10  11 

4 5 16  24  50 
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Table  4.10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network  (cont.) 

ARC-ARC  TRANSITION  HATRIX  ROB  SECTOR  464 


451 

455 

465 

466 
470 
472 
475 
PHL 
BAR 
AdE 
TOT 


000  451  454  455  457  461  465  466  470  472 


0 0 0 0 0 

0 0 0 0 0 

C 0 2 2 2 

0 0 0 0 0 

0 0 10  0 
1 0 C C 0 

0 0 17  11  0 

0 1 0 C 0 

0 0 0 1 0 

0 0 0 0 0 

1 1 20  14  2 


1 0 0 3 4 
0 0 0 1 0 
0 0 0 0 0 
0 0 0 0 5 
0 0 0 0 0 

0 0 3 0 0 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 1 
0 2 10  2 

1 2 4 4 12 


475  460  ABE 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 1 

0 0 0 

0 0 1 

0 0 0 

0 0 0 

1 2 0 

1 2 2 


TOT 

e 

1 

6 

5 

2 

4 

29 

1 

2 

8 

66 


ARC-ARC  TRANSITION  HATRIX  FOR  SECTOR  465 


000 

464 

466 

470 

BAR 

TOT 


000  464  466  470  TOT 

2 4 0 3 9 

2 0 0 0 2 

3 0 0 4 7 

5 0 10  6 

1 0 0 0 1 

13  4 1 7 25 


ARC-ARC  TRANSITION  HATRIX  FOR  SECTOR  466 


000 

451 

464 

465 

PHL 

ILG 

ICT 

000 

0 

0 

5 

3 

6 

1 

15 

4 51 

3 

1 

0 

1 

0 

0 

5 

453 

1 

0 

0 

0 

0 

0 

1 

464 

3 

0 

0 

0 

0 

0 

3 

465 

1 

0 

0 

C 

0 

0 

1 

475 

1 

0 

0 

0 

0 

0 

1 

4 79 

1 

1 

0 

0 

14 

0 

16 

BAR 

5 

1 

1 

1 

0 

0 

8 

TOT 

15 

3 

6 

5 

20 

1 

50 
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Table  4.10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network  (cont.) 


ABC-ABC  TBANSITIOB  HATBIX  FOB  SECTOB  467 


COO 

453 

456 

467 

468 
477 
PHL 
IL3 
TOT 


000  453  457  458  467  468  PHL  ACY  ILG  TOT 


5 12  0 1 
1 0 0 0 0 

1 1 C 0 0 

1 0 0 0 0 

1 2 0 0 0 

1 4 0 C 0 

5 0 0 4 0 

2 0 0 0 0 

17  8 2 4 1 


1 21  3 16 

0 0 2 0 3 

0 2 0 2 6 

0 0 0 0 1 

0 4 0 0 7 

0 0 0 0 5 

3 0 0 0 12 

0 0 0 0 2 

4 8 3 5 52 


ABC-ABC  TBAMSITICS  HATRIX  FOB  SECTOR  468 


000 

000  6 
458  4 

467  4 

477  2 

482  1 

ACY  3 
TOT  20 


458 

0 

0 

0 

0 

0 

1 

1 


467  477  ACY  IlG  TOT 
621  1 16 
0 0 0 0 4 

0 0 0 0 4 

0 0 0 0 2 

C 0 0 0 1 

2 C 0 0 6 

8 2 1 1 33 


ABC-ABC  TRAHSITICB  MATRIX  FOB  SECTOR  469 


000 

000  1 
455  1 

458  3 

460  0 

483  6 

JPK  2 
TOT  13 


458  460  483  TOT 

3 0 0 4 

0 C 0 1 

0 0 14 

10  0 1 

1 2 C 9 

0 0 2 4 

5 2 3 23 
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Table  4.10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network  (cont.) 

&BC-ARC  TR&ISITICi  H&TRIX  FOB  SECTOR  470 


000 

464 

465 

471 

472 

AVP 

ELM 

TOT 

000 

2 

1 

0 

6 

1 

0 

0 

10 

464 

1 

0 

0 

2 

0 

1 

0 

4 

465 

3 

0 

0 

2 

0 

0 

2 

7 

471 

3 

2 

2 

0 

0 

0 

0 

7 

472 

1 

0 

3 

0 

0 

n 

0 

4 

AVP 

1 

0 

1 

0 

1 

0 

0 

3 

ELM 

4 

0 

0 

0 

0 

0 

0 

4 

TOT 

15 

3 

6 

10 

2 

1 

2 

39 

ARC-ABC  IRAHSITIOR  MATRIX  FOR  SECTOR  471 


000 

470 

472 

473 

48C 

BGB 

ELM 

TOT 

OOC 

2 

5 

1 

3 

1 

1 

1 

14 

461 

0 

0 

G 

0 

0 

2 

0 

2 

463 

2 

0 

C 

0 

0 

1 

0 

3 

470 

4 

0 

0 

3 

0 

0 

0 

7 

472 

1 

0 

0 

0 

0 

0 

0 

1 

473 

1 

2 

0 

0 

0 

1 

0 

4 

480 

0 

0 

1 

0 

0 

0 

0 

1 

BGH 

0 

2 

0 

1 

0 

0 

0 

3 

ELM 

1 

0 

1 

0 

0 

0 

0 

2 

TOT 

11 

9 

3 

7 

1 

5 

1 

37 

ABC-ARC  TRAHSITIOH  MATRIX  FOR  SECTOR  472 


000 
COO  2 
454  0 

461  0 

464  0 

470  0 

471  0 

473  0 

475  0 

480  0 

EMR  0 
TOT  2 


462  463  464  470  471 

0 C 0 0 0 

1 0 C 0 0 

0 0 111 
0 4 0 0 0 

0 0 0 0 0 

0 C 0 0 fl 

0 C 1 0 0 

0 3 0 0 0 

0 10  0 0 

0 0 3 3 0 

18  5 4 1 


473  475  480  EBB  TOT 

0 0 0 0 2 

0 0 0 0 1 

0 1 6 0 10 

9 0 0 0 13 

0 0 0 2 2 

0 0 0 3 3 

0 0 0 0 1 

1 0 0 0 4 

0 0 0 0 1 

0 3 14  0 23 

10  4 20  5 60 
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Table  4.10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network  (cont.) 

ARC-lfiC  TBAMSITIOR  IliTBIX  FOB  SECTOB  473 


000  462  471  472  473  474  481  HPK  TOT 
000  101100104 

459  1 0 1 1 0 0 0 0 3 

461  200000C02 

462  1 0 1 0 0 3 0 0 5 

463  40000435  16 

471  1 0 0 0 0 1 5 0 7 

472  2 0 0 0 0 6 1 1 10 

473  0 0 0 0 0 1 0 0 1 

474  022200006 

480  COC010203 

TOT  12  2 5 4 1 15  12  6 57 


ABC-ABC  TBAHSITICI  HATBIX  FOP  SECTOR  474 


000  463  473  480  481  HPB  TCT 
000  2 0 8 2 5 3 20 

459  4 0 0 0 0 0 4 

473  14  1 0 0 0 0 15 

481  2 0 0 0 0 0 2 

TOT  22  1 8 2 5 3 41 


AHC-ARC  TRANSITION  MATRIX  FOB  SECTOB  475 


000  464  466  472  475  478  479  480  TOT 
000  0 25  0 2 0 0 0 1 28 

464  100C00001 

472  1C0004005 

475  1 0 0 0 0 0 0 0 1 

478  000100012 

479  000100023 

480  7 1 2 0 1 7 1 C 19 

TOT  10  26  ? 4 1 11  1 4 59 
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Table  A. 10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network  (cont.) 

ABC-AFC  TRANSITION  HATP.IX  FOR  SECTOR  476 


000 

453  455 

456 

477 

482 

TOT 

000 

0 

12  1 

4 

4 

1 

22 

456 

3 

0 0 

0 

4 

0 

7 

477 

C 

c c 

6 

0 

0 

6 

TOT 

3 

12  1 

10 

8 

1 

35 

ARC- 

ARC 

TRANSITION 

NATRIX 

FOR  SECTOR 

000 

467 

468 

476 

479  48 

3 

TOT 

OOC 

2 

3 

0 

6 

0 

0 

11 

458 

6 

0 

1 

0 

11 

0 

18 

467 

1 

0 

C 

C 

0 

1 

2 

468 

C 

0 

0 

0 

0 

2 

2 

476 

2 

0 

0 

0 

0 

4 

6 

483 

C 

1 

3 

0 

1 

0 

5 

TOT 

11 

4 

4 

6 

12 

7 

44 

ABC- 

ABC 

TRANSITION  MATRIX  FOR  SECTOR 

000 

456 

475 

TOT 

OOC 

1 

17 

1 

19 

461 

1 

0 

0 

1 

475 

1 1 

0 

c 

11 

479 

0 

n 

1 

1 

TOT 

13 

17 

2 

32 
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Table  4.10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network  (cont.) 


ARC'&RC  TBiHSITIOB  HiTRIX  ?0B  SECTOR  479 


000 

466 

475 

478 

479 

TOT 

000 

3 

14 

3 

1 

1 

22 

451 

8 

0 

0 

0 

0 

8 

475 

1 

0 

C 

0 

0 

1 

477 

11 

0 

0 

0 

0 

11 

479 

1 

0 

0 

0 

0 

1 

TOT 

24 

14 

3 

1 

1 

43 

ARC-ARC  THAHSITIC*  MATRIX  FOB  SECTOR  480 


000 

459 

471 

472 

473 

475 

481 

TOT 

000 

3 

0 

1 

0 

2 

5 

4 

15 

459 

0 

0 

0 

0 

0 

1 

0 

1 

463 

17 

0 

0 

0 

0 

2 

0 

19 

464 

2 

0 

0 

0 

0 

0 

0 

2 

471 

0 

0 

0 

0 

0 

1 

0 

1 

472 

8 

0 

0 

0 

0 

9 

0 

17 

474 

2 

0 

0 

0 

0 

0 

0 

2 

475 

0 

1 

0 

1 

1 

0 

1 

4 

481 

0 

0 

0 

0 

0 

2 

0 

2 

TOT 

32 

1 

1 

1 

3 

20 

5 

6 3 

ARC-ARC  TRARSITIOH  MATRIX  FOR  SECTOR  481 


000 

473 

474 
480 
TOT 


000  460  462 

1 3 7 

14  0 

0 1 0 

1 3 0 

3 117 


3 474  480  483 

4 0 2 1 

0 10  0 

0 10  0 

0 0 0 1 

4 2 2 2 


LG  A 

HPM 

TOT 

8 

0 

26 

7 

2 

15 

3 

0 

5 

0 

0 

5 

18 

2 

51 
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Table  4.10  Arc-arc  Transition  Matrices  with 
Counts  in  New  York  Enroute  Network  (cont.) 

ARC-APC  TRANSITION  BATRIX  FOR  SECTOR  482 


000  458  483  TOT 
COC  1 1 6 e 

4 58  5 0 0 5 

476  0 0 1 1 

TOT  6 1 7 14 


ARC-ARC  TRANSITION  BATRIX  FOR  SECTOR  483 


OCC 

458 

469 

477 

481 

482 
TOT 


COO  460  469  477  TOT 

0 2 4 6 12 

1 0 0 0 1 

3 0 C 0 3 

2 0 6 0 8 

0 10  0 1 
1 3 C C 4 

7 6 10  6 29 
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. th 

1 row  is 


Then 


E 


^’ijk 


E "ijt 
il=l 


where  c Is  the  number  of  columns  in  the 


, th 


matrix. 


Tables  4.6,  4.7,  and  4.10  thus  completely  describe  the  sequence 
of  sectors  through  which  aircraft  move.  Given  X for  all  sources, 
for  the  probability  of  aircraft  generated  at  each  source  entering  ' 

the  network  over  arc  (s,Nj^),  and  the  probabilities  p. for  moving  air- 
craft through  the  network,  the  traffic  is  comnictel;-  Specified. 


4.5.6  Relating  Sink  Attraction  Rates  to  Source  Generation  Rates 

The  problem  described  earlier  for  determining  source  generation 
rates  from  sink  attractions  can  now  be  solved  under  the  assumption  that 
class  (ill)  traffic  is  insignificant.  (Class  (iil)  is  the  traffic  that 
both  departs  and  arrives  at  airports  in  the  region.)  In  the  data  sample, 
this  class  represented  less  than  3 per  cent  of  the  traffic  through  the  network. 


If  it  is  assumed  that  class  (iil)  traffic  is  not  significant, 
then  all  traffic  attracted  to  the  airports  in  the  region  can  be  con- 
sidered to  enter  the  sector  through  the  enroute  source.  Thus 


X 


e 


\ + 


m 

t 


(4.53) 


where  X is  the  rate  of  through-traffic  passing  through  the  network. 
If  X ana  the  ij's  are  used  to  determine  X , equation  (4.32)  will  be 
consistent . 


The  final  problem  involves  determining  the  conditional  proba- 
bility Pj^i  for  entry  to  the  network  from  the  enroute  source.  The 

' Q 

solution  proposed  is  to  specify  the  sink  attraction  rates  p.  and  then 
backsolve  using  the  transitional  probabilities  p to  dete'imine  the 
arrival  components  of  P^  i ^ • The  tiirough-components  would  be  tlie  weighted 


e 

difference  between  the  arrival  components  and  the  estimates  of 
determined  from  all  traffic  at  the  enroute  source. 


s 

e 
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A. 5. 7 Summary 


A 


The  formulation  presented  in  this  chapter  could  be  applied  to  a 

I general  ATC  network  to  determine  the  flow  patterns  over  ti.e  network 

given  specified  rates  of  traffic  arriving  and  departing  from  sources 
I and  sinks  to  t^>e  network.  An  extensive  data  base  Is  npaded  If  all  transi- 

I tion  probabilities  are  to  be  estimated,  but  for  specified  probabilities 

the  techniques  are  easily  applied. 

; It  should  be  recognized  that  the  aspects  of  the  network  studied 

! in  this  report  only  give  a microscopic  description  of  traffic  flow. 

Important  variables  such  as  transit  times  and  communication  measures 
have  yet  to  be  studied.  The  attempt  here  has  been  to  formulate  the  ATC 
system  in  network  terms  and  to  describe  the  general  patterns  of  flow  In 
such  a way  that  a simulation  model  could  be  easily  constructed. 

The  New  York  enroute  network  example  was  perhaps  a more  compli- 
cated network  than  would  usually  be  studied.  For  oceanic  systems,  the 
number  of  sectors  will  be  much  smaller  and  the  traffic  patterns  much 
simpler.  The  same  sort  of  analysis  as  was  applied  to  the  N.Y.  enroute 
netv.'ork  in  this  report  could  thus  be  applied  with  much  less  effort  to 
these  simpler  networks,  yet  the  possible  conclusions  would  be  at  least 
as  powerful. 

Finally,  it  is  felt  that  while  the  network  model  at  this  level 
of  sophistication  could  produce  Interesting  conclusions  about  the  locations 
of  heavy  traffic  flow  in  networks,  stronger  statements  will  be  possible 
after  transit  times  and  communication  requirements  are  introduced. 
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APPENDIX  4-A 

TABULATION  OF  ENROUTE  SECTOR  TRAFFIC 

Technical  report  PU-42  includes  a list  of  all  traffic  observed 
through  the  N.Y.  enroute  sectors  during  a 2 -hour  sample  period.  The 
list,  which  provided  the  raw  data  for  use  in  this  chapter,  includes 
the  following  information: 


•k 

1  Classification  of  aircraft  type 


2  Aircraft  identification  code 


3  Entry  and  exit  movement  patterns 


4  Sample  entry  and  exit  times  in  seconds  (0-7200) 


5  Traffic  class 

a.  ARV  - arrival  from  an  airport  in  the  region 

b.  DEP  - departure  from  an  airport  in  the  region 

c.  (blank)  - through-traffic 

d.  -----  _ unidentifiable  since  sector 

sequence  is  cut  by  end  of  sample  period 


6  Sector  sequence  in  order  of  entry  times. 


For  description,  see  Hunter  et  al.(1973) 


The  following  sectors  were  held  to  indicate  arrivals  and 
departures  from  airports  in  the  region: 


EWR-ARV 

546,547,548 

EWR-DEP 

534 

JFK-ARV 

549,550,551 

JFK-DEP 

535 

LGA-ARV 

552,553,554 

LGA-DEP 

536 

PHL-ARV 

555,556 

PHL-DEP 

537,538 

ACY-ARV 

545 

ACY-DEP 

533 

I LG- ARV 

511,543 

ILG-DEP 

511,543 

AVP-ARV 

525,539 

AVP-DEP 

525,539 

BGM-ARV 

526,540 

BGM-DEP 

526,540 

HAR-ARV 

509,541,542 

HAR-DEP 

509,541,542 

ABE-ARV 

524,529 

ABE-DEP 

524,529 

ELM-ARV 

527,530 

ELM-DEP 

527,530 

HPN-ARV 

510,531 

HPN-DEP 

510,531 
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APPENDIX  4-B 


LIST  OF  NOTATIONS 


Notation 


m 

G 

N 

A 


IJ 


t 
D 
d, 

F 

(X.Y) 

f(X,Y) 


s . 
1 


m 

s 


A . 
1 


u . 

I 


X 

s 


//  of  sectors  (nodes)  in  a network  G 

sector  k (k=l,2,...,  m) 

network 

set  of  nodes  in  G 
set  of  arcs  in  G 
super-source 
super-sink 

node-node  incidence  matrix 

element  in  i*"^  row  and  j*"*'  column  of  D 

traffic  flow  rate  from  s to  t 

set  of  arcs  from  nodes  in  set  X to  nodes  in 
set  Y 

traffic  flow  rate  over  arcs  in  (X,Y) 

source  i (i  = 1,2,  ...,  m^) 

//  of  sources  with  arcs  to  s 

sink  j (j  = 1,  2,  . . . , m^) 

of  sinks  with  arcs  from  t 

traffic  generation  rate  at  source  i 

traffic  attraction  rate  at  sink  j 

prob,  that  aircraft  enters  network  at  sector  k 

prob.  that  aircraft  generated  at  source  1 
enters  network  at  sector  k 

traffic  flow  rate  through  s 
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(s^(k). 


P 


(m^x  1)  vector  of  source  generation  rates 
(m^xm)  matrix  of  source-node  entry  probabilities 


Ik 


k=0,l,2,  ...}  ' 
m 

n 

n(q) 


ij 


■ij 


"jk’  "ijk 


^ijk 

'^jk’  ’’ijk 
X 


k s 


= source-node  entry  count  matrix 

= element  of  C equal  to  the  number  of  aircraft 
which  were  generated  at  source  1 and  entered 
the  network  over  arc  (s,N|^) 

= number  of  aircraft  generated  at  source  i 
during  the  sample  period 

= number  of  time  units  In  the  sample  period 

th 

= sequence  of  sectors  traversed  by  n aircraft 

th 

= number  of  network  sectors  traversed  by  n air- 
craft 

= q-step  transition  matrix 

= prob.  of  going  from  sector  i to  sector  i In 
q-steps 

= probability  of  an  arrival  from  source  i ending 
up  at  sink  j 

= counts  of  observed  transitions  between  sectors  In 
a sample  of  sector  sequences 

(x,y,z)  = indicator  functions  for  sector  sequences 

= state  transition  probabilities 

= flow  rate  of  through-traffic  in  C 

= arrival  rate  of  enroute  traffic 

= conditional  probability  of  entering  the  network 
over  arc  (s,Nj^)  given  generation  at  the  enroute 
source. 
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CHAPTER  5 


f 


I 


I 


t, 


I 


MISCELLANEOUS  RELATED  STUDIES 


5.1  Introduction 

As  has  been  explained  in  previous  chapters,  in  the  ATC  simulation 
model,  seven  input  parameters  plus  a hlstorica],  histogram  are  used  to 
characterize  the  entire  input  system.  Questions  have  been  asked  as  to  how 
much  each  of  these  variables  contributes  to  the  expected  utilization  of  the 
communications  channel.  To  answer  this  question  a mathematical  analysis  has 
been  carried  out  and  reported  in  section  5.2. 

In  practical  applications  of  the  simulation  model,  particularly  for 
measuring  the  operational  stability  of  the  ATC  system,  indices  which  are 
easy  to  understand  and  convenient  for  use  are  of  great  help.  A proposal  on 
the  compilation  of  a triplet  of  three  operational  Indices,  each  representing 
a major  response,  is  given  in  section  5.3. 

Continuing  effort  has  been  devoted  to  the  exploration  of  the 
coupling  of  responses  observed  in  adjacent  sectors.  Transfer  function 
models  which  link  the  aircraft  loading  series  at  two  functionally  related 
sectors  have  been  constructed  using  Box-Jenkins  modeling  techniques,  and 
reported  in  section  5.4. 

A statistical  investigation  on  four  sets  of  S -hour  data 
collected  from  New  York  center  (Sectors  454LT,  515LC,  536UP  and  554AR)  was 
completed  and  results  reported  in  section  5.5.  This  study  will  be  useful 
for  constructing  a simulation  model  with  unstable  input  flows,  observed  over 
long  periods  of  time. 

To  pursue  a statistical  model  satisfactorily  describing  the  queuing 
time  processes  generated  from  a simulation  model,  a transformation  tecliniqiie 
was  developed  to  make  the  series  easier  for  parametric  time  series  analysis. 

The  resulting  statistical  model  is  useful  for  queuing  time  forecasting  in  a 
nonhomogeneous  arrival  system  such  as  ATC  communications.  The  potential 
usefulness  of  this  study  to  ATC  network  simulations  is  elucidated  and 
reported  upon  in  section  5.b. 

5.2  Sensitivity  Analysis  of  the  Expected  ATC  Communications  Channel 

Utilization 

In  the  simulation  model  which  has  been  developed  for  ATC  systems 
to  date,  seven  initiating  parameters  plus  one  historical  histogram  are  used 
to  characterize  the  entire  input  system  to  the  simulation  framework. 

(Questions  have  been  asked  about  how  Important,  relatively,  each  of  these 
input  parameters  is  to  the  overall  expected  channel  utilization.  To  answer 
these  questions  a sensitivity  analysis  is  developed  and  presented  in  tins 
sect  ion. 

1 


5-1 


I 


5.2.1  Derivation  of  the  Sensitivity  Equations 

As  earlier  demonstrated  in  the  report,  volume  III,  parameters  of 
major  importance  in  the  simulation  model  are: 

0 : number  of  aircraft  arrivals  per  hour; 

a, A : parameters  of  a gamma  distribution  fitted  to  the 

transmission  lengths; 

p,  k : the  two  parameters  in  a negative  binomial  distribu- 

tion fitted  to  the  number  of  CT's  per  aircraft 
(shifting  the  origin  by  two  for  the  data  of  Enroute 
sectors,  and  by  one  for  the  remaining); 

a^.a^^  : regression  coefficients  for  generating  intercom- 

munication gap  lengths,  depending  on  #CT/AC; 

T : the  expected  number  of  transmissions  per  communi- 

cation transaction. 

Important  properties  and  relationships  among  these  parameters  are  briefly 
reviewed  and  displayed  in  the  following  mathematical  equations: 

(i)  K(‘i'R)  = expected  transmission  length  (seconds) 

= a • A = b -b^0  = 3.70  - 0.01586  (for  all  except 

° TX  Sectors) 

= 2.97  - 0.00420  (for  LC  Sectors) 

(ii)  E(#TR/CT)  = expected  number  of  transmissions  per  trans- 
action = T 

(iii)  E(#CT/AC)  = expected  number  of  transactions  per  aircraft 

= k(—  - l)+n  v/here  m=2,  for  Enroute  sectors,  and 
P 

m=l,  otherwise. 


If  we  now  assume  the  three  variables  TR,  //TR/CT  and  //CT/AC  are  mutually 
independent  and  that  the  values  of  the  parameters  a^  and  a^  are  moderate, 
it  can  be  shown  that: 
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E(CU)  = expected  channel  utilization.  In  seconds  per  hour 
= 0 • E(TR)*E(#TR/CT)E(//CT/AC) 

= 0 • (b^-bj^O) 'T*  [k(-^  -l)+m] 

= b • 0(l-b,0/b  )*T*fk(^  -l)+in] 
o 1 o P 


Let  y = E(CU)  and  x be  any  variable  on  the  right  side  of  Equation  (5.1). 
The  elasticity  n , that  is,  the  percentage  change  in  y resulting  from  a 
one  percent  change  in  x,  is  then  defined  as 


h . = 


8 log  y _ ^ 


31og  X 


dy 


= ly  / 

~ y ' X 


ih'.  , dx 
3x  dlog  X 


From  Equation  (5.1),  taking  logarithms  on  both  sides,  we  find  that 

^1  1 

log  y = log  b + log  0 + log  (1  - —0)  + log  T + log  [k(-  -l)+m] 
o b p 

o 


To  reduce  notation,  let  w = [k(p  - 1)  +m],  then  the  elasticities  of  the 
expected  channel  utilization  with  respect  to  t,  w and  0 are 


and 


, = 14°^  = 1 

T 3 log  T 


„ = 44°^  = 1 

w 3 log  w 


= 1 - 


0 


(b^/b^)-0 


(^o^^l)-20 

(b^/bj^)-0 


(5.2) 

(5.3) 

(5.4) 


Using  Equation  (5.2)  and  (5.3)  we  note  that  a 1-  percent  change  in  either 
E(//TR/CT)=  T or  E(#CT/AC)  = w,  but  not  both,  will  generate  a 1-per  cent 
change  in  the  value  of  E(CU),  given  all  the  other  elements  in  the  equation 
being  held  constant.  On  the  other  hand.  Equation  (5.4)  tells  us  that  the 
relative  change  in  E(CU)  is  always  smaller  than  unity,  when  the  traffic 
density,  9,  is  changed.  In  other  words,  when  the  traffic  d-nsity  is  reduced 
(increased)  by  50  per  cent,  th<;  channel  utilization  will  drop  (rise)  by  less 
than  50  per  cent. 


Furthermore,  to  pursue  the  elasticity  measure  p one  more  step,  we 

find 

n,  = k •[k(-  - l)+m]"^  (i  - 1)  (3.5) 

k p p 

Ip  = - 1)+'"]  ^ (k/p)  (3.0'/ 
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Though  the  physical  meanings  of  rij^  and  n are  difficult  to  express,  a 
numerical  illustration  will  be  illuminating.  Using  the  N.Y,  HI  sector 
function  as  an  example  we  calculate  the  elasticities  n , n,  and  n for 
values  of  0,  k and  p,  close  to  those  observed  in  the  historical  dSta.  The 
resulting  curves  are  displayed  in  figures  5. 1-5. 3.  (In  figures  5.2  and  5.3 
since  the  estimates  of  (k,p)  for  the  HI  sector  function  are  (k=4.03,  p=0.58), 
we  computed  at  p=0.58  and  computed  n at  k=4.03).  Computations  for  other 
sector  functions  can  be  similarly  performed. 


5.2,2  Concluding  Remarks  on  the  Sensitivity  Analysis 

In  the  discussion  presented  above  we  have  ignored  the  influence  of 
the  driving  variables  upon  queuing.  It  is  naive  to  judge  the  influence  of 
a variable  solely  based  on  the  expected  channel  utilization.  A sensitivity 
analysis  for  the  queuing  response  requires  large  scale  simulation  experi- 
ments and  is  postponed  for  a future  study. 


5.3  Construction  of  Operational  Indicies  for  ATC  Communication  Per- 

formance 


A characterization  of  various  major  responses  in  ATC  communication 
channels  has  been  presented  in  earlier  reports.  In  volume  III  sophisti- 
cated statistical  techniques  are  used:  to  construct  a simulation  model,  to 
validate  the  simulated  responses,  and  to  detect  operational  instability 
under  heavy  traffic  situations.  However,  much  of  this  material  is  mathe- 
matical, containing  expressions  usually  not  of  great  value  to  those  with 
the  responsibility  of  actually  operating  ATC  processes.  In  this  study  we 
explore  some  equivalent,  but  simple,  indicators  of  the  operational  sound- 
ness (stability)  of  ATC  communications  performance.  The  objective  is  to 
stimulate  interest  in  experimenting  with  the  ATC  communication  simulation 
model  by  ATC  personnel. 

As  discussed  in  chapter  3 of  this  volume,  communication  capacities 
of  ATC  communications  functions  were  measured  by  the  average  responses 
generated  from  a series  of  simulation  exercises,  along  with  a nonstatlon- 
arlty  detector  defined  as  P (C>1),  imposed  upon  the  aircraft  loading  time 
series.  To  express  the  operational  stability  or  soundness  of  an  ATC 
channel  in  a simple  and  Informative  manner,  ideas  that  produced  the 
measurements  of  many  physical  phenomena  such  as  temperature,  humidity, 
etc.,  can  be  of  great  use.  For  Instance,  the  boiling  temperature  of  water 
under  the  normal  pressure  of  gravity  was  set  equal  to  100°  in  the 
Celsius  scale.  In  this  case,  the  number  "100"  is  used  to  Indicate  a crit- 
ical condition  of  water  transformations.  Similarly,  the  operation  of  an 
ATC  system  can  be  expressed  and  quantified  in  the  way  temperature  is 
represented  by  a thermometer.  \Jsing  such  a system  the  communication  per- 
formance can  be  read  and  understood  at  the  instant  a measurement  is  made. 
However,  the  determination  of  a critical  point,  such  as  the  100° 
design  for  a Celsius  thermometer,  for  an  ATC  system  depends  on  the 
system  conditions,  the  ability  of  a controller,  etc.  More  information  and 
understanding  is  needed  to  formulate  these  measurements.  Some  behavior 
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inherent  to  the  ATC  systems  is  first  described  in  the  first  subsection  witli 
some  analogies  from  physical  sciences.  Numerical  illustrations  are  then 
provided  in  the  subsequent  subsection  for  examples  previously  studied. 


5.3.1  Some  Basic  Considerations 

Physics  tells  us  that  when  the  concentration  of  a solvent  in  the 
solution  tends  closer  to  its  saturation  point,  an  exponential  increase  in 
time  or  pressure,  is  required  to  achieve  each  incremental  gain  in  con- 
centration. Further,  many  physical  phenomena  are  multivariate  in  aspect, 
requiring  the  simultaneous  monitoring  of  several  responses.  Similarly,  in 
a communications  channel,  when  the  percentage  channel  utilization  approaches 
its  celling,  say  85  per  cent,  the  demands  made  upon  a controller  are  nor  pro- 
portionally but.  Instead,  exponentially  related  to  incremental  changes  in 
channel  utilization.  And  further,  variables  in  addition  to  channel 
utilization  are  necessary  to  characterize  sector  communications  loading. 

Following  this  reasoning,  we  have  selected  exponential  models  for 
the  construction  of  three  separate  operation  indices  to  represent  communi- 
cations loading.  An  index  is  compiled  for  each  of  three  major  responses: 
aircraft  loading,  channel  utilization  and  comir.unications  queuing.  For  tlie 
reasons  just  stated  in  an  analogy  of  Celsius  thermometer  the  number  100  is 
selected  to  be  used  as  a critical  value  for  each  index.  Thus,  when  a 
response  yields  an  index  close  to  100,  that  response  is  about  to  become 
unstable  and  unsatisfactory  operation  of  the  channel  destined  to  ensue. 


5.3.2  Operational  Indices  for  ATC  Communications 

The  simulation  model  for  ATC  communications  within  a sector, 
described  in  volume  III,  can  be  exercised  by  varying  any  of  the  following 
parameters  previously  defined,  either  singly  or  in  combination:  6,  p,  k, 

a,  X,  So,  a^,  and  the  distribution  of  #TR/CT.  Further,  by  using  the  mes- 
sage element  dictionaries,  additional  modifications  in  the  number  and  lengt';.; 
of  transmissions  are  possible.  For  any  particular  setting  of  these  input 
parameters  three  time  series  responses  are  generated  in  the  simulation 
representing  the  burden  of  air  traffic  and  communications  in  the  sector: 

n^,  : number  of  aircraft  present  per  minute, 

c^  : channel  utllizatioh  per  minute, 

q^  : number  of  aircraft  in  communications  queue  per  minute. 

A series  equivalent  to  (tlj.)  Is  : 

Q : communications  queuing  times  per  CT,  wliere  i means  "for 

the  i*^h  Qx  attempted” 

A second  order  autoregressive  time  series  model  was  found  adequate  for 
describing  the  series  {n  ).  Furthermore,  overall  averages  of  c'nannel 
utilization  or  queuing  time  can  be  calculnt'nl  from  entire  si'ries  of  (c  } 
or  {Q^}.  5_7 


In  volume  III,  chapter  4,  the  aircraft  arrival  rate  S was  experi- 
mentally varied  to  determine  the  effect  upon  the  three  responses.  Much 
valuable  Information  was  gained  from  a close  study  of  the  three  time  series 
responses.  The  Information  can  now  be  encapsulated  by  the  three  indices 


Let  the  multivariate  aspects  of  the  sector  communications  loading 

be  written  as  a triplet  of  three  indices  I , I,  , I 

a D c 

where  I : an  index  which  indicates  che  stability  of  the  aircraft 

^ loading  time  series; 

I^  : an  index  determined  from  the  recorded  channel  utilization, 
as  a fraction; 

I : an  index  on  communication  queuing  times  generated  from 
the  simulation  model. 


Useful  functional  forms  which  express  these  three  indices  are: 


I 

a 


I 

c 


fj^(x)  = 100  e 


(1.645  - x) 


.0.20 

1-y 


1) 


f^(y)  = 100  e 


f^(z)  = 100 


(z/50  -1) 


(5.7) 


Although  the  constants  in  these  index  models  are  somewhat  arbitrary,  tlic 
results  of  indexing  applied  to  some  data  generate^  by  simulations,  appears 
consistent  with  our  a priori  judgment.  The  quantities  x,v,  and  z are 

X = where  f.=  l-  <(>2  ^nd  where  4>,  and  are  the 

autoregressive  parameters  of  the  aircraft  loadin. 
time  series  generated  from  the  simulation  model 
and  a ^ is  the  associated  standard  deviation  ; 

y = overall  channel  utilization  with  the  value  be- 

tween zero  and  one; 

7-  = average  queuing  time  (seconds)  per  communications 

transaction. 


Plots  of  the  indices  I^,  Ij^  and  If.  as  functions  of  x,  y,  and  z are 
displayed  in  figures  5.4  - 5.6.  Numerical  illustrations  using  various 
examples  demonstrated  in  volume  III  are  provided  in  tables  5.1  and  5.2. 
More  specifically,  table  5.1  presents  the  movements  of  the  indices  at 
different  traffic  densities  while  table  5.2  provides  a set  of  indices 
for  various  functions  exercised  at  the  historical  traffic  densities. 
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Figure  5.6  Operational  Index  1^ 
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Table  5.1  ATCC  Operational  Indices  for  N.  Y.  Sector  475  HI 
at  Various  Experimental  Traffic  Densities 


Traffic 

density 

X 

y 

z 

(sec.) 

('a’  ^c) 

50% 

4.4 

0.3290 

3.861 

( 6,  50,  40) 

100% 

4.0 

0.5631 

10.428 

( 9,  58,  45) 

125% 

3.6 

0.6976 

20.526 

( 14,  71,  55) 

150% 

1.6 

0.8353 

56.818 

(105,124,115) 

Traffic  density  100%  = 33  aircraft  arrivals/hour. 
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Table  5.2  ATCC  Operational  Indices  for  Eight  N.  Y.  Sectors  at 
Historical  Traffic  Densities 


Sector  ID 

X 

y 

z 

(seconds) 

( ^a’ 

I ) 
c 

453  LT 

2.100 

0.53 

12.530 

( 63, 

56, 

47) 

475  HI 

3.965 

0.56 

10.430 

( 9, 

58, 

45) 

504  GN 

2.700 

0.47 

15.405 

( 35, 

54, 

50) 

510  LC 

2.415 

0.47 

7.215 

( ^-6, 

54, 

42) 

524  LG 

1.895 

0.34 

9.021 

( 78, 

50, 

44) 

534  DP 

2.515 

0.47 

6.782 

( 42, 

54, 

42) 

543  AD 

2.475 

0.46 

16.691 

( 44, 

53, 

51) 

553  AR 

1.940 

0.59 

12.680 

( 74, 

60, 

47) 
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'5.3.3  Concluding  Remarks  on  Operational  Indices 


1 

I 


! 


;) 


;i 

'i 


A composite  index  can, of  course, be  computed  by  simply  averagin' 
the  three  sub-indices,  or  by  using  a weighting  function  for  averaging. 
Future  experience  should  develop  the  indices  discussed  in  this  section 
into  a more  mature  form. 

It  is  hoped  that  experimentation  with  the  simulation  model  will 
be  stimulated  through  the  use  of  these  indices  of  ATC  communications 
performance . 


5.4  A Transfer  Function  Model  for  Aircraft  Flows  at  Adjacent  ATC 

Sectors 

In  extending  the  ATC  sector  communications  simulation  model  to 
describe  a general  ATC  network,  it  is  important  to  characterize  the  flow 
of  traffic  in  a sequence  of  functionally  related  sectors.  This  section 
considers  a transfer  function  model  which  links  the  time  series  of  air- 
craft loading  indices  in  a sector  with  that  in  the  next,  in  a single 
stream  of  traffic.  Box  and  Jenkins'  methods  for  Identification  and 
estimation  of  transfer  function  models,  discussed  in  Time  Series  Analys_i_s, 
Forecasting  and  Control , are  applied  to  the  modeling  of  responses  gen- 
erated by  a traffic  flow  from  a local  control  sector  to  a radar  departure 
control  sector,  both  handling  traffic  leaving  JFK  International  Airport. 


5.4.1  Transfer  Functions  Representing  Dynamic  Relationships 
in  a Sequence  of  Traffic  Flow 


The  dynamic  relationship  of  responses  in  two  adjacent  sectors  of 
a single  stream  of  traffic  can  be  represented  by  figure  5.7,  where  the 
"input"  series  X measures  the  response  in  the  first  sector,  and  the  "out- 
put" series  Y measures  the  response  in  the  second  sector.  The  "dynamic 
system"  is  the  mechanism  by  which  response  in  one  sector  is  translated 
into  response  in  the  next. 


While  the  traffic  stream  is  inherently  continuous  in  time, 
observations  are  generally  available  t)nly  at  discrete  intervals.  The 
relationship  between  Y and  X miglit  tims  be  approximated  by  a linear 
filter  of  the  form 


Y 

t 


^ "l^-l  "2^-2  ^ 


= (V. 


+ v^B  + 


V2« 


+ ...)X|.  = v(B)X^ 


(5.8) 
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Figure  5.8  Application  of  the  Impulse  Response  Fuhctlon 
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where  B is  the  backwards  operator  sucli  that  BX.  = X ..  Vhe  operator 
v(B)  is  called  the  "transfer  function"  of  the  filter,  and  tlie  Individual 
weiglits  Vq,  form  the  "impulse  response  function".  If  there  are 

b whole  periods  of  delay  before  a deviation  in  the  input  series  affects 
the  output,  then  v^,  v^,  ...,  v^  will  be  equal  to  zero.  Figure  5.8  shows 
how  the  Impulse  response  function  translates  the  input  series  X into  the 
output  series  ^ 

The  superposition  of  several  traffic  streams  entering  a given 
sector  could  be  represented  as  a sum  oi  several  transfer  functions.  For 
example,  if  the  output  sector  was  fed  by  two  traffic  streams,  then  the 
coiiibined  effect  Y*  would  be 

In  general,  an  output  series  can  be  represented  by 


Y*  = v^(B)X^_^  -f 


where  N is  noise  which  may  be  autocorrelated  and  is  usually  assumed  to  Ih 

independent  of  X, 

l.t 


5.4.2  Varameterization  of  the  Model 

Just  as  continuous  dynamic  systems  are  often  representable  by 
linear  differential  equations,  discrete  dynamic  systems  can  often  be 
represented  by  linear  difference  equations  of  the  form 

(1  + + ...+  C^V‘')Y^  = g(l  + T)^V  + ...  + 

where  V is  the  backward  difference  operator,  g is  the  steady-state  gain  ot 
the  system,  and  b is  the  number  of  time  units  of  delay  between  input  and 
output.  In  terms  of  tiie  backwards  operator  B=l-V,  Equation  (59)  can  be 
written  as 

(1  - A B - . . . - <S  b'^)Y  = ( aw  - o).B  - ...  - .i  ) B^X  (5.in) 

’ r t 0 J I 

Equation  (5.9)  or  (5.10)  defines  a transfer  function  model  of  order  (r,s). 
The  transfer  function  v(B)  is  given  by 

v(B)  = a"*  (B)  (-i(B)b'’  (5.11) 


where  A(B)  and  m(B)  are  poivnonials  in  B. 
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The  problem  of  fitting  the  transfer  function  model  Includes  deter 
mining  the  order  of  the  polynomials,  identifying  the  delay  parameter  b, 
estimating  the  parameters.  Box  and  Jenkins  have  suggested  methods  useful 
for  these  purposes,  and  the  corresponding  computer  programs  are  availabli' 
for  use.  In  the  next  subsection,  the  Box-Jenkins  techniques  are  applied  to 
the  analysis  of  responses  in  two  adjacent  ATC  sectors  through  which  a 
stream  of  traffic  is  passed. 

5. A. 3 An  Example  of  Modeling 

The  traffic  stream  selected  for  this  example  is  that  of  depart  in; 
aircraft  from  JFK  during  a 2 -hour  sample  period.  At  that  time,  two 
local  control  sectors  (LC)  were  in  operation,  one  of  which  was  handling 
solely  departures  and  fed  all  traffic  to  a radar  departure  controller  (I)P'  . 
The  data  represent  about  48  aircraft  in  the  2 -hour  period.  The  departure 
controller  also  handled  a small  amount  of  additional  traffic  from  other 
sources. 


The  model  appropriate  to  this  problem  is  thus 

= v(B)X^  + N|.  (5.12) 

where  shocks  are  structured  and 

= (f)“^(B)  0 (B)a^  (5.13) 

vhcrc  a^  instead  are  independent  identically  distributed  normal  variables 

t/ith  mean  zero  and  variance  equal  to  o'"  ( see  Box  and  Jenkins  for  refereare  . 

The  series  Y^  is  the  deviation  of  the  aircraft  loading  series  in  the  DP 
sector  about  its  mean  level,  is  a similar  quantity  for  the  LC  sector, 
and  Nj,  represents  deviations  due  to  additional  traffic  and  random  noise. 

Both  aircraft  loading  series  are  shown  in  figure  5.9.  The  following 
results  from  the  analysis  were  obtained: 

(1)  Both  the  input  and  output  series,  when  modeled  separately, 
were  well  fitted  by  second  order  autoregressive  models. 

Input : 

X = 0.867X  , - 0.176X  , + a X = 1.153 

t t-1  t-2  t 

Variance  of  original  scries  = 0.513 
Residual  mean  square  = 0.235 
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Figure  5.9  Aircraft  Loading  Tire  Series  for  N.Y.  I.C  Sector  513  and  DP  Sector  535 


Output : 


= 1.172Y^  ^ - 0.322Y^_2  + Y = 2.463 

Variance  of  original  series  = 1.803 
Residual  mean  square  = 0.366 

(2)  After  applying  the  prewhitening  transformation  estimated  from 
the  above  model  for  the  input  series,  significant  cross- 
correlation was  observed  between  the  prewhitened  input  and 
output  series  beginning  at  a lag  of  2 or  3 minutes.  The 
pattern  of  the  estimated  impulse  response  weights  suggested 

a model  of  order  (2,2)  or  lower;  tlie  exact  order  is  discussed 
below.  In  addition,  the  noise  remaining  in  the  output  series 
after  applying  the  estimated  transfer  function  was  v;ell  cii.  r 
acterized  by  an  AR  (2)  model. 

(3)  In  experimenting  with  various  models,  the  following  was 
observed : 

(a)  Models  with  r=2  gave  adequate  transfer  functions  for  b=2 
and  b=3.  Estimated  parameters  for  s>0  did  not  appear 
to  be  significant  in  either  case.  The  following  model 
appeared  to  be  appropriate: 

Y = P_-2,y X _ + ^ a 

(1-1.  58B+.715b"^)  ^ (1-.802B+.118B  ) ^ 

Residual  mean  square  = 0.263 
aiain  = 1.71) 

The  Impulse  response  function  for  this  model  has  weights 
equal  to  zero  at  lags  0 and  1 and  then  decays  as  a damped 
sine  wave,  eventually  satisfying 

V.  - 1.58u.  , + .715y.  = 0 j>2. 

J J-l  1-2 

(b)  Models  with  r=l  and  b=2  gave  good  results.  However,  with 

r=l  and  b='J,  tlie  estimation  of  the  parameters  in  the  non-linear 
least  squares  routine  was  liighly  unstable,  suggesting  tliai 
a probable  factorization  existed  to  a simpler  model  with  r=0. 

(c)  Models  with  r=0  gave  adequate  transfer  functions  for  s=3 

beginning  at  lag  3 and  for  s=4  at  lag  2.  These  models  have 
impulse  response  functions  witli  zero  values  at  lags  up  to 
b-1  and  s+1  values  V|,u]  , •••!  b+s  which  follow  no  fi:;..J 

pattern.  The  fitted  model  for  h=2  is. 
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= (.205+.389B+.  359B^+.513b\.384B^)X^_2 

^ L__.  a 

(l-.866B+.164B‘^) 

Residual  mean  square  = 0.286 
(Gain  = 1.85) 

(d)  The  parameter  estimates  for  the  AR(2)  noise  model  were 

relatively  insensitive  to  the  transfer  function  selected, 
at  least  for  those  models  which  appeared  to  be  adequate. 

On  the  basis  of  this  analysis,  it  appears  that  series  of  longer  length 
may  be  required  to  discriminate  among  the  various  possible  transfer  fuiu 
tion  models.  Analysis  of  traffic-  in  sectors  for  which  the  transit  time 
is  longer  might  also  help  to  clarify  the  appropriate  form  of  model.  Tin- 
main  advantages  of  models  with  r 0 over  those  with  r=0  are:  (i)  they  need 
fewer  parameters;  (il)  the  impulse  response  decays  smootlily,  while  the 
weights  in  the  models  with  r~0  appear  to  be  very  irregular.  Confirmation 
of  the  values  of  s selected  is  also  needed. 


5.5  Analysis  of  Nj|w_  York  Eij^ht-Hou r_  Data 

The  analysis  of  the  ATC  verbal  communications  system  r'^norted 
in  this  and  earlier  volumes  has  been  largely  based  on  a 2-hoiir  sample'  of 
101  air  traffic  control  positions  in  the  New  York  metroplex.  For  more 
thorough  studies  on  longer  periods  of  time,  additional  6 ho\jrs  of  data 
were  obtained  for  four  of  these  positions.  With  this  additional  data,  il 
is  possible  to  Investigate  more  closely  some  of  the  behavior  Inlierent 
to  the  input  stream  and  to  obtain  more  reliable  estimates  of  system 
parameters.  (For  the  purpose  of  illustration  only  the  data  and  analysis 
of  the  traffic,  flow  of  an  LC  sector  are  included  In  this  section.) 

TTie  analysis  of  data  over  the  8-hour  sample  periods,  as 
expected,  generates  new  problems  for  study.  A basic  assumption  in  all  of 
the  modeling  to  date  is  that  the  system  was  approximately  stationary 
Lliroughout  llie  period  of  observation.  This  assumption  is  not  tenable  tor 
the  longer  sample  period,  as  is  immediat c I \’  evident  from  a plot  of  the 
8-hour  aircraft  loading  series  for  N.Y.  sector  515LC  shown  in  figure  5.10. 
As  has  been  observed  bv  nan^'  researchers  in  the  transportation  field, 
transportation  systems  simply  do  not  remain  stationary  (even  approximatclv) 
for  more  than  a couple  of  hours  at  a time. 

From  visual  inspection  of  the  aircraft  loading  series  in  the-  four 
ATC  sectors,  there  appear  to  he  both  trends  and  step  shifts  in  the  Ic'vel 
of  the  process.  The  most  likely  explanation  for  such  nonstationary 
behavior  lies  in  the  arrival  and  departure  rates  of  aircraft  f r('m  l lie 
sectors.  In  particular,  the  assumption  of  homogeneous  Poisson  strc'.ims  ol 
aircraft  entering  and  leaving  the  sectors  is  extremely  dubious  for  the 
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entire  8-hour  periods.  Since  these  streams  are  the  fundamental  driving 
variables  in  the  system,  a caretui  analysis  of  their  structure  over  i ime 
is  of  primary  importance. 

Models  of  traffic  streams,  expressed  as  point  processes  have  been 
extensively  studied  over  the  last  years.  As  early  as  1936,  Adams 
(1936,  see  appendix  5-A  for  references),postulated  a Poisson  model  for 
the  distribution  of  counts  of  vehicles  on  a highway  passing  a point  in 
fixed  intervals  of  time.  Brelman  (1969)  formulated  a general  model  for  a 
traffic  stream  in  the  time-space  domain.  He  considered  the  paths  of 
vehicles  in  a stream  over  time  and  space  as  a "trajectory  process"  and 
defined  the  process  to  be  homogeneous  in  the  direction  of  a unit  vector  c 
if  its  stochastic  properties  were  invariant  under  translation  along  tlie 
vector. 


Now  since  most  observations  of  traffic  streams  are  performed  by 
collecting  data  at  a specific  point  in  space  or  at  a specific  point  in 
time,  streams  are  commonly  studied  as  point  processes  along  lines 
parallel  to  either  the  time  or  space  axis.  It  is  then  only  necessary  to 
assume  homogeneity  in  time  or  space,  depending  upon  the  perspective  of 
the  observer.  In  the  literature,  much  analysis  of  the  distribution  of 
vehicles  in  the  stream  is  based  on  assuming  that  the  observed  point 
process  is  stationary. 

In  the  case  of  streams  of  aircraft  througii  ATC  sectors,  the  con- 
cept of  trajectories  in  the  time-space  domain  does  not  have  as  precise  a 
physical  meaning  as  it  does  in  hirVway  traffic.  In  particular,  it  is 
necessary  to  adopt  the  meaning  of  "space"  so  tliat  all  trajectories  of 
aircraft  through  the  sector  are  considered  to  originate  at  a common  "point 
of  arrival"  and  to  terminate  at  a common  "point  of  departure."  The 
slope  of  individual  trajectories,  considered  to  be  straight  lines 
connecting  an  aircraft's  point  of  arrival  with  its  point  of  departur(>, 
represents  the  aircraft's  rate  of  passage  through  the  sect<5r  rather  than 
an  actual  physical  velocity. 

In  studying  each  8-hour  stream  of  aircraft,  two  point  proc- 
esses will  be  considered.  An  "arrival  stream"  will  be  considered  to  he 
a point  process  evolving  in  time  where  an  event  occurs  at  the  moment  an 
aircraft  initiates  the  first  communication  with  the  sector  controller. 

The  "departure  stream"  will  be  considered  to  be  a point  process  evolving 
in  time  where  an  event  occurs  at  the  moment  an  aircraft  concludes  its 
final  communication  with  the  sector  controller.  These  two  strc'ams  will 
have  the  same  stochastic  properties  if  tlie  underlving  trajectorv  process 
is  homogeneous  in  space. 

There  are  many  ways  of  representing  an  observed  point  proci'ss 
graphically.  Two  of  the  most  useful  are  shown  in  figures  111-  3.14 
for  the  arrival  and  departure  streams  observed  in  a local  control  sector. 
In  each  pair  of  figures  displayed  on  the  same  page,  ‘.'.i.  » oji  graph  plots 
the  cumulative  number  of  events  from  the  start  of  the  sample  period 
against  time.  The  bottom  graph  shows  counts  of  aircraft  arrivals  in  fixed 
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Figure  5.13  Plot  of  Cumulative  Departure  Count  versus 
Time  for  LC  Sector  515 
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Figure  5.14  Plot  of  Departure  Count  In  Twenty-Minute 
Intervals  for  LC  Sector  515 
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Intervals  of  20  minutes.  The  slope  of  the  cumulative  count  curve 

and  the  level  of  the  counts  in  the  bottom  graph  represent  the  rate  at 
which  aircraft  are  arriving  or  departing  from  the  sectors. 

Superimposed  on  the  cumulative  count  graphs  is  a dotted  line 
which  gives  the  theoretical  line  for  a homogeneous  process  with  an  expected 
number  of  counts  equal  to  the  total  number  observed.  Systematic  departuri' 
of  the  observed  counts  from  this  theoretical  curve  indicates  nonhomogeneity 
in  the  arrival  (departure)  rate.  An  investigation  of  the  graplis  displaycti 
is  being  conducted  to  determine  if  there  is  a statistically  significant 
change  in  rate  over  the  8— hour  periods  and  to  formulate  mathematical 
models  for  that  change. 

Another  way  of  studying  point  processes  is  to  consider  the  dis 
tance  in  time  or  space  between  successive  events.  Such  distances  are 
commonly  termed  "headways"  and  in  this  analysis  represent  the  times  between 
the  successive  arrivals  or  departures  of  aircraft  for  an  ATC  sector.  The 
consideration  of  headways  does  not  involve  the  arbitrary  selection  of  an 
interval  length  which  complicates  the  analysis  of  counts. 

Many  models  for  the  distributions  of  headways  have  been  proposed 
all  of  which  assume  that  the  underlying  point  process  is  stationary.  Edie 
(1974)  has  given  a general  summary  of  many  of  these  models,  which  includes 
renewal  models  and  semi-Markov  processes.  Taylor  et  al.  (1974)  l)ave 
recently  compared  several  Poisson  clustering  models,  which  are  a speci='l 
case  of  branching  renewal  processes  on  one  hand  and  doubly  stochastic 
Poisson  processes  on  the  other. 

In  studying  headways,  it  is  important  to  determine  whether 
successive  headways  are  independent.  Significant  serial  correlation 
can  sometimes  be  detected  by  visual  inspection  of  a plot  depicting  tlie 
location  of  paired  values  of  successive  headways.  Such  plots  for  the 
aircraft  traffic  streams  are  given  in  figures  5.15  and  5.17.  No  signifi- 
cant patterns  of  serial  correlation  are  easily  detectable,  but  a 
thorough  statistical  investigation  has  not  been  completed. 

An  important  function  related  to  the  distribution  of  headways 
is  the  survivor  function  R(x)  defined  by  R(x)=l-F(x),  where  F(x)  is  the 
distribution  function  of  tiie  headways.  For  a homogeneous  Poisson  proi- 
ess  F(x)  is  the  negative  exponential  distribution  and  log  R(x)  is 
linear  in  x.  Log  survivor  functions  which  are  convex  correspond  to 
k headway  distributions  which  are  underdispersed  with  respect  to  the  neg- 

ative exponential,  while  concave  log  survivor  functions  correspond  t ■ 
k overdispersed  headway  distributions.  The  graphs  in  figures  5.1b  and  lb 

^ show  the  empirical  log  survivor  functions  for  the  aircraft  arrival  aiui 

[•  departure  streams. 

II 

i It  is  very  important  to  note  that  the  analysis  of  alternative 

f models  for  headways,  and  for  other  facets  of  the  ATC  requires  that  tlu- 

f underlying  process  be  stationary,  at  least  to  tlie  second  order.  It  is 

■ thus  important  either  to  be  able  to  isolate  periods  in  wliich  the  trattlc 
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Figure  5.15  Scatter  Plot  of  Pairs  of  Consecutive 

Headways  In  LC  Sector  515  Arrival  Stream 


Cj 


Figure  5.16  Log  Survivor  Function  for  Headways 
In  LC  Sector  515  Arrival  Stream 
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intervals  of  20  minutes.  The  slope  of  the  cumulative  count  curve 
and  the  level  of  the  counts  in  the  bottom  graph  represent  the  rate  at 
which  aircraft  are  arriving  or  departing  from  the  sectors. 


Superimposed  on  the  cumulative  count  graphs  is  a dotted  1 ine 
which  gives  the  theoretical  line  for  a homogeneous  process  with  an  expected 
number  of  counts  equal  to  the  total  number  observed.  Systematic  departuri 
of  the  observed  counts  from  this  theoretical  curve  indicates  nonhomogene i tv 
in  the  arrival  (departure)  rate.  An  investigation  of  tlie  graphs  displa’d 
is  being  conducted  to  determine  if  there  is  a statistically  signifitan' 
change  in  rate  over  the  8— hour  periods  and  to  formulate  mathematical 
models  for  that  change. 

Another  way  of  studying  point  processes  is  to  consider  the  dis 
tance  in  time  or  space  between  successive  events.  Sucli  distanc-s  are 
commonly  termed  "headways"  and  in  this  analysis  represent  the  times  between 
the  successive  arrivals  or  departures  of  aircraft  for  an  ATC  sector.  Tlie 
consideration  of  headways  does  not  involve  the  arbitrary  selection  of  an 
interval  length  which  complicates  the  analysis  of  counts. 

Many  models  for  the  distributions  of  headways  have  been  |)ruposed 
all  of  which  assume  that  the  underlying  point  process  is  stationary.  Kdie 
(1974)  has  given  a g-^neral  summary  of  many  of  these  models,  whicli  ini  ludis 
renewal  models  and  semi-Markov  processes.  Taylor  et  al.  (1974)  have 
recently  compared  several  Poisson  clustering  models,  which  are  a spec'-’l 
case  of  branching  renewal  processes  on  one  hand  and  doubly  stochastic 
Poisson  processes  on  the  other. 

In  studying  headways,  it  is  important  to  determine  whether 
succf'ssive  headways  are  independent.  Significant  serial  correlation 
can  sometimes  be  detected  by  visual  inspection  of  a plot  depicting  tlie 
location  of  paired  values  of  successive  headways.  Sucii  plots  for  the 
aircraft  traffic  streams  are  given  in  figures  5.15  and  5.17.  No  sign!  M- 
cant  patterns  of  serial  correlation  are  easily  detectable,  but  a 
thorough  statistical  investigation  has  not  been  completed. 

An  important  function  related  to  the  distribution  of  lu'adways 
is  the  survivor  function  R(x)  defined  by  R(x)=l-F(x),  where  F(x)  is  tlu 
distribution  function  of  tlie  headways.  For  a homogeneous  Poisson  proc- 
ess F(x)  is  the  negative  exponential  distribution  and  log  R(x)  is 
linear  in  x.  Log  survivor  functions  which  are  convex  correspond  to 
headway  distributions  which  are  underdispersed  with  respect  to  the  neg- 
ative exponential,  while  concave  log  survivor  functions  correspond  to 
overdispersed  headway  distributions.  The  graphs  in  figures  5.1b  and  5.18 
show  the  empirical  log  survivor  functions  for  the  aircraft  arrival  and 
departure  streams. 

It  is  very  important  to  note  that  the  analysis  of  alti’rnative 
models  for  headways,  and  for  other  facets  of  the  ATC  requires  that  the 
underlying  process  be  stationary,  at  least  to  the  second  order.  It  is 
thus  imoortant  either  to  be  able  to  isolate  periods  in  which  the  traltii 
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Figure  5.15  Scatter  Plot  of  Pairs  of  Consecutive 

Headways  In  LC  Sector  515  Arrival  Stream 
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Figure  5.17 
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Log  Survivor  Function  for 
Headways  In  LC  Sector  515 
Departure  Stream 

streams  are  approximately  stationary  or  to  manipulate  the  point  process 
in  such  a way  as  to  achieve  stationarity  in  a transformed  metric.  Current 
research  is  being  conducted  in  both  of  these  areas.  Results  from  the 
studies  of  the  8 -hour  sample  data  can  provide  much  valuable  infor- 
mation to  the  operation  of  ATC  simulations  for  a longer  period  of  time  in 
the  future. 


5.6  (Queuing  Time  Analysis  and  Forecasting  for  ATC  Communications 

In  section  3.2  we  have  pointed  out  some  applications,  as  well  as 
the  limitations  of  the  usefulness,  of  the  conventional  queuing  theory  to 
the  analysis  of  ATC  communications.  It  was  found  that  the  expected  queuing 
time  per  communication  transaction  derived  from  the  theoretical  results 
can  only  provide  an  approximate  lower  bound  to  the  value  realized  in  an 
ATC  system.  Although  information  obtained  from  analytic  solutions  is 
helpful  for  selecting  the  range  for  simulation  experimentation,  it  is  not 
adequate  for  providing  final  estimations. 

To  pursue  further  the  analysis  of  the  ATC  queuing,  and  to  perfo'm 
forecasting  which  may  be  useful  for  network  simulations  in  the  future,  in 
this  section  a transformation  technique  useful  for  modelling  queuing  timi 
processes  is  first  Introduced,  Applications  of  this  technique  are  then 
illustrated  for  two  simulated  ATC  communications  responseso 


A portion  of  the  material  presented  in  this  section  has  been 
described  in  more  detail  in  an  article  by  Hsu  (D.  -A.  Hsu,  "Queuing 
Time  Forecasting  Civil  Engineering  Research  Report  //75-TR-7,  Princetoi, 
University,  July  1975).  Readers  are  referred  to  tliis  article  for  tech- 
nical references  and  discussions. 


5.6.1  Times  in  System  and  a Transformation  Technique 

Although  the  series  of  queuing  times,  recorded  in  the  order 
of  arrivals,  of  each  communication  transaction,  is  an  important  variable 
for  study,  some  difficulties  are  inherent  to  the  analysis  of  this  variable. 
One  of  the  causes  is  that  the  queuing  time  variable  has  a nonzero  prob- 
ability at  the  point  of  zero.  This  makes  the  distribution  function 
difficult  to  express  and  unsuitable  for  analysis.  In  order  to  avoid  tills 
problem,  we  select  another  variable  wliich  is  closely  related  to  the 
variable  in  question  for  analysis.  That  is  the  process  of  times  in 
system  of  each  communication  transaction,  defined  as  the  sum  of  the  queuing 
time  and  the  CT  length  associated  with  each  communication  transaction. 

Let  this  series  be  denoted  by  i W* 1 , where  t indicates  the  arriving  order 
of  the  CT's, 

A preliminary  analysis  of  general  processes  of  times  in  system 
suggests  that  a first  order  autoregressive  model  AR(1)  is  a natural  process 
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for  an  independent  homogeneous  arrival  system,  given  certain  transtormations 
being  made.  However,  as  has  been  pointed  out  in  an  earlier  analysis  (see 
subsection  3.2.1),  the  CT  arrivals,  in  general,  do  not  conform  to  the 
independency  requirement  Imposed.  Therefore,  a general  class  of  models, 
such  as  those  of  autoregressive  or  moving  average  or  the  mixture  of  botli, 
suggested  by  Box  and  Jenkins,  should  be  considered  in  the  modelling 
of  the  {W*}  series  just  mentioned.  Nevertheless,  in  a recursive  trans- 
formation''procedure,  an  AR(1)  model  is  still  a good  starting  model.  The 
model  can  be  refined  after  the  transformation  having  reached  an  intermediate 
stage.  In  mathematical  form,  an  AR(1)  model  is  expressed  as 

(W*  - p ) - - u)  = (5.1/) 

where  p is  the  overall  mean  of  the  series  W*  and  a is  an  independent, 
identically  distributed  normal  variable. 

A transformation  procedure,  useful  for  the  analysis  of  the  series 
of  {W*}  generated  from  the  ATC  simulation  model,  consists  of  the  following 
stages^ 

(a)  Fit  the  time  series  model  specified  by  Equation  (5.14),  using 
the  least  square  estimator  for  ij>,  to  the  data  generated  from  a simulation 
experiment  (the  simulation  run  provides  a sufficiently  large  number  of 
observations,  with  a given  set  of  input  parameters^  and  compute  the 
residuals  using  the  following  formula. 


= (Wj.  - p) 


- - k) 


t=2,  3, 


‘:5.15) 


where 


N 

= 

t=i 

N 

y 

= t=2 


W*/N 


(Wj.  - .)(W|._^  - p) 


N-1 

I 

t=l 


* ■'2 
(W^  - p)^ 


and  N = the  number  of  observations  contained  in  tlic  series. 

(b)  Separate  the  set  of  data  {W^l  into,  say,  10  equal-size  groups 
in  order  of  their  magnitude  (i.e.  divide  the  data  at  the  1 x 10%  quantiles 
for  i = 1,2,  .,.,10)  and  label  these  observations  by  the  ordinal  number  of 
the  group  to  which  the  observations  belong.  Then  classify  the  a *s  into 
the  groups  according  to  the  label  assigned  to  W*.  (For  instance^  if  W*^ 
is  labeled  "5",  the  a^^  should  be  assigned  to  the  fifth  group.)  Then 
compute  the  following: 
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(5.16) 


(5.17) 


where  1.  contains  the  elements  which  belong  to  the  i group.  Let  us  further 
define  that 


R.  = n /u 
1 i a 


(5.18) 


Then,  denoting  the  transformed  values  as  W^,  the  following  transformation 
is  suggested: 


"t  ■ \'h 


for  tr,  I, 


(5.19) 


- * , * . 
= max  {W^-i  + - max 


t'el.  , 
1-1 


i = 2,3,... 


(c)  Again  fit  the  first  order  autoregressive  model  to  the  series 
of  Wj.  and  obtain  a new  set  of  residuals.  The  distribution  of  residuals  is 
then  examined  using  the  conventional  measurements  of  skewness  and  kurtosis, 

~ ^ where  = 1,2, 3, 4,  are  the  i*^*’  central 

moments  (with  respect  to  the  mean)  computed  from  the  data  observed.  If  the 
value  of  is  significantly  different  from  xero  (tables  and  graphs  compiled 
for  the  tests  of  and  are  available  in  Snedecor  and  Cochran,  Statistical 
Metl^d,  1970i  the  procedure  specified  by  F.qtiations  (5.16)  - ( 5. 1 9)  a'lid'Tli’en' ' 
again  (5.15)  using  the  last  set  of  transformed  data  is  suggested  to  he 
Recursively  operated  until  the  value  of  is  brought  to  a satisfactory  level 

To  examine  the  adequacy  of  the  expression  specified  in  Equation 
(5.14),  the  methods  for  diagnostic  checking  suggested  by  Box  and  .lenkins  in 
time  series  modeling  can  be  employed  to  investigate  the  nonexistence  of  the 
serial  correlation  in  the  final  residuals,  (In  the  case  where  a is  far 
away  from  normality  these  methods  should  be  used  with  caution  Further 
study  should  be  performed  to  investigate  the  influence  ol  nonnormality  upott 
the  testing  statistics.)  The  test  statistics  available  are  (i)  ^ ~ 

r (a  )/  I'n,  where  r (a  ) is  the  sample  autocorrelation  function  of  the 
fina^  Residuals  at  ^ag  k and  n is  the  relevant  sample  size;  and  (ii) 

Q = n ) r ^(a  ),  where  K is  a constant  which  can  he  adequately  deter- 

k=l  ''  ' 

mined  for  each  specific  case.  The  first  test  statistic  useful  for 

testing  individual  autocorrelation  coefficients  is  normally  (i  i si  r ibut  «'d 
with  mean  zero  and  variance  one,  if  tlie  coefficient  is  virtually  zero.  In 


2 

r^^  (a^),  where 


is  a constant  which  can  he  adequately  deter- 
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addition,  the  second  statistic  Q is  distributed,  in  our  case  of  first  order 
autoregressive  process,  as  a (K-1)  variable,  where  K-1  represents  the 
degree  of  freedom,  and  useful  as  a portmanteau  test  of  all  relevant 
autocorrelation  coefficients. 

Model  modification  can  then  be  proceeded  using  Box-denkins 
technifiues.  Diagnostic  instruments  described  in  the  last  paragraph  can 
again  be  used  to  provide  information  helpful  for  identifying  the  lack  of 
fit  ana  suggesting  the  direction  of  modifications. 


5.6.2  Examples  from  ATC  Simulations 


Two  examples  are  provided  in  this  subsection  to  illustrate  the 
modeling  procedure  described  earlier  in  this  section.  Extensions  of  tlie 
resulting  models  to  the  queuing  time  forecasting  are  finally  mentioned. 


^x^^l^e  j:  The  series  of  times  in  system  for  each  of  tlie 

communication  transactions,  denoted  by  fW* 1 was  collected  from  a simulation 
exercise  of  2-hour  period  for  New  York  .Sector  524  LC  (ABE  airport)  at  the 
historical  traffic  density  (0  = 25  aircraft  arrivals/hour).  Totally  210 
observations  were  used  in  the  following  modeling  procedure. 


Stage  (i):  The  entire  series  of  fW  1 were  separated  into  10 

groups,  each  with  size  23,  according  to  their  magnitudes  in  seconds.  The 
10  partition  points  are:  (4.69,  7.53,  9.53,  12.01,  14.36,  17.03,  21.06, 
27.94,  43.23,  98.19),  where  the  last  value  is  the  maximum  value  observed 
from  the  sample.  An  AR(1^  model  was  then  fitted  to  this  data.  Results  of 
fitting  are:  = 0.8453.  tV  = 96.  1309,  u =^20.0514,  x^(a  ) = 0.  3766  and 

^2^^t^~  1*6135.  Both  ^ 1^2^  a ) s ign i f icanF 1 y depart  from 

normality,  A transformation  is  therefore  needed. 


Stage  (ii):  From  Equations  (5.16)  - (5.18),  values  of  R.  wi-re 

calculated  and  found:  (R  , R , = (0  .8327,  1 . 1‘)78,  1 , 4 508 , ' 1 . 265  7 , 

1.7133,  1.6064,  1.4860,  0,8126,  0.8446,  0.6058).  Tliese  values  indicat^ 
that  the  transformation  function  is  convex  at  the  lower  range  of  the  W 
value  and  is  concave  at  the  higher  range.  Following  tin*  transformation 
specified  by  Equation  (5.19),  we  transformed  {17  1 int('  {W  ' and  again 
fitted  the  AR(1)  model  into  the  new  set  of  data.  It  was  K>und  that  tlie 
new  results  are:  = 0.7595,  o^  = 110.4477,  n = 21.  2753,  >j(a  ) = 0.0623 

and  X ,(a^)  = -0.0617.  This  time  bntii  Y-[(a^)  and  r^(a^)  conform  to  the 
conditions  for  normality. 


Stage  (lii):  Further  examination  on  the  sample  autocorrelation 

function  calculated  from  the  residuals  indicated  that  the  value  of  r*(a,) 
is  equal  to  3.34  and  is  significantly  different  from  zero  at  the  0,01  level. 
In  addition  the  value  of  Q,  based  on  the  first  24  autocorrelation  co- 

24 

efficients,  is  Q(a  ) = 230  ) r (a  ) = 40.23  being  gre.iter  than  tin-  .025 

'■  k-i  " ^ 

significance  level  of  a x^  vari.abli-  with  23  di'grees  ot  freedom.  A mixed 
autoregressive  and  moving  average  model,  with  the  moving  average  parameter 
being  at  tlie  seventh  order,  was  then  trii'd  on  I lu'  data  of  . I'lie  ni'W 
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estimates  of  the  parameters  are;  4i  = 0.7400,  G = 0.1847,  = 107.4800, 

u = 21.0410,  Yj^(a  ) = 0.0511  and  y ) = -0.0655.  The  estimate  0,  is 
significantly  different  from  zero  at^'the  0.01  level.  A further  check  on 
the  residual  autocorrelation  function  indicated  that  the  value  of  r*(a  ) 
has  been  reduced  to  0.61  and  the  value  of  fj  has  been  brovight  down  to  2*5.81 
which  is  well  below  a critical  point.  Thp  transformation  function  and  the 
last  set  of  parameter  estimates  are  the  essential  information  for  con- 
structing a predictive  proced\ire  for  qtieulng  times,  as  will  be  pointed  out  in 
the  concluding  remarks  of  this  section. 

Example  2 : A similar  series  of  W was  also  collected  from  a 

simulation  run  for  New  York  sector  553  AR  (La  Guardia  airport)  at  a traffic 
density  equal  to  21  aircraft  arrivals  per  tiour.  Results  of  modeling  on  the 
total  of  370  observations  are  briefly  reported  below: 

Stage  (i):  The  10  partition  points  are  (6.43,  9.96,  13.39,  16.26, 

19.81,  23.32,  27.69,^35.76,  45.11,  131.61).  The  first  fitting  of  the  AR(1) 
model  gave  results:  <j>  = 0.7590,  = 159.3731,  p = 24.4667,  y-^  (a^ ) = - 0. 0942 , 

and  Y?(a  ) = 2.4327.  The  value  of  y^(a  ) is  much  too  hlgli  to  be  from  a 
normal  distribution.  Transformations  were  therefore  pursued. 

Stage  (ii):  Estimates  of  R.  are:  (R,,  R„,  ...,  Rf())  = (0.6923, 
1.2136,  1.4175,  1.4746,  1.2576,  1.44k,  1.28lk  173633,  0.9306,  0-  5921), 

The  new  set  of  estimates  of  parameters  in  the  AR(1)  model,  for  the  trans- 
formed data,  are  : $ = 0.6984,  o“  = 185.3221,  n = 25.9928,  y (a. )=  -0.1811 
- - 1 ^ 
and  ~ 0.  3540.  The  values  of  been  substantially  rediu  ed 

and  is  now  well  below  ,i  critical  point. 

Stage  (^ii):  Checks  on  the  residual  autocorrelation  function 
showed  that  all  r (a  )'s  for  k = 1,  2,  ...  24  are  within  -2.0  and  2.0,  and 
0 = 15.58.  No  model  modification  is  proved  necessary  from  these  evidences. 


5.6.3  Concluding  Remarks  on  Queuing  Time  Analysis 

Given  the  parameter  values  of  an  autoregressive  model  or  a moving 
average  model  or  the  mixture  of  botii,  tlie  forecasting  of  variables  generated 
from  the  process  is  a st raightfciward  matter.  Confidence  limits  for  fore- 
casts obtained  from  Box-.Ienkins  procedure  can  be  transformed  hack  to  an 
original  scale  using  the  inverse  transformation  function  determined  by  R . ‘ , 
Applications  of  the  procedure  illustrated  above  to  queuing  network  simu- 
lations are  natural  extensions  of  this  work. 
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CHAPTER  6 


CONCLUSIONS  AND  RECOMMENDATIONS 


Discussions  on  many  subjects  of  interest  concerning  the  properties 
and  the  applications  of  the  simulation  model  of  Air  Traffic  Control 
Communications  have  been  provided  in  detail  in  the  previous  chapters. 

However,  conclusions  scattered  over  individual  chapters  or  sections  need  to 
be  gathered  to  form  a complete  piece  of  picture  which  represents  our  effort 
devoted  in  this  volume.  Further,  recommendations  adequately  drawn  from  the 
research  results  are  also  helpful  in  suggesting  future  pursuit  on  the  subject 
of  ATC  communications.  For  these  reasons  we  present  in  the  following 
sections  the  overview  and  the  outlined  recommendations. 


6 . 1 Conclusions 

In  this  volume  we  have  demonstrated  tlie  general  applicability  of 
the  simulation  model  of  ATC  communications  in  both  distinct  geographical 
areas  and  different  types  of  studies.  We  have  compared  the  communications 
data  of  New  York  versus  Houston  ARTCC  and  found  that  only  one  or  two  of  the 
Sensitive  input  variables,  most  likely  the  number  of  communication  trans- 
actions per  aircraft,  need  to  be  adjusted  in  order  to  perform  communications 
simulations  in  different  geographical  areas.  Further,  we  have  illustrated 
the  competence  of  the  model  in  answering  many  questions  of  practical  interest. 
The  questions  Investigated  Include  the  evaluations  of:  (1)  the  communication 
capacity  of  a typical  sector  performing  a particular  function;  (ii)  the 
effects  of  a tone-burst  scheme  on  delaying  communications  performance;  (iii) 
the  influences  of  reducing  the  number  of  transmissions  per  communication 
transaction  upon  the  communication  responses.  On  the  first  subject,  we  havi- 
determined  the  communication  capacities  of  nine  sector  functions  in  the  New 
York  area  through  the  simulation  experiments  and  a special  statistical  tech- 
nique in  detecting  the  onset  of  unstable  performance.  For  tlie  second  problem 
we  have  evaluated  the  reduction  of  communication  capacity  and  the  induced 
extra  communications  delay  at  various  combined  values  of  the  tone  duration 
and  the  expected  number  of  aircraft  arrivals.  Finally,  for  tlie  third  question, 
we  are  able  to  provide  numerical  answers  to  the  relative  changes  of  the 
communications  responses  at  various  rates  of  reduction  on  the  number  of 
transmissions  per  communication  transaction. 

In  addition  to  the  researcli  effort  devoted  to  t Ite  studies  of  the 
communications  in  a single  sector,  a netVvtork  formulation  of  c-nroute  control 
sectors  has  been  proposed  and  a second  order  Markov  chain  model  is  used  to 
describe  the  transferring  of  an  aircraft  from  one  sector  to  the  next  for  the 
New  York  ARTCC.  Estimations  of  the  aircraft  arrival  rates,  the  transition 
patterns,  and  the  flow  rates  have  been  studied  and  empirically  illustrated 
using  the  New  York  data.  This  network  modelling  eflort  has  prepared  lor  the 
simvilation  of  ATC  communications  of  a netw(>rk  as  a whole. 
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Othor  subiects  also  explored  in  this  volume  are:  (1)  a sensitivitv 
analysis  of  the  expected  ATC  communications  channel  utilization  with  resjunt 
to  various  input  variables;  (ii)  the  construction  of  operational  indices  for 
ATC  communications  performance  that  simplify  the  presentation  of  various 
responses;  (iii)  a transfer  function  model  for  aircraft  flows  at  adiaeent 
ATC  sectors  describing  the  dynamic  relationship  betv^^een  the  aircraft  loadings 
in  tv;o  consecutive  sectors;  (iv)  an  analvsis  of  New  York  air  tralfic  dat;i 
of  an  H-hour  period  in  which  the  traffic  density  is  apparently  nonhomogeneous ; 
and  (v)  a queuing-time  analysis  and  ft)recasting  for  ATC  communications  with 
potential  applications  to  the  simulation  of  network  communications. 


6.2  Recommendations 


As  logical  results  of  the  conclusions  presented 
recommendations  concerning  potential  applications  of  the 
future  investigations  of  ATC  communications  are  provided 


above,  several 
simulation  model  in 
as  the  following: 


(1)  As  has  been  pointed  out  earlier  in  section  5.3  of  volume  ITT 
and  section  4.1  of  this  volume,  to  Investigate  the  ATC  communications 
performance  in  a metropcjlltan  area  such  as  Nevvc  York  as  a whole,  which 
consists  of  a network  of  control  sectors,  is  a natural  extension  of  our 
effort  to  date.  Simulations  of  a sequence  of  sectors,  each  simulation  being 
executed  using  the  constructed  model,  can  be  accomplished  in  accordance  with 
the  geographical  or  functional  adjacencies  of  the  sectors.  In  othor  words, 
the  aircraft  departure  stream  of  one  sector  is  to  be  used  as  inputs  to  the 

next  sector  (or  several  sectors)  v^hlch  is  (are)  functionally  related.  Sucli  a 

network  study  would  be  useful  in  the  evaluation  of  the  present  performance  of 

a facility  and  the  prediction  of  performance  in  hypothetical  environments. 

Thus,  experimental  results  generated  from  network  simulations  could  serve  as 
valuable  aids  in  decision  making  processes  concerned  with  changes  in  tlie  rules 
of  communication  operations,  the  Introduction  of  new  communlcat ion  technimies, 
expansions  or  constructions  of  airports,  etc. 


(li)  Air  Traffic  Control  performance  in  distinct  geographical  areas 
may  be  different  in  one  aspect  or  another,  as  has  been  shown  in  a comparison 
of  New  York  versus  Houston  data  presented  in  chapter  2 of  this  volume. 
UTiether  these  differences  warrant  discriminatory  requirements  on  the  AT(’ 
performance  is  a question  requiring  further  studies.  Here  we  recommend  tliat 
the  simulation  model  can  help  answer  the  question  empirically  and  comprehen- 
sibly. 


(iii)  As  has  been  illustrated  in  section  3.4  of  tins  volume,  tlie 
simulation  model  provides  certain  insights  about  the  effectiveness  of  a change 
in  a performance  parameter  upon  various  responses  of  concern.  This  tvpe  of 
study,  to  a certain  extent,  provides  a basis  to  the  cost  and  benefit  analvsis 
of  an  improvement  measure  and  may  serve  as  a concrete  motive  to  certain  govern- 
ment actions. 


(iv)  Questions  concerning  the  effect  of  access  control  techniques, 
channel  allocation  methods,  etc.,  upon  air/ground  communications  are  being 
raised  by  those  engaged  in  the  study  of  satellite  mobile  communications/ 
surveillance  systems  Investigations  of  the  efficiency,  flexlbilltv,  and 
capacity  of  proposed  alternatives  can  be  supported  by  simulation  experiments 
using  methodology  developed  in  the  present  work. 


It  is  the 

gations  presented  in 
extensions  can  serve 
associated  with  many 


belief  of  those  who  have  participated  in  the  invest  i- 
the  four  volumes  that  the  simulation  model  and  related 
as  a very  powerful  tool  in  analyzing  various  outcomes 
future  governmental  decisions  of  concern. 
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